[R] nlme model specification
Martin Henry H. Stevens
hstevens at muohio.edu
Wed May 17 22:51:55 CEST 2006
Hi folks,
I am tearing my hair out on this one.
I am using an example from Pinheiro and Bates.
### this works
data(Orange)
mod.lis <- nlsList(circumference ~ SSlogis(age, Asymp, xmid, scal),
data=Orange )
### This works
mod <- nlme(circumference ~ SSlogis(age, Asymp, xmid, scal),
data=Orange,
fixed = Asymp + xmid + scal ~ 1,
start = fixef(mod.lis) )
### I try a slightly different model specification for the fixed
effects, and it does not work.
### fixed = list(Asymp ~ 1, xmid ~ 1, scal ~ 1)
### I tried following the example on page 355.
mod <- nlme(circumference ~ SSlogis(age, Asymp, xmid, scal),
data=Orange,
fixed = list(Asymp ~ 1, xmid ~ 1, scal ~ 1),
start = fixef(mod.lis) )
### I get
Error in parse(file, n, text, prompt) : syntax error in "~ "
nlme version 3.1-71
> version
_
platform powerpc-apple-darwin8.6.0
arch powerpc
os darwin8.6.0
system powerpc, darwin8.6.0
status
major 2
minor 3.0
year 2006
month 04
day 24
svn rev 37909
language R
version.string Version 2.3.0 (2006-04-24)
>
Dr. M. Hank H. Stevens, Assistant Professor
338 Pearson Hall
Botany Department
Miami University
Oxford, OH 45056
Office: (513) 529-4206
Lab: (513) 529-4262
FAX: (513) 529-4243
http://www.cas.muohio.edu/~stevenmh/
http://www.muohio.edu/ecology/
http://www.muohio.edu/botany/
"E Pluribus Unum"
More information about the R-help
mailing list