r - R - 最小二乘均值对比单向方差分析
问题描述
我正在审查一种方法方差分析并尝试整合最小二乘均值。这是 mtcars 的一个例子。
mtcars.mod <- mutate(mtcars, cyl.chr = case_when(
cyl == 4 ~ "A",
cyl == 6 ~ "B",
cyl == 8 ~ "C"
))
library(lsmeans)
model <- lm(mpg ~ cyl.chr, data = mtcars.mod)
lsmeans(model,
~ cyl.chr,
adjust = "sidak")
我的输出是这样的:
cyl.chr lsmean SE df lower.CL upper.CL
A 26.7 0.972 29 24.2 29.1
B 19.7 1.218 29 16.7 22.8
C 15.1 0.861 29 12.9 17.3
我试图得到看起来像这样的东西(值不能反映真实数据;它们是来自https://rcompanion.org/handbook/G_06.html的填充物/示例):
$contrasts
contrast estimate SE df z.ratio p.value
A - B 4.943822 1.3764706 NA 3.5916658 0.0010
A - C 0.633731 0.9055691 NA 0.6998152 0.7636
B - C -4.310091 1.3173294 NA -3.2718403 0.0031
P value adjustment: tukey method for comparing a family of 3 estimates
### Remember to ignore “estimate” and “SE” of differences with CLM,
### as well as “lsmeans” in the output not shown here
我错过了什么?
解决方案
这就是功能的工作tukeyHSD
:
TukeyHSD(aov(model))
Tukey multiple comparisons of means
95% family-wise confidence level
Fit: aov(formula = model)
$cyl.chr
diff lwr upr p adj
B-A -6.920779 -10.769350 -3.0722086 0.0003424
C-A -11.563636 -14.770779 -8.3564942 0.0000000
C-B -4.642857 -8.327583 -0.9581313 0.0112287
推荐阅读
- c++ - 基于 CLRS 的合并排序 C++ 上的算法介绍,带倒数计数
- math - 十六进制(数字)所需的数学指导
- python - 在熊猫数据框中为给定的列值创建顺序 ID
- reactjs - 在 React 中对 props 的属性进行 Props 验证
- django - 如何在外部 python 脚本 django 中导入模块
- craftcms - 如何添加资产字段
- php - PHP - 在 3 维数组中循环提供有线输出
- google-apps-script - 从一个电子表格更新到另一个电子表格时附加多行
- angular - 寻找使用 ngrx 8 沉浸式生产的示例
- javascript - 未填充 Google Analytics(分析)自定义维度