データセット (サンプルサイズが等しくない、分散が等しくない) を使用して、r で双方向分散分析を実行しました: 3 つの種で測定された 1 つの変数 (各種のオスとメス)。これは種間で有意な結果をもたらすため、どのペアワイズ比較が有意性を生み出すかを知りたい. Rで事後テストを実行するためのパッケージに関数があることを私は知っています：例えば
http://www.uwlax.edu/faculty/toribio/math305_fall09/multiple.txtからの Dunnett の事後検定。必要なパッケージ: 「multcomp」、「mvtnorm」、「survival」、「splines」
library(multcomp)
test.dunnett=glht(anova_results,linfct=mcp(method="Dunnett"))
confint(test.dunnett)
plot(test.dunnett)
*注: glht は「multcomp」に記述されています
ただし、Dunett の検定は、すべてのグループを対照と比較するように設計されています。代わりに、すべてのグループを相互に比較したいと思います。Dunnett C です。Dunnett C を実行するパッケージを知っているか、それをコーディングする方法を知っている人はいますか? (式: http://pic.dhe.ibm.com/infocenter/spssstat/v21r0m0/index.jsp?topic=%2Fcom.ibm.spss.statistics.help%2Falg_posthoc_unequalvar.htm )