tt=function(x,a) (mean(x)-a)/(sd(x)/sqrt(length(x))) n=16 N=500 átlagh0=100 átlagh1=108 szórás=15 x1=átlagh1+szórás*rnorm(n*N) y1=matrix(x1, nrow=N) hist(y1,breaks=20) mintaátlagok=apply(y1,1,mean) hist(mintaátlagok,breaks=20) mean(mintaátlagok) sd(mintaátlagok) z=apply(y1,1,tt,átlagh0) hist(z, breaks=20) ecdf(z,datadensity="hist",q=c(.25,.5,.75))