在进行大的ANOVA时，有时我们需要使用Tukey HSD测试来验证配对显着性。由于这是一项耗时的任务，想请教一下如何进行部分测试。我的意思是，我知道摘要后的重要项目（aov），我只想对主要项目进行Tukey HSD，并摆脱其余部分。以下是参考代码：
seed(100)
# create the DF for aov
X <- rep(letters[1:3],each=6)
Y <- rep(sample(100,3),each=6)
Z <- rep(rep(LETTERS[3:4],9))
m <- data.frame(X,Y,Z)
# do aov and summary
fit <- aov(Y~X+Z, data=m)
summary(fit)
# Df Sum Sq Mean Sq F value Pr(>F)
# X 2 7756 3878 1.329e+30 <2e-16 ***
# Z 1 0 0 7.260e-01 0.409
# Residuals 14 0 0
# ---
# Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
TukeyHSD(fit)
# Tukey multiple comparisons of means
# 95% family-wise confidence level
# Fit: aov(formula = Y ~ X + Z, data = m)
# $`X`
# diff lwr upr p adj
# b-a 4 4 4 0
# c-a -1 -1 -1 0
# c-b -5 -5 -5 0
# $Z
# diff lwr upr p adj
# D-C 0 -2.363621e-15 2.363621e-15 1
我可以做X而不是Z吗？类似于TukeyHSD(fit[["X"]])？
$'X'，但不是$ Z？非常感谢。