Displaying 2 results from an estimated 2 matches for "fm1co2".
2006 Jul 18
2
Using corStruct in nlme
...e in the nlme package. The one major
difference is that in my case the 'conc' steps are not the same for each '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.560...
2005 Nov 02
1
NLME
...odels. R package version 3.1-65.
on a WINDOWS 2000 machine
I am trying to run the scripts from the Mixed Models book and am running into some difficulty - could any one tell me why this doesn't work (from pages 369 - 370)
CO2 #FINE
plot(CO2,outer= ~Treatment*Type,layout=c(4,1)) #FINE
fm1CO2.lis <- nlsList( SSasympOff, CO2 ) #FINE
fm1CO2.lis #FINE
fm1CO2.nlme <- nlme( fm1CO2.lis ) #COPIED DIRECTLY FROM BOOK!!??
Gives the error message:-
Error in nlme.formula(model = uptake ~ SSasympOff(conc, Asym, lrc, c0), :
Maximum number of iterations reached without convergen...