r - 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)))