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

