Tengo el siguiente problema
Sea X1 , . . . , Xn una muestra aleatoria de una población con distribución N(μ,σ^2).
Si σ^2 es desconocida. Encuentre el valor mínimo de n que garantice que, con probabilidad 0.90, el intervalo del 95 % de confianza para μ tendrá longitud no mayor que σ/4.
Para lo cual escribí el siguiente código:
for(n in 1:1000){
z <- 64*qchisq(0.9, n-1)*qt(0.975, n-1)^2
if((n*(n-1))>=z)
print(n)
break
}
O sea con un bucle for recorror las n y con un if pongo la condición para encontrar a la indicada, pero al momento de correrlo me lanza el error
Warning message: NaNs produced
Error: missing value where TRUE/FALSE needed
>
¿me pide un booleano?, ¿en donde debería ir?