对于MANOVA具有变量的 a,我想在个级别之间进行成对比较。
多重比较时,采用什么合适的方法?
多个测试以及校正或 FDR/pFDR 是否合适?FDR/pFDR q 值更可取,因为值在这里很重要。
HotellingBonferroniR对包做同样的任何建议?(特别是对于 MANOVA 事后多重比较)如何在多重比较的等价检验中而不是
编辑
根据rvl的回答和进一步评论,我能够探索并提出以下建议。
library(lsmeans)
# Use the `oranges` dataset in `lsmeans` package.
# multivariate linear model
oranges.mlm <- lm(cbind(sales1,sales2) ~ price1 + price2 + day + store,
data = oranges)
# Get the least square means
oranges.Vlsm <- lsmeans(oranges.mlm, "store")
# Multiple comparisons with fdr p value adjustment
test(contrast(oranges.Vlsm, "pairwise"), side = "=", adjust = "fdr")
# With threshold spcified
test(contrast(oranges.Vlsm, "pairwise"), side = "=", adjust = "fdr", delta = 0.25)