R-lmer以随机效应降低未使用的等级

Question

我在lFormula（lmer 1.1-9）中调用lme4函数来获得随机效果（Zt）的稀疏模型矩阵的转置。问题是输出矩阵中未使用的级别被丢弃。有没有办法解决这个问题？

我意识到这个问题可能看起来很尴尬，所以here就是GitHub Gist的一个例子，试图让多成员模型适应一个比观察更多级别的因子。