s2a <- 6.81
s2e <- 5.82
t <- 5
alpha <- 0.05
r_max <- 30

r_power <- matrix(rep(0,5*(r_max-1)),ncol=5)



trtdf <- t-1

for (r in 2:r_max) {
errdf <- t*(r-1)
f_alpha <- qf(1-alpha,trtdf,errdf)
power <- 1-pf((s2e/(s2e+r*s2a))*f_alpha,trtdf,errdf)
r_power[r-1,] <- cbind(r,trtdf,errdf,f_alpha,power)
}

power.names <- c("r","trtdf","errdf","f_0.05","power")
colnames(r_power) <- power.names

r_power