以下に部分的に示すデータセットがあります。私は anova を使用して、共変量として地域を持つ値の性別間の違いを調べました。ただし、複数の比較の修正に TukeyHSD を使用すると、必要以上の比較が表示されることがわかりました (以下を参照)。私が興味を持っている唯一の部分は、リージョン間ではなく、最初の行 (つまり、各リージョン内の M 対 F) です。anovaモデルまたはTukeyでこれを指定する方法はありますか? ありがとう！
ID Age Sex Hemi Region value
1 62 M R a 1.81
2 62 M R a 1.90
3 72 M R a 2.25
1 61 M L a 1.58
2 57 F L a 2.66
3 62 M L a 2.19
1 72 M R b 1.93
2 64 F R b 1.07
3 65 F R b 1.37
1 64 M L b 0.97
2 77 F L b 1.59
3 27 M L b 1.84
model=aov(value~Sex*Region,data=data
TukeyHSD(model)
$Sex:as.factor(Region)
diff lwr upr p adj
M:a-F:a 0.1125161943 -0.2989139 0.5239462408 1.0000000
F:b-F:a -0.0383423077 -0.4866878 0.4100031910 1.0000000
M:b-F:a 0.1273767206 -0.2840533 0.5388067671 1.0000000