#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

sse.mat <- matrix(rep(0,82),ncol=2)
row <- 0
for (i in seq(-2,2,0.1)) {
row <- row+1
if (i == 0) ty <- log(tcount)
else ty <- (tcount^i-1)/i

e <- residuals(lm(ty~site))

sse.mat[row,1] <- i
sse.mat[row,2] <- sum(e^2)

}

sse.names <- c("power","sse")
colnames(sse.mat) <- sse.names

sse.mat

dev.off()