R lme4故障中的提取器功能？

我试图理解

R

中

lmer

函数背后的线性代数运算，我在

lme4

软件包的创建者道格拉斯·贝茨（Douglas Bates）的一篇文章中在线找到了一个很好的资源

该示例处理数据集染料，并调用混合效果模型，如下所示：

fm1 <- lmer(Yield ~ 1 + (1 | Batch), Dyestuff)

或者

执行

？expand

会将此函数标识为仍然存在，似乎是为了生成矩阵分解，如LU或RQ

在网上搜索时，我发现道格拉斯现在使用的是Julia（下一种统计语言能有一个不那么令人难以置信的名字吗？不，不是“Pied Pier”！对不起，我离题了…）

---

我做错了什么？

lme4

现在是孤立的并且正在衰退吗？幻灯片中有打字错误吗？

您看到的是非常旧的幻灯片

• 包装的线性混合模型功能的完整描述，包括内部表示的详细信息，见J.Stat.软件中的press（参考此信息也在中）。（我们仍在撰写描述GLMM功能的论文。）

• getME（）

  函数是当前推荐的访问模型信息的方法

• lme4

  仍处于相当活跃的开发阶段
• 道格·贝茨确实更感兴趣，但在某种程度上，他仍然参与了

  lme4

  维护

表示法/内部表示法有所改变，但从内部信息重构方差协方差矩阵的方法如下（内部

Lambdat

相当于旧表示法中的

t（t%*%S）

）

库（“lme4”）
fm1
efm1 <- expand(fm1)
efm1$S
6 x 6 diagonal matrix of class "ddiMatrix"

#       [,1]      [,2]    [,3]    [,4]    [,5]    [,6]
[1,]  0.84823       .       .       .       .       .
[2,]       .  0.84823       .       .       .       .
[3,]       .       .  0.84823       .       .       .
[4,]       .       .       .  0.84823       .       .
[5,]       .       .       .       .  0.84823       .
[6,]       .       .       .       .       .  0.84823

efm1$T6 x 6 sparse Matrix of class "dtCMatrix"
[1,] 1 . . . . .
[2,] . 1 . . . .
[3,] . . 1 . . .
[4,] . . . 1 . .
[5,] . . . . 1 .
[6,] . . . . . 1

(fm1S <- tcrossprod(efm1$T %*% efm1$S))

6 x 6 sparse Matrix of class "dsCMatrix"

[1,] 0.71949 .       .       .       .       .
[2,] .       0.71949 .       .       .       .
[3,] .       .       0.71949 .       .       .
[4,] .       .       .       0.71949 .       .
[5,] .       .       .       .       0.71949 .
[6,] .       .       .       .       .       0.71949

efm1 <- expand(fm1)

Error in (function (classes, fdef, mtable)  : 
  unable to find an inherited method for function ‘expand’ for signature ‘"lmerMod"’

efm1$S
Error: object 'efm1' not found

library("lme4")
fm1 <- lmer(Yield ~ 1 + (1 | Batch), Dyestuff)
crossprod(getME(fm1,"Lambdat"))