Displaying 1 result from an estimated 1 matches for "nrechazo".
Did you mean:
rechazo
2012 Nov 05
1
Potencia de la prueba T
...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+Rech...