Con este código (adjunto bibliotecas para quien necesite):
library(rcompanion)
library(psych)
library(MASS)
library(nortest)
library(ggplot2)
n <-1e3
LímiteMenorMuestra=20
LímiteMayorMuestra=150
m <- seq(from=LímiteMenorMuestra, to=LímiteMayorMuestra)
te2 <- for(i in m){
t2err <- 0
for(j in 1:n){
x <- rnorm(i, 0, 1)
if(t.test(x, mu = 0.5)$p.value > 0.05){
t2err <- t2err + 1
}
}
print(i)
print(1-(t2err/n))
}
plot(i, (1-t2err/n))
se obtienen resultados de error tipo 2 para muestras de 20 a 150. Y si bien me imprime cada uno de los valores de Error para cada muestra, tales como
[1] 20
[1] 0.577
[1] 21
[1] 0.561
[1] 22
[1] 0.589
[1] 23
[1] 0.633
A la hora de graficar sólo obtengo el ultimísimo resultado (para i = 150, t2err=1). ¿Alguna forma de hacer que salgan todos y cada uno de los valores de t2err? Mi gratitud por adelantado.