R语言模拟置信区间估计

1.1总体均值的区间估计

 方差已知，大样本

代码：

 attach(faithful) ##获取火山灰数据

population <- sample(eruptions,1000,replace = T) #做1000次有放回取样

N <- length(population) #总体数量

mu <- mean(population) #总体均值

#sd计算样本方差（注：样本方差除以（n-1））

sigma <- sd(population)*sqrt((N-1)/N) #总体标准差

 

layout(matrix(1:1, 1, 1))

experimentes=100

n <- 30

sample_sd <- sigma/sqrt(n)

z<-qnorm(0.975) #计算置信水平为95%的z值

z

 

plot(mu,experimentes,type="n",xlim=c(mu-6*sample_sd,mu+6*sample_sd),

ylim=c(0,experimentes))

abline(v=mu) #x=mu的垂直竖线作为参考线