r - 以数据帧 + r + lme 的形式访问 Intervals.lme 的结果

Question

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

data(mtcars)

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

然后我使用 - 访问置信区间

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

另一种可能性是使用broom.mixed::tidy()。使用@TarJae 的示例:

broom.mixed::tidy(mlme, conf.int = TRUE, effects = "fixed") |> 
     dplyr::select(term, estimate, conf.low, conf.high)

# A tibble: 2 × 4
  term        estimate conf.low conf.high
  <chr>          <dbl>    <dbl>     <dbl>
1 (Intercept)    31.5     25.7      37.3 
2 wt             -3.52    -4.99     -2.04