Hola,
¿Puedes explicar un par de iteraciones?
Tengo estas dudas:
- En T y T2 ya utilizas "n", pero "n" lo declaras después
n=seq(10,900,70)... ??.
- Por otro lado en el bucle iteras desde 1:Nr (1 hasta 1000) y más
adelante cuando declaras "x" (x <- rnorm(n[i], Mx, sigma)
escoges el valor
de "n[i]" pero "n" tiene una longitud de 13 elementos no
de 1000
elementos... ??.
Saludos,
Carlos Ortega
www.qualityexcellence.es
El 5 de noviembre de 2012 20:32, Lex B
<bernalm.alexander@gmail.com>escribió:
> Buen día para todos los que ven este correo.
>
> En estos días estoy estudiando la potencia de una prueba en inferencia,
> pero no entiendo bien como
> correr un codigo que lo haga pertinentemente, el codigo que corro es el
> siguiente, la idea es que a medida que crece
> el tamaña muestral pues la potencia de la prueba aumente, pero no se como
> hacer para que, repita la prueba variando el tamaño de la muestra
> pues no logro acomodar el elemento del vector para que genere la muestra
> aleatoria de datos normales.
>
>
> Mx=170
> M0=168
> alpha=0.05
> sigma=8
> T=qt((1-alpha/2),n-1)
> T2=qt(alpha/2,n-1)
> Nr=1000
> Nrechazo=0
> n=seq(10,900, 70)
> dF <- data.frame(PromX = rep(0,length(n)), VarX = rep(0,length(n)))
> for(i in 1:Nr){
> x <- rnorm(n[i], Mx, sigma)
> dF$PromX[i] <- mean(x)
> dF$VarX[i] <- var(x)
> Tobs=(mean(x)-M0)/(sqrt(var(x)/n[i]))
> Rechazo=ifelse(Tobs>=T| Tobs<=T2 ,1,0)
> Nrechazo=Nrechazo+Rechazo
> }
> Nrechazo/Nr
>
> la idea es que cada que haga las nr´repeticiones de la sentencia for,
> genere una nueva muestra normal, con el siguiente valor del vector n.
>
> muchas gracias!
> --
>
> ALEXANDER BERNAL MUÑOZ
> U SANTO TOMAS - ESTADÍSTICA - BOGOTÁ
>
> [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>
--
Saludos,
Carlos Ortega
www.qualityexcellence.es
[[alternative HTML version deleted]]