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]]

