pdf("kuehl_06_05.pdf") temp <- rep(seq(25,35,10),each=18) density <- rep(seq(80,160,80),each=9,times=2) saline <- rep(seq(10,40,15),each=3,times=4) wtgain <- c(86,52,73,544,371,482,390,290,397, 53,73,86,393,398,208,249,265,243, 439,436,349,249,245,330,247,277,205, 324,305,364,352,267,316,188,223,281) shrimp <- data.frame(temp,density,saline,wtgain) attach(shrimp) temp <- factor(temp) density <- factor(density) saline <- factor(saline) tapply(wtgain,list(temp,density,saline),mean) tapply(wtgain,list(temp,density),mean) tapply(wtgain,list(temp,saline),mean) tapply(wtgain,list(density,saline),mean) shrimp.mod1 <- aov(wtgain ~ temp*saline*density) summary(shrimp.mod1) dev.off()