从 R 中的摘要对象中提取 df

Question

我在下面的 R 代码中想知道如何从 ems$emmeans 对象中提取 df 值？

library(lme4)
library(emmeans)
dat <- read.csv('https://raw.githubusercontent.com/hkil/m/master/z.csv')
dat$year <- as.factor(dat$year)

m1 <- lmer(y~year*group + (1|stid), data = dat)

ems <- emmeans(m1, pairwise ~ group*year, infer = c(T, T))

# How to extract 'df' values from ems$emmeans?

Answer 1

'emmGrid'

有as.data.frame方法

methods(class = 'emmGrid')
#[1] [             +             as.data.frame as.glht       as.list       CLD           coef          confint      
#[9] contrast      pairs         plot          predict       print         rbind         show          str          
#[17] summary       test          update        vcov

因此，我们可以将 'emmeans' 的输出转换为 data.frame 并使用 $ 或 [[

提取 'df'

as.data.frame(ems$emmeans)$df

有一个Note暗示`D.f。计算被禁用，因为观测值超过 3000，要启用，请添加 'pbkrtest.limit = 3435'

ems <- emmeans(m1, pairwise ~ group*year, infer = c(T, T), pbkrtest.limit = 3435)
as.data.frame(ems$emmeans)$df
#[1] 3090.821 2994.815 3328.634 3164.587 3380.889 3371.980

从 R 中的摘要对象中提取 df

Extract df from a summary object in R

statistics

regression

r

dataframe

lm