Hola, Manuel.
Gracias por atender mi consulta.
Lamento la ambigüedad que se presenta en la misma. Sucede que estaba creando un código para generar las 500 muestras, que el ítem por el cual buscaba apoyo (la ambigüedad se presentó al subir todo el punto); sin embargo, estaba olvidando una parte del ejercicio: la distribución de e; al no ponerla, los histogramas no cuadraban y, por ende, el modelo no se ajustaba a la recta propuesta.
Un compañero me hizo ver el error y, gracias a él, se generó este código que tiene una estructura similar al que buscaba:
x_9 <- seq(1,10)
y_9 <- matrix(data=NA,nrow=10,ncol=500)
###simulaciones
for (i in 1:500) {
y_9[,i] <- 50+(10*x_9)+(rnorm(10,mean = 0, sd=4))
}
#Parte a
estim<- matrix(data=NA,nrow=2,ncol=500)
###estimaciones
for (i in 1:500) {
estim[,i] <- coef(lm(y_9[,i]~x_9))
}
###histogramas
hist(estim[1,], xlab="b_0",main="Valores muestrales de b_0")
hist(estim[2,], xlab="b_1",main="Valores muestrales de b_1")
summary(estim[1,])
summary(estim[2,])
Con esto, ya solucioné el inconveniente que presentaba y te agradezco mucho el que me hayas hecho ver la ambigüedad de la pregunta. Espero que, por favor, me puedas apoyar en cualquier inquietud posterior que esté relacionada con el código.
Saludos.