reps
data points from this at timepoints obs
simulategp.perturbed(N = 100, reps = 3, obs = c(0, 5, 10, 15, 20), gps.ctrl = NULL, filename = NULL, xs = seq(0, 20, 0.2), l.noise = 12, l.shape = 2, l.scale = 2.5, sigmaf.noise = 0.25, sigmaf.shape = 4, sigmaf.scale = 0.5)
simulategp