[R-sig-ME] including covariance pattern in a nlme model
ONKELINX, Thierry
Thierry.ONKELINX at inbo.be
Thu Jun 30 10:46:24 CEST 2011
I think you need some extra brackets.
update(LVB.nlme, corr=corAR1())
> -----Oorspronkelijk bericht-----
> Van: r-sig-mixed-models-bounces at r-project.org [mailto:r-sig-mixed-models-
> bounces at r-project.org] Namens Mike Gibson
> Verzonden: woensdag 29 juni 2011 18:01
> Aan: r-sig-mixed-models at r-project.org
> Onderwerp: [R-sig-ME] including covariance pattern in a nlme model
>
>
> I am trying to include a defined covariance pattern in my nlme model. I followed
> an example made by Bates online. I am close but I keep getting an error
> message. Any help would be greatly appreciated. The details are below.
>
>
> Here is my code: I already defined an equation for LVB (von Bert growth
> equation).
>
> LVB.nlme=nlme(L~LVB(Age,t0,Lmax,K), data=datagr,
> fixed=list(t0~1,Lmax~1,K~1), random=t0+Lmax+K~1, start=list(fixed=c(t0=-1.0,
> Lmax=700, K=0.25)), control=nlmeControl(maxIter=10000, tolerance=1e-2,
> msVerbose=TRUE))
>
> #Everthing is fine up to this point. Next, I define the covariance structure.
>
> LVB2.nlme<-update(LVB.nlme, corr=corAR1)
>
> #This is where I am stuck. I get the following error message:
>
> ERROR: object of type 'closure' is not subsettable
>
>
> I would greatly appreciate any help.
>
> Mike
>
>
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
More information about the R-sig-mixed-models
mailing list