Displaying 1 result from an estimated 1 matches for "pdconstruct".
Did you mean:
construct
2004 Mar 01
0
question about mixed effects model
...x1 - fixed factor, x2 - random factor). Fixed
factor x1 is repeated twice for each x2 value.
I believe that the correct way to fit a mixed model would be
ddd <- groupedData(y ~ x1 | x2, dd)
l1 <- lme(y ~ x1, random = pdBlocked(list(pdIdent(~x2-1))), data = ddd)
But I get an error
Error in pdConstruct.pdBlocked(object, form = form, nam = nam, data =
data, : None of the arguments specify more than one block
I would also have to test for significance of each variable (x1 and x2)
and the way to do this seem to be to fit a model without x1 and compare
it to the original l1 model (and then the...