pdf("C:\\Rmisc\\graphs\\kuehl_04_03.pdf") site <- rep(1:6,each=25) count <- c(0,0,22,3, 17,0,0,7,11,11,73,33,0,65,13,44,20,27,48,104,233,81,22,9,2, 415,466,6,14,12,0,3,1,16,55,142,10,2,145,6,4,5,124,24,204,0,0,56,0,8, 0,0,4,13,5,1,1,4,4,36,407,0,0,18,4,14,0,24,52,314,245,107,5,6,2, 0,0,0,4,2,2,5,4,2,1,0,12,1,30,0,3,28,2,21,8,82,12,10,2,0, 0,1,1,2,2,1,2,29,2,2,0,13,0,19,1,3,26,30,5,4,94,1,9,3,0, 0,0,0,2,3,0,0,4,0,5,4,22,0,64,4,4,43,3,16,19,95,6,22,0,0) crab <- data.frame(site,count) attach(crab) site <- factor(site) tcount <- count + 1/6 ty <- ((tcount^-0.5)-1)/(-0.5) tcrab.mod <- aov(ty~site) anova(tcrab.mod) e <- residuals(tcrab.mod) ybar.sort <- sort(predict(tcrab.mod)) ybar.sort i <- 1:150 ff_i <- (i-0.5)/150 nqf_i <- qnorm(ff_i,0,1) par(mfrow=c(1,2)) plot(nqf_i,ybar.sort-mean(ty),ylim=c(-4,4)) qqnorm(e,ylim=c(-4,4)) qqline(e) dev.off()