Displaying 1 result from an estimated 1 matches for "massout".
Did you mean:
massot
2010 Oct 18
1
Crossed random effects in lme
...hannel (ch)
as random factors in a model trying to explain variation in oxygen
consumption. From Pinheiro & Bates (2000), I've gathered that this model
could be fit making use of pdBlocked and pdIdent, so I've tried fitting the
below model:
m1.bmr<-with(bmred.df,lme(mlmin~treat*year+massout,random=pdBlocked(list(pdIdent(~hack-1),pdIdent(~ch-1)))
))
However, my model fails with the following error message:
Error in getGroups.data.frame(dataMix, groups) :
Invalid formula for groups
I would much appreciate any input on this!
Kind regards,
Andreas Nord
Sweden
--
View this me...