比较来自世界各地的卖家的域名和 IT 服务价格

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

我处理不平衡的设计 / 样本和最初研究过
aov//

. 现在我知道这是我的测试 ANOVA 我需要使用平方和 III, 其中包括使用配件使用
lm//

, 而不是使用
aov//

.

问题是接收特殊测试。 /特别是, tukey's HSD/ 使用
lm//

. 所有研究都花在了,谈到了使用
simint

在包装中
multcomp

它将起作用,但现在它是更新的,这个团队似乎无法使用。 他显然也依赖于段落
aov//

计算。

基本上所有测试都是 Tukey HSD, 我找到了哪个 R, 建议你正在使用
aov//

相比,不是
lm//

. 获得像这样的平方和 III, 我需要使用不平衡的设计,我必须使用:


mod<-lm/Snavg~StudentEthnicity*StudentGender/

Anova/mod, type="III"/


我如何使用测试 Tukey HSD 用我的mod使用
lm//

? 或者,相反,计算我的 ANOVA, 使用类型 III, 然而,可以开始测试 Tukey 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!

要回复问题请先登录注册