pdf("F:\\sta4211\\KNNL_TA23_01.pdf") hormone <- read.table("http://www.stat.ufl.edu/~winner/sta4211/datasets/CH23TA01.txt",header=F, col.names=c("growrate","gender","bonedev","obsnum")) attach(hormone) gender <- factor(gender) levels(gender) <- c("Male", "Female") bonedev <- factor(bonedev) levels(bonedev) <- c("Severely Depressed", "Moderately Depressed", "Mildly Depressed") interaction.plot(bonedev,gender,growrate,xlab="Bone Development",ylab="Change in Growth Rate") tapply(growrate,list(gender,bonedev),mean) tapply(growrate,list(gender,bonedev),length) dev.off()