Have you tried GLMM in lme4? Doug Bates is the primary architect
of both nlme and lme4. Therefore, I would think that a spatial
dependence model that works in nlme might also work in GLMM.
hope this helps. spencer graves
Martin Henry H. Stevens wrote:
Running Mac OS 10.3.5 and R 2.0
Does glmmPQL use the same spatial dependence models as gls in nmle? It
does not seem to - I get the following, for example:
> m2 <- glmmPQL(S.Early ~ fertilized*watered, data=geodat,
family="poisson", random=~1|col)
iteration 1
iteration 2
> plot(Variogram(m2, form=~col+row))
> m3 <- update(m2,correlation=corSpher(c(5,.5), form=~col+row, nugget=T))
iteration 1
Error in model.frame(formula, rownames, variables, varnames, extras,
extranames, :
invalid variable type
>
Any leads on how to specify exponential or spherical models would be
appreciated (e.g., pageg or chapter in Venables and Ripley 2002).
Cheers,
Hank
Dr. Martin Henry 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/botany/bot/henry.html
http://www.muohio.edu/ecology/
http://www.muohio.edu/botany/
"E Pluribus Unum"
______________________________________________
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html
--
Spencer Graves, PhD, Senior Development Engineer
O: (408)938-4420; mobile: (408)655-4567
______________________________________________
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html