Displaying 1 result from an estimated 1 matches for "promx".
Did you mean:
promo
2012 Nov 05
1
Potencia de la prueba T
...acer 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 h...