[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