add stopifnot() thresholds to each of the replication vignettes

[ajdamico]

something like stopifnot( abs( estimate - actual ) / estimate < 0.01 ) to break the build if the numbers ever fall outside of some range