Displaying 1 result from an estimated 1 matches for "fm2co2".
Did you mean:
fm1co2
2006 Jul 18
2
Using corStruct in nlme
...ach 'Plant'. I have replicated the problem using the CO2 data in nlme (based off of the Ch08.R script).
This works (when 'conc' is the same for each 'Plant':
(fm1CO2.lis <- nlsList(SSasympOff, CO2))
(fm1CO2.nlme <- nlme(fm1CO2.lis, control = list(tolerance = 1e-2)))
(fm2CO2.nlme <- update(fm1CO2.nlme, random = Asym + lrc ~ 1))
CO2.nlme.var <- update(fm2CO2.nlme,
fixed = list(Asym ~ Type * Treatment, lrc + c0 ~ 1),
start = c(32.412, 0, 0, 0, -4.5603, 49.344),
weights=varConstPower(fixed=list(const=0.1, power=1)), verbose=T)
CO2.nlme.CAR<-update(CO2.nlme.v...