data(Loblolly)
fm1 <- nlme(height ~ SSasymp(age, Asym, R0, lrc),
data = Loblolly,
fixed = Asym + R0 + lrc ~ 1,
random = Asym ~ 1,
start = c(Asym = 103, R0 = -8.5, lrc = -3.3))
fm2 <- update(fm1, random = pdDiag(Asym + lrc ~ 1))
summary(fm1)
summary(fm2)
Run the code above in your browser using DataLab