如何制作Tyuki测试 HSD 使用分散分析团队 /与汽车的服务包/

尝试

HSD.test

在

agricolae

library/agricolae/

data/sweetpotato/

model<-lm/yield~virus, data=sweetpotato/

comparison <- HSD.test/model,"virus", group=TRUE,

main="Yield of sweetpotato\nDealt with different virus"/

出口

Study: Yield of sweetpotato

Dealt with different virus



HSD Test for yield 



Mean Square Error: 22.48917 



virus, means



 yield std.err replication

cc 24.40000 2.084067 3

fc 12.86667 1.246774 3

ff 36.33333 4.233727 3

oo 36.90000 2.482606 3



alpha: 0.05 ; Df Error: 8 

Critical Value of Studentized Range: 4.52881 



Honestly Significant Difference: 12.39967 



Means with the same letter are not significantly different.



Groups, Treatments and means

a oo 36.9 

ab ff 36.33333 

 bc cc 24.4 

 c fc 12.86667

我发现

HSD.test//

也非常恪守地指的是你如何建造模型

lm//

或者

aov//

, 你为她使用了。

没有结论来自

HSD.test//

使用我的数据时，我使用以下编码想法

lm//

:

model<-lm/sweetpotato$yield ~ sweetpotato$virus/ 

 out <- HSD.test/model,"virus", group=TRUE, console=TRUE/

输出只有一个：

Name: virus 

 sweetpotato$virus

使用相同的逻辑时，结果同样差

aov//

model<-aov/sweetpotato$yield ~ sweetpotato$virus/

得出结论

HSD.test//

lm//

/或者也在使用时

aov//

对于模型 /
它应该严格使用响应中提出的逻辑构建 MYaseen208:

model <- lm/yield~virus, data=sweetpotato/

我希望这将有助于某人没有收到的结果

HSD.test//

.

作为初始通知，如果尚未更改以获得类型的平方和的正确结果 iii, 您需要建立因子变量的对比度编码。 它可以在呼叫内完成

lm

或者在帮助

options

. 使用以下示例

options

.

我会小心使用

HSD.test

除非文档涉及在这些情况下，否则具有不平衡结构的类似功能。 在文件中

TukeyHSD

提到它被调整了设计。 "mildly unbalanced". 我不知道如果

HSD.test

不同地联系事情。 您必须检查此功能给出的附加包文档或源链接。

作为一个符号，通过结束整个功能

HSD.test

在括号中，您将强制打印结果。 请参见下面的示例。

一般来说，我建议放弃

Tukey.HSD

和类似的功能和使用灵活包

emmeans

/nee.

lsmeans

/ 或者

multcomp

满足HOC比较的所有需求。

emmeans

特别有用
https://rcompanion.org/rcompanion/d_08.html
或者
https://rcompanion.org/rcompanion/h_01.html
.

使用不平衡设计，您可以报告ANEM E.M. /或者 L.S./ 而不是平均算术。 厘米
http://rcompanion.org/handbook/G_05.html
. 在下面的示例中注意，该报告的限制平均值

emmeans

, 与那些报告不同

HSD.test

.

还要注意 "Tukey" 在

glht

与比较无关 Tukey HSD 或者 tukey-adjusted; 它只是设置对比所有对测试的对比，正如他们在输出中所说的那样。

但是，功能

adjust="tukey"

在

emmeans

真的意味着使用比较与Tyuki的修正，因为他们在产出中说。

以下示例是部分借来的
https://rcompanion.org/rcompanion/d_05.html
分析 .

if/!require/car//{install.packages/"car"/}

library/car/

data/mtcars/

mtcars$cyl.f = factor/mtcars$cyl/

mtcars$carb.f = factor/mtcars$carb/



options/contrasts = c/"contr.sum", "contr.poly"//



model = lm/mpg ~ cyl.f + carb.f, data=mtcars/



library/car/

Anova/model, type="III"/



if/!require/agricolae//{install.packages/"agricolae"/}

library/agricolae/

/HSD.test/model, "cyl"/$groups/



if/!require/emmeans//{install.packages/"emmeans"/}

if/!require/multcompView//{install.packages/"multcompView"/}

library/emmeans/

library/multcompView/



marginal = emmeans/model,

 ~ cyl.f/



pairs/marginal, adjust="tukey"/



cld/marginal, adjust="tukey", Letters=letters/





if/!require/multcomp//{install.packages/"multcomp"/}

library/multcomp/



mc = glht/model,

 mcp/cyl.f = "Tukey"//



summary/mc, test=adjusted/"single-step"//



multcomp::cld/mc/

我遇到了同样的问题 HSD.test, 没什么打字。 你需要放置

console=TRUE

在功能内，以便自动打印它。

例如：

HSD.test/alturacrit.anova, "fator", console=TRUE/.

Hope it helps!