统计计算——随机模拟

star_and_sun

于 2022-12-23 17:21:04 发布

阅读量272

点赞数

CC 4.0 BY-SA版权

分类专栏： R语言文章标签： r语言算法

本文链接：https://blog.csdn.net/star_and_sun/article/details/128421682

R语言专栏收录该内容

9 篇文章

订阅专栏

这篇博客探讨了随机模拟在估计覆盖率方面的应用，使用了Wilson方法来计算置信区间。通过函数`confp.sim1`，模拟了多次实验，计算了置信区间实际包含真实参数的概率，并输出了平均长度、标准差以及覆盖率。该内容涉及统计推断和概率模拟。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

随机模拟
覆盖率：置信区间实际包含真实参数的概率
部分函数如下：

confp.sim1 <- function(n=30, p=0.3, gam=0.95, M=100){
    lam <- qnorm((1+gam)/2)  #分位数
    hatp <- rbinom(M, n, p) / n  #M个p的估计值
    rad1 <- lam * sqrt(hatp*(1-hatp)) / sqrt(n) #λ*根号下p(1-p)
    tildep <- 0.5*lam^2 / (n + lam^2) + n*hatp / (n + lam^2) 
    #上面是Wilson下p的估计
    rad2 <- (lam / sqrt(n) * sqrt( hatp * (1 - hatp) + lam^2/(4*n))/(1 + lam^2/n))
    cov1 <- mean(abs(p - hatp) <= rad1)# 判断p在区间里面的概率
    ncov1 <- M*cov1
    cov2 <- mean(abs(p - tildep) <= rad2)
    ncov2 <- M*cov2
    len1 <- 2*mean(rad1)
    sd1 <- 2*sd(rad1)
    len2 <- 2*mean(rad2)
    sd2 <- 2*sd(rad2)
    cov1
    cov2

}