pdf("E:\\sta6208\\poker_skill.pdf") poker <- read.table("http://www.stat.ufl.edu/~winner/data/poker_skill.dat", header=F,col.names=c("skill","hand","limit","cash")) attach(poker) skill <- factor(skill); levels(skill) <- c("Expert", "Novice") hand <- factor(hand); levels(hand) <- c("Bad","Neutral","Good") limit <- factor(limit); levels(limit) <- c("Fixed","None") (tapply(cash,list(skill,hand,limit),mean)); (tapply(cash,list(skill,hand),mean)) (tapply(cash,list(skill,limit),mean)); (tapply(cash,list(hand,limit),mean)) (tapply(cash,skill,mean)); (tapply(cash,hand,mean)); (tapply(cash,limit,mean)) par(mfrow=c(2,2)) interaction.plot(hand[skill=="Expert"],limit[skill=="Expert"],cash[skill=="Expert"], ylim=c(0,18),xlab="hand", ylab="Cash Balance",main="Experts") interaction.plot(hand[skill=="Novice"],limit[skill=="Novice"],cash[skill=="Novice"], ylim=c(0,18),xlab="hand", ylab="Cash Balance",main="Novices") interaction.plot(hand[limit=="Fixed"],skill[limit=="Fixed"],cash[limit=="Fixed"], ylim=c(0,18),xlab="hand", ylab="Cash Balance",main="Fixed Limit") interaction.plot(hand[limit=="None"],skill[limit=="None"],cash[limit=="None"], ylim=c(0,18),xlab="hand", ylab="Cash Balance",main="No Limit") options(contrasts=c("contr.sum","contr.poly")) poker.aov1 <- aov(cash ~ skill*hand*limit) anova(poker.aov1) summary.lm(poker.aov1) dev.off()