Foro de debate

Simulación

Re: Simulación

de Luis Benites -
Número de respuestas: 0
Prueba este código, solo agregue un bucle, pero no he conseguido hacer funcionar a la funcion "line" y la función "point"

a<-function(n)
{
par(mfrow=c(4,5),mai=c(0.2,0.2,0.2,0))
for(i in 1:20)
{
poblacion<- seq(50,90,length=1000)
muestra<-sample(poblacion,n,replace=T)
print(muestra)
x=poblacion
et<-sd(poblacion)
yest=dnorm(x,mean(poblacion),et)
plot(x,yest, type="l",col = "gray", lwd=2)
abline(v=mean(poblacion),col="red",lwd=2,lty=2)
l<- mean(sample(poblacion, n,replace=T))
li<- l - qnorm(0.975)*et
ls<- l + qnorm(0.975)*et
c(li,ls)
lines(c(li,ls),dnorm(c(0,0),mean(poblacion),et), col="orange",lwd=2)
points(l,0,col="blue", lwd=1.999,bg="blue",cex=1)
}

}