lme中用户指定的Z矩阵

Question

我一直在寻找如何在R中做到这一点并且找不到任何东西！基本上，我想要使用LMM缩小预测变量。所以我有一组固定的效果，X，我有一组预测变量Z，我想要随机效果，所以模型是

Y=X*beta+Z*u+e

其中u~N（0，sigma_u ^ 2 * I）和e~N（0，sigma_e ^ 2 * I）。我以为我可以用lme做到这一点

fit <- lme(Y~X,random=pdIdent(~-1+Z))

但我只收到错误：

Error in getGroups.data.frame(dataMix, groups) : 
invalid formula for groups

非常感谢有关此问题的任何帮助。

Answer 1

你试过了吗？

N =样本量

group <- rep(1, N)
fit <- lme(Y~X, random=list(group=pdIdent(~-1+Z)))