rm(list=ls()) iconf <- function(x, alpha){ ee = qt(1-alpha/2, df = length(x)-1) * sqrt(var(x)/length(x)) return(c(mean(x)-ee, mean(x)+ee)) } Utilizando la función construida anteriormente averiguar gráficamente como el tamaño del intervalo de confianza baja cuando sube el nivel de confianza. Utilizar una muestra de tamaño 40 desde 1000 valores simulados desde una N(5; 1).

nivel <- c(0.1, 0.05, 0.01) 
muestra <- rnorm(1000, 5, 1) 
intervalo <- rep(0, length(nivel)) 
 
for(i in 1:length(nivel)){ 
    for(j in 1:1000){ 
        aux <- sample(muestra, 40) 
        aux <- iconf(aux, nivel[i]) 
        intervalo[i] <- intervalo[i] + diff(aux) 
    } 
    intervalo[i] <- intervalo[i]/1000 
} 
 
plot(nivel, intervalo)