我使用包中的lmer()函数构建了一个混合效果模型lme4。由于lme4某些良好的哲学原因，该程序包不输出系数的p值。但是，我仍然需要p值才能在出版物中进行报告。我知道有多种方法可以使用lmer()，例如here创建的模型来计算p值。
我的问题是：我想使用包中的tidy()函数提取p值broom。在这里，我确实要坚持下去，tidy()因为我想维护以下管道：
data_frame %>% group_by(grouping variables) %>% do(tidy(fitted_model))
一种选择是创建一个自定义函数，并将其附加到管道中。但是，该broom软件包的手册页（http://rpackages.ianhowson.com/cran/broom/man/lme4_tidiers.html）表示：
"p.value P-value computed from t-statistic (may be missing/NA)".
这样，我假设已经在扫帚中实现了根据lmer模型给出的t值计算p值的函数。因此，我不愿意重新发明轮子。
问题是我根本没有得到名称为p.value的列。我原本以为最坏的情况是NA称为p.value的列。
代码：
library(lme4)
library(broom)
lme <- lmer(Reaction ~ Days + (1 + Days | Subject), sleepstudy)
tidy(lme)
tidy(lme, effects = "fixed")
输出：
> tidy(lme)
term estimate std.error statistic group
1 (Intercept) 251.40510485 6.824557 36.838306 fixed
2 Days 10.46728596 1.545789 6.771485 fixed
3 sd_(Intercept).Subject 24.74045195 NA NA Subject
4 sd_Days.Subject 5.92213312 NA NA Subject
5 cor_(Intercept).Days.Subject 0.06555113 NA NA Subject
6 sd_Observation.Residual 25.59181564 NA NA Residual
> tidy(lme, effects = "fixed")
term estimate std.error statistic
1 (Intercept) 251.40510 6.824557 36.838306
2 Days 10.46729 1.545789 6.771485
您将需要包装lmerTest以获得p值。tidy将不适用于该lme对象，因此您需要将其附加到您的格式中。
attach(mtcars)
lme <- lmer(mpg ~ cyl + (1 + cyl | carb), mtcars)
summary(lme)
本文收集自互联网，转载请注明来源。
如有侵权，请联系[email protected] 删除。
我来说两句