[R-sig-ME] NLME model syntax

nrm2010 nrm2010 at zoho.com
Thu Feb 24 23:16:18 CET 2011


Hello, 

I'm trying to use NLME for a non-linear model with one random effect 
that is the same as the one grouping level (chamber). Using the following 
statement: 

test <- nlme(flux ~ b0 * exp((b1*soiltemp) - (b2*soiltemp^2)) 
* vpd^b3 * mstsoil2^b4 * airtemp^b5, data = ac5, 
fixed = list(b0 + b1 + b2 + b3 + b4 + b5 ~ 1), 
random = chamber ~ 1, 
start = c(b0 = 570, b1 = 0.08, b2 = 0.002, 
b3 = -0.06, b4 = 3, b5 = 0.3), subset = airtemp > 0) 

I always get the error message 
"Error in eval(expr, envir, enclos) : object 'chamber' not found" 

Data "ac5" is a groupedData object with chamber as the grouping factor: 

> head(ac5) 
Grouped Data: flux ~ soiltemp + vpd + mstsoil2 + airtemp | chamber 

I've tried many variations such as 
random = list(chamber ~ 1)  or random = ~1|chamber
or adding the "groups" statement or trying to put chamber on the 
RHS of the model statement, after the fixed effects, 
but I always get the same error message. 

What is the correct way to formulate the call to nlme? Or is this not a syntax error? 

Thank you in advance. 

Toby Gass 

tobygass at cnr dot colostate dot edu




More information about the R-sig-mixed-models mailing list