3.4 分层抽样法

> rm(list = ls())
> set.seed(1)
> n <- 1000
> U <- runif(n)
> h <- function(U){exp(-U)/(1+U^2)}
> theta.e <- mean(h(U))  #积分估计值
> sigma2.e <- sum((h(U)-theta.e)^2)/(n*(n-1))  #积分估计量方差
> theta.e
[1] 0.5251122
> sigma2.e
[1] 5.929583e-05

采用分层等比例抽样方法，模拟量n = 1000，分四层，即k = 4，按照等比例方式，传媒曾抽取样本250个，每层的概率Pi = 0.25.这意味着分层变量Y的分布为：

> k <- 4  #层数
> p <- 1/k  #每层概率
> ni <- n/k  #每层抽取样本量
> theta.i <- c()  #每层样本均值
> sigma2.i <- c()  #每层样本方差
> for(j in 1:k){
+   U <