[R-sig-ME] nlme vs nlmer
Christof Kluß
ckluss at email.uni-kiel.de
Fri Sep 28 17:53:41 CEST 2012
Hi
can you give me a hint how to convert a nlme(..) like
nlme(Y ~ fun(X,a,b,c),
fixed = list(a ~ 1, b ~ 1, c ~ 1),
random = list(a ~ 1, b ~ 1, c ~ 1),
groups = ~ location,
data = measurements,
start = c(a = 350, b = 0.2, c = 120))
to an equivalent nlmer(...)? The following seems to be wrong
nlmer(Y ~ fun(X,a,b,c) ~ (a + b + c | location),
data = measurements,
start = c(a = 350, b = 0.2, c = 120))
it throws "gradient attribute of evaluated model must be a numeric matrix"
thx
Christof
More information about the R-sig-mixed-models
mailing list