【问题标题】:Approximation of results结果的近似值
【发布时间】:2014-10-26 06:44:40
【问题描述】:

我只是在学习如何使用 R。我正在练习一些统计数据,如正态分布、泊松等。

当我尝试计算概率并且答案是一个非常接近零 (0) 的数字时,程序显示结果为 0,因此我看不到完整答案,我需要完整答案。总有一个概率,哪怕是很小的概率!!

我的问题是:我可以关闭自逼近或者我可以使用哪个代码来获得完整的答案?

例子:

1-pbinom(q =10, size = 10,prob = 0.8) 

结果:

0

【问题讨论】:

    标签: r approximation


    【解决方案1】:

    pbinom 函数给出了累积密度函数。 i 值小于或等于特定值的概率。因此,对于离散分布,例如具有 10 次绘制的二项式分布

    pbinom(10, 10, .8)
    # [1] 1
    

    告诉您有 100% 的变化,您将观察到 10 次或更少的成功。

    也许您正在考虑概率密度函数(或概率质量函数,因为这是一个离散分布)dbinom

    dbinom(10, 10, .8)
    # [1] 0.1073742
    

    意味着您的所有抽奖都有大约 11% 的机会成功。这也是真的

    sum(dbinom(0:10, 10, .8))
    # [1] 1
    

    0 通过的概率之和正好是 1。

    因此,通过这些案例,您将得到准确的答案。 R 会根据 options(digits=) 设置在控制台中对值进行舍入,但这不是这里发生的事情。

    【讨论】:

      【解决方案2】:

      pbinom 是二项分布的分布函数,它是离散的,因此可以恰好为 1(如您的示例中所示)。您可能一直在考虑 连续 分布,例如正态分布或伽马分布。在这种情况下,四舍五入可能会导致您的结果被截断,例如

      > 1 - pnorm(10, 0, 1)
      [1] 0
      

      但是,p[dist] 函数有一个参数 lower.tail=FALSE 旨在解决这个问题:

      > pnorm(10, 0, 1, lower.tail=FALSE)
      [1] 7.619853e-24
      

      【讨论】:

        猜你喜欢
        • 1970-01-01
        • 1970-01-01
        • 2017-07-31
        • 2019-03-14
        • 2023-03-05
        • 2013-12-08
        • 1970-01-01
        • 1970-01-01
        • 1970-01-01
        相关资源
        最近更新 更多