> (tapply(cash,list(skill,hand,limit),mean)); (tapply(cash,list(skill,hand),mean)) , , Fixed Bad Neutral Good Expert 7.3292 9.85 12.1192 Novice 5.0504 9.80 15.7996 , , None Bad Neutral Good Expert 6.8892 11.8488 12.3912 Novice 3.8396 9.4500 11.6300 Bad Neutral Good Expert 7.1092 10.8494 12.2552 Novice 4.4450 9.6250 13.7148 > (tapply(cash,list(skill,limit),mean)); (tapply(cash,list(hand,limit),mean)) Fixed None Expert 9.766133 10.376400 Novice 10.216667 8.306533 Fixed None Bad 6.1898 5.3644 Neutral 9.8250 10.6494 Good 13.9594 12.0106 > (tapply(cash,skill,mean)); (tapply(cash,hand,mean)); (tapply(cash,limit,mean)) Expert Novice 10.07127 9.26160 Bad Neutral Good 5.7771 10.2372 12.9850 Fixed None 9.991400 9.341467 > poker.aov1 <- aov(cash ~ skill*hand*limit) > anova(poker.aov1) Analysis of Variance Table Response: cash Df Sum Sq Mean Sq F value Pr(>F) skill 1 49.2 49.17 2.8391 0.093077 . hand 2 2646.6 1323.28 76.4125 < 2.2e-16 *** limit 1 31.7 31.68 1.8294 0.177258 skill:hand 2 219.0 109.51 6.3237 0.002052 ** skill:limit 1 119.1 119.11 6.8779 0.009191 ** hand:limit 2 97.3 48.64 2.8089 0.061922 . skill:hand:limit 2 42.4 21.19 1.2237 0.295650 Residuals 288 4987.5 17.32 ---