t <- seq(0, 30, by = 6)
k <- 0.3
n <- 40
set.seed(100)
y <- FOMTm(t, k, n) * (1 + rnorm(length(t), 0, 0.05))
nlsFOMT <- nls(y ~ FOMTm(t, k, n),
data = list(y = y, t = t),
start = par_est_FOMT(t, y)
)
summary(nlsFOMT)
Run the code above in your browser using DataLab