假设X~Fx,感兴趣的参数为0.则对j=1,·,m: 1.产生第个随机样本，…，日. 2.计算基于第j个样本的置信区间C 3.计算=I(0∈C) 4.计算经验的置信水平可=品∑严1 例5：置信水平的Monte Carlo估计在上一个例子中，我们使用了for循环 来实现计算a的m个置信上界.我们也可以使用replicatei函数： 下Cde n<-20 alpha <-.05 UCL <replicate(1000,expr = x <-rnorm(n,mean 0,sd =2) (n-1)*var(x)/qchisq(alpha,df n-1) ]) Previous Next First Last Back Forward 13bX ∼ FX, a,ÎÍèθ. KÈj = 1, · · · , m: 1. )1jáëÅx (j) 1 , · · · , x (j) n . 2. Oéƒu1jáò&´mCj 3. Oéyj = I(θ ∈ Cj ) 4. Oé²ò&Y²y¯ = 1 m Pm j=1 yj . ~ 5: ò&Y²Monte CarloO 3˛òá~f•, ·Ç¶^ for ÃÇ 5¢yOéσ 2máò&˛.. ·Çè屶^replicateºÍ: ↑Code n <- 20 alpha <- .05 UCL <- replicate(1000, expr = { x <- rnorm(n, mean = 0, sd = 2) (n-1) * var(x) / qchisq(alpha, df = n-1) } ) Previous Next First Last Back Forward 13