<br>
limite.central&lt;-function(r=runif, distpar=c(0,1), m=.5, s=1/sqrt(12), n=c(1,3,10,30), N=1000) {<br>
for(i in n){<br>
if(length(distpar)==2){<br>
x &lt;- matrix(r(i*N, distpar[1], distpar[2]), nc=i)<br>
}else{<br>
x &lt;- matrix(r(i*N, distpar), nc=i)<br>
}<br>
<br>
x&lt;-(apply(x, 1, sum) - i*m )/(sqrt(i)*s)<br>
hist(x, col='light blue', probability=T,main=paste("n=",i),<br>
ylim=c(0,max(.4, density(x)$y)))<br>
lines(density(x), col='red', lwd=3)<br>
curve(dnorm(x), col='blue', lwd=3, lty=3, add=T)<br>
<br>
if( N&gt;100 ){<br>
rug(sample(x,100))<br>
}else{<br>
rug(x)<br>
}<br>
}<br>
}



雷达卡





京公网安备 11010802022788号







