我在 lmerlme4 1.1-9)中调用 lFormula 函数来获得随机效应稀疏模型矩阵的转置( Zt ) . 问题是输出矩阵中未使用的级别被丢弃 . 有没有办法解决这个问题?

我意识到这个问题可能看起来很尴尬,所以here是一个关于GitHub Gist的例子,试图使用一个比观察更多级别的因子来拟合多成员模型 .