作为数据框访问 intervals.lme 的结果 + r + lme

Question

我正在使用线性混合效应模型来确定 2 个变量随时间变化的关系，模型 -

mlme <- lme(RBS.Score.Endorsed ~ DairyTotal, random = ~ Period|ID, data = df)

然后我使用 -

访问置信区间

conf <- intervals(mlme, level = 0.95, which = 'fixed')
conf

返回结果-

Approximate 95% confidence intervals

 Fixed effects:
               lower      est.    upper
(Intercept) 7.431921 18.416639 29.40136
DairyTotal  2.001397  6.716849 11.43230

我想访问下限和上限的值，以便将它们用作变量，但找不到这样做的方法。例如，如果您想访问系数，您可以使用

coef(summary(mlme))

并将其转换为数据帧。

我尝试使用

data.frame(coef)

将结果转换为数据框，但出现错误：

> data.frame(conf)
Error in as.data.frame.default(x[[i]], optional = TRUE, stringsAsFactors = stringsAsFactors) : 
  cannot coerce class ‘"intervals.lme"’ to a data.frame

谁能推荐一种方法来访问以这种格式返回的变量？

Answer 1

这里是 mtcars 数据集的一个例子，我们如何访问较低和较高的 ci 作为变量：

library(nlme)

data(mtcars)

# linear mixed effects model
mlme <- lme(mpg ~ wt, random = ~ 1|cyl, data = mtcars)

# Confidence intervals
conf <- intervals(mlme, level = 0.95, which = 'fixed')

# Lower and upper confidence intervals as variables
lower <- conf$fixed["(Intercept)", "lower"]
upper <- conf$fixed["(Intercept)", "upper"]

[1] 25.7204
[1] 37.32919