Criando uma Função

No R(), pode-se criar uma função que realiza este cálculo para um experimento em bolcos completos casualizados:

function(n2,s,t,d,ri,alpha=0.95)
{
## Cálculo do número de repetições para um experimento em BCC

## n2 = graus de liberdade do resíduo de um experimento anterior
## s  = desvio padrão da variável resposta
## t  = número de tratamentos
## d  = diferença significativa
## ri = número inicial de repetições 

df.res<-(ri*t-1)-(ri-1)-(t-1)
q<-qtukey(alpha,t,df.res)
f<-qf(alpha,df.res,n2)
r<-(q^2*s^2*f)/(d^2)
return(r,df.res,q,f)
}

Para um experimento conduzido em outro delineamento, baste redefinir df.res.