
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]]


