如何从lmer / lme结果计算可靠性估计值？

Question

免责声明：很抱歉包含图片。我尝试以降价格式编写公式但未通过

我来自商业MLM（HLM7）软件，想在R中复制一些数字。

具体来说，我正在寻找一个函数或公式来计算 J 2级单位 的每个1级系数的最小二乘估计的可靠性p>

以下是基于简单sleepstudy数据的示例。我正在寻找的是一种计算可靠性值的方法，不仅在这个例子中，而且在有更多level1变量的情况下。

从HLM7手册（Raudenbush，Bryk（2002），第11页）给出了可靠性的定义：

分层线性模型中的公式3.58（第2版）

之后是以下概念：

我使用sleepstudy包中的lme4数据来计算lme4::lmer的随机拦截和斜率模型：

library(lme4)
m <- lmer(Reaction ~ Days + (Days|Subject), data = sleepstudy)
summary(m)

使用HLM7软件

固定和随机效应估计非常相似（舍入发生差异），但HLM7也将提供其可靠性估计值：

 ----------------------------------------------------
  Random level-1 coefficient   Reliability estimate
 ----------------------------------------------------
  INTRCPT1, G0                        0.730
      DAYS, G1                        0.815
 ----------------------------------------------------

这是我希望能够从lmer()结果获得的。

有什么想法吗？

非常感谢