Cómo escribir programas en lenguaje R
Primero, obtenga la estimación de máxima verosimilitud de la media y la varianza:
mu = media(x)
Varianza = sd(x)^2*(n-1)/n
Crear un conjunto de datos
mu_boot = c()
Crear un bucle
For (i en 1:1000){ p >
y = rnorm(n,mu,sqrt(varianza))
mu_boot[i] = media(y)
Variance_boot[i] = sd(y ) ^2*(n-1)/n
}
Encontrar el intervalo de confianza
Cuantil(mu_boot, c(0.025, 0.975))
Cuantil(varianza_start, c(0,025, 0,975))