R Levene和#x27的多重比较事后测试；s检验

我想在R中对Levene的测试做一个成对比较事后测试。我知道如何使用PROC GLM在SAS中进行测试，但我似乎不知道如何在R中进行测试。有人知道吗？在下面的例子中，我希望能够测试所有级别的“cat”之间方差的同质性，即A-B、A-C、A-d、B-C、B-d、C-d。我发现的最好的方法是将我的数据子集到这些对中的每一对，然后对每个子集运行Levene测试，最后进行Bonferroni校正。然而，当我的因素数量变大时，这不是一个实际的解决方案

library(car)
dat <- rnorm(100, mean=50, sd=10)
cat <- rep(c("A", "B", "C","D"), each=25)
df <- data.frame(cat,dat)
df$cat <- as.factor(df$cat)

LT <- leveneTest(dat ~ cat, data = df)

库（车）
dat由于Levene检验只是对样本方差（残差）而非样本均值进行的方差分析，因此您可以手动计算残差，然后使用TukeyHSD检验作为事后检验运行方差分析

首先，正如标题所示，多重比较：
以您的示例为例，使用额外的因子（cat2），我们也可以进行交互：

df <- df %>% group_by(cat, cat2) %>% 
  mutate(dat.med = ifelse(dat,median(Ctmax, na.rm=TRUE), ifelse(dat==NA, NA)))

对于两级因子的成对比较（使用包PairedData）：

您能否在
PROC GLM
/中向我们展示您将如何在SAS在线文档中描述此选项？我实际上没有SAS可以运行它，因此需要R。但是，我相信编辑我的标题时，您已经忽略了我的问题与R相关的事实。。。不是任何其他语言。我之所以提到SAS，是因为在一些旧的注释中，我有如下内容：PROC GLM data=dataset；表示classvar/HOVTEST=LEVENE/tukey；我实际上没有SAS来运行它，因此需要R。然而，我相信编辑我的标题，你已经从我的问题与R有关的事实中拿走了。。。不是任何其他语言。我之所以提到SAS，是因为在一些旧的注释中，我有如下内容：PROC GLM data=dataset；类变量；模型Y=classvar；表示classvar/HOVTEST=LEVENE/tukey；另外，在编辑我的代码时，我觉得你让它变得不那么可读。。。将来，如果你不介意的话，我不会认为对你来说最简单的方式对每个人来说都是最简单的。（1）这是标准的，因此标签中的信息（例如“在R中”）不应该也出现在标题中：（2）缩进：抱歉，欢迎你回滚我的更改。。。
df$dat.med.res<-abs(df$dat-df$dat.med)

# Then we run an ANOVA, and post-hoc if necessary:
levene.dat.aov<-aov(dat.med.res~cat*cat2,df)
summary(levene.dat.aov)
TukeyHSD(levene.dat.aov)

 aov(dat.med.res~cat+Error(Subject/(cat)),df)

levene.var.test(df$dat[df$cat=="A"], df$dat[df$cat=="B"],location=c("median"))