I have an ANOVA with a numerical score as the DP and a category as an IV:
model = aov(score~group, data = dat) And then I run a Tukey:
TukeyHSD(model). I get the output:
diff lwr upr p adj no-maybe 6.4962406 -1.667561 14.660042 0.1473910 yes-maybe 5.6575851 1.191148 10.124022 0.0087644 yes-no -0.8386555 -7.959717 6.282406 0.9582525
But i need to report mean and SD, how can I get this information?
Edit: I already tried that agricolae package but it won't load for me