# 2005 年 06 月 15 日作成
# number of repeated sample
no.rep<-5000
# function of sampling median from population with sample size n
sample.median<-function(n) median(rnorm(n,0,1))
# defining objects to store results
res.5<-c(); res.25<-c(); res.100<-c()
for (i in 1:no.rep) res.5[i]<-sample.median(5)
for (i in 1:no.rep) res.25[i]<-sample.median(25)
for (i in 1:no.rep) res.100[i]<-sample.median(100)
plot(density(res.100),xlim=range(res.5),type="l",main="", 
xlab="sampling distribution of median for n = 5, 25, 100")
lines(density(res.25))
lines(density(res.5))