R语言实现常用多重比较方法

R语⾔实现常⽤多重⽐较⽅法

⽂章⽬录

在单因素⽅差分析 ANOVA中，如果该因素影响⽐较显著，那么需要进⼀步利⽤多重⽐较⽅法⽐较该因素不同⽔平的影响，确定不同⽔平下该因素的影响是否显著。常见的多重⽐较⽅法主要有两种， LSD法和 Tukey HSD法。下⾯对R语⾔中，这两种多重⽐较⽅法的实现进⾏举例。

前期数据如下，影响因素为group，指标为value：

> head(tarD)

value group sample time

A0522W11NC1 0.0002053745 normal A0522W11NC1 11week

A0522W11NC2 0.0031773712 normal A0522W11NC2 11week

我爱土豆

A0522W11NC3 0.0060378288 normal A0522W11NC3 11week

极限空间下载

A0522W11NC4 0.0017626931 normal A0522W11NC4 11week

A0522W11NC5 0.0018035261 normal A0522W11NC5 11week

A0522W11NC6 0.0036690067 normal A0522W11NC6 11week

> tmp <- aov(value ~ group, tarD)

注意，⽆论是LSD还是Tukey进⾏多重⽐较，都需要先进⾏ANOVA分析，也就是上例中的tmp变量

最⼩显著差数检验法(LSD法)

> library(agricolae)荷花淀教学反思

> res <- st(tmp, 'group', p.adj = 'bonferroni')

> print(res$groups)

trt means M

1 normal 2.576910e-03 a

2 drug

3 7.552555e-0

4 bwikipedia

3 drug2 7.269247e-05 b

掌握学习理论4 high_fat 6.220610e-0

5 b

5 drug1 2.954733e-05 b

Tukey⽒固定差距检验法(Tukey HSD)

> TukeyHSD(tmp)

Tukey multiple comparisons of means

95% family-wise confidence level

Fit: aov(formula = value ~ group, data = tarD)

$group

diff lwr upr p adj

drug2-drug1 4.314514e-05 -0.0015468705 0.0016331608 0.9999916 drug3-drug1 7.257082e-04 -0.0008643074 0.0023157239 0.6929965 high_fat-drug1 3.265877e-05 -0.0015149488 0.0015802664 0.9999969 normal-drug1 2.547362e-03 0.0009997549 0.0040949700 0.0002613 drug3-drug2 6.825631e-04 -0.0009487586 0.0023138847 0.7563196 high_fat-drug2 -1.048637e-05 -0.0016005020 0.0015795293 1.0000000 normal-drug2 2.504217e-03 0.0009142017 0.0040942330 0.0004945 high_fat-drug3 -6.930494e-04 -0.0022830651 0.0008969662 0.7277757 normal-drug3 1.821654e-03 0.0002316386 0.0034116699 0.0175538 normal-high_fat 2.514704e-03 0.0009670961 0.0040623113 0.0003161

> TukeyHSD(tmp)$group

diff lwr upr p adj

drug2-drug1 4.314514e-05 -0.0015468705 0.0016331608 0.9999915820 drug3-drug1 7.257082e-04 -0.0008643074 0.0023157239 0.6929965170 high_fat-drug1 3.265877e-05 -0.00151

Obama speech in Shanghai49488 0.0015802664 0.9999969171 normal-drug1 2.547362e-03 0.0009997549 0.0040949700 0.0002612744 drug3-drug2 6.825631e-04 -0.0009487586 0.0023138847 0.7563195891 high_fat-drug2 -1.048637e-05 -0.0016005020 0.0015795293 0.9999999705 normal-drug2 2.504217e-03 0.0009142017 0.0040942330 0.0004944674 high_fat-drug3 -6.930494e-04 -0.0022830651 0.0008969662 0.7277757202 normal-drug3 1.821654e-03 0.0002316386 0.0034116699 0.0175537862 normal-high_fat 2.514704e-03 0.0009670961 0.0040623113 0.0003161003

参考

本文发布于:2024-09-21 04:22:57，感谢您对本站的认可！

本文链接：https://www.17tex.com/xueshu/40863.html

上一篇：肾上腺髓样脂肪瘤与肾上腺皮质腺瘤的影像学诊断比较

下一篇：针灸联合神经肌肉电刺激、导管球囊

标签：因素影响检验法分析

留言与评论（共有 0 条评论）