更改

删除452字节 、 2021年8月8日 (日) 16:57
第529行: 第529行:     
由于<math> X \sim B(n, p) </math>和<math> Y \sim B(X, q) </math>,由全概率公式,
 
由于<math> X \sim B(n, p) </math>和<math> Y \sim B(X, q) </math>,由全概率公式,
  −
:<math>\Pr[Y = m] &= \sum_{k = m}^{n} \Pr[Y = m \mid X = k] \Pr[X = k] \\[2pt]</math>
      
<math>P(p;\alpha,\beta) =\frac{p^{\alpha-1}(1-p)^{\beta-1}}{\mathrm{B}(\alpha,\beta)}.</math>
 
<math>P(p;\alpha,\beta) =\frac{p^{\alpha-1}(1-p)^{\beta-1}}{\mathrm{B}(\alpha,\beta)}.</math>
第541行: 第539行:     
由于<math>\tbinom{n}{k} \tbinom{k}{m} = \tbinom{n}{m} \tbinom{n-m}{k-m},</math>上式可表示为
 
由于<math>\tbinom{n}{k} \tbinom{k}{m} = \tbinom{n}{m} \tbinom{n-m}{k-m},</math>上式可表示为
  −
:<math> \Pr[Y = m] = \sum_{k=m}^{n} \binom{n}{m} \binom{n-m}{k-m} p^k q^m (1-p)^{n-k} (1-q)^{k-m} </math>
      
将 <math> p^k = p^m p^{k-m} </math> 进行分解,并将所有不依赖于 <math> k </math> 的项从总和中抽出,即可得到
 
将 <math> p^k = p^m p^{k-m} </math> 进行分解,并将所有不依赖于 <math> k </math> 的项从总和中抽出,即可得到
    
<font color="#ff8000">边缘分布 marginal distribution </font>是二项分布较完善的随机数产生方法。
 
<font color="#ff8000">边缘分布 marginal distribution </font>是二项分布较完善的随机数产生方法。
  −
<math>\Pr[Y = m] &= \binom{n}{m} p^m q^m \left( \sum_{k=m}^n \binom{n-m}{k-m} p^{k-m} (1-p)^{n-k} (1-q)^{k-m} \right) \\[2pt]</math>
      
一种从二项分布中产生随机样本的方法是使用<font color="#ff8000">反演算法 inversion algorithm </font>。要做到这一点,我们必须计算从到的所有值的概率。(为了包含整个样本空间,这些概率的和应该接近于1。)然后,通过使用伪随机数生成器来生成介于0和1之间的样本,可以使用在第一步计算出的概率将计算出的样本转换成离散数。
 
一种从二项分布中产生随机样本的方法是使用<font color="#ff8000">反演算法 inversion algorithm </font>。要做到这一点,我们必须计算从到的所有值的概率。(为了包含整个样本空间,这些概率的和应该接近于1。)然后,通过使用伪随机数生成器来生成介于0和1之间的样本,可以使用在第一步计算出的概率将计算出的样本转换成离散数。
第554行: 第548行:     
将 <math> i = k - m </math> 代入上述表达式后,我们得到了
 
将 <math> i = k - m </math> 代入上述表达式后,我们得到了
  −
:<math> \Pr[Y = m] = \binom{n}{m} (pq)^m \left( \sum_{i=0}^{n-m} \binom{n-m}{i} (p - pq)^i (1-p)^{n-m - i} \right) </math>
      
这个分布是由雅各布伯努利 Jacob Bernoulli推导出来的。他考虑了p = r/(r&nbsp;+&nbsp;s)的情形,其中 p 是成功的概率,r 和 s 是正整数。早些时候,布莱斯 · 帕斯卡 Blaise Pascal考虑过p&nbsp;=&nbsp;1/2的情况。
 
这个分布是由雅各布伯努利 Jacob Bernoulli推导出来的。他考虑了p = r/(r&nbsp;+&nbsp;s)的情形,其中 p 是成功的概率,r 和 s 是正整数。早些时候,布莱斯 · 帕斯卡 Blaise Pascal考虑过p&nbsp;=&nbsp;1/2的情况。
863

个编辑