add replicate.R