> anovakun(dat, "sA", 3, auto=T) [ sA-Type Design ] This output was generated via anovakun 4.3.3 at R version 2.15.2. It was executed on Thu Jun 06 15:24:17 2013. << DESCRIPTIVE STATISTICS >> -------------------------------------------- A N Mean S.D. -------------------------------------------- a1 291 8.9794 3.0114 a2 291 23.0997 4.0938 a3 291 15.0344 5.4822 -------------------------------------------- << SPHERICITY INDICES >> == Mendoza's Multisample Sphericity Test and Epsilons == ---------------------------------------------------------------------------------- Effect Lambda approx.Chi df p LB GG HF ---------------------------------------------------------------------------------- A NaN NaN 2 0.5000 0.8504 0.8549 ---------------------------------------------------------------------------------- LB = lower.bound, GG = Greenhouse-Geisser, HF = Huynh-Feldt << ANOVA TABLE >> == Adjusted by Greenhouse-Geisser's Epsilon for Suggested Violation == ----------------------------------------------------------------------------------------------------- Source SS df MS F-ratio p-value ----------------------------------------------------------------------------------------------------- s 8480.4192 290 29.2428 ----------------------------------------------------------------------------------------------------- A 29206.1100 2 14603.0550 1096.3789 0.0000 *** sxA 7725.2234 580 13.3194 ----------------------------------------------------------------------------------------------------- Total 45411.7526 872 +p < .10, *p < .05, **p < .01, ***p < .001 << POST ANALYSES >> < MULTIPLE COMPARISON for FACTOR A > == Shaffer's Modified Sequentially Rejective Bonferroni Procedure == == The factor < A > is analysed as dependent means. == == Alpha level is 0.05. == -------------------------------------------- A N Mean S.D. -------------------------------------------- 1 291 8.9794 3.0114 2 291 23.0997 4.0938 3 291 15.0344 5.4822 -------------------------------------------- ------------------------------------------------------------------------------------------------- Pair Interval t-value df p adj.p ------------------------------------------------------------------------------------------------- a1-a2 -14.1203 61.2380 290 0.0000 0.0000 a1 < a2 * a2-a3 8.0653 24.3587 290 0.0000 0.0000 a2 > a3 * a1-a3 -6.0550 18.1067 290 0.0000 0.0000 a1 < a3 * ------------------------------------------------------------------------------------------------- output is over --------------------///