数学建模社区-数学中国

标题: 离散函数的数字特征及其R语言的应用 [打印本页]

作者: 杨利霞 时间: 2021-6-24 16:22
标题: 离散函数的数字特征及其R语言的应用

离散函数的数字特征及其R语言的应用
目录
0引言
本文结构
理论公式
1、几何分布
2、负二项分布
3、帕斯卡分布
4、泊松分布
5、参考链接
0引言
本文结构
在文章统计学基础——负二项分布的数字特征1中介绍了负二项分布，在博客2中介绍了离散分布的数字特征。
本文计算一些离散分布的：密度函数、分布函数、均值、方差、偏度、峰度、特征函数、矩母函数

理论公式
为了方便先给出计算公式：

– 密度函数:f ( x ) f(x)f(x)

– 分布函数：F ( x ) = ∫ − ∞ x f ( x ) d x F(x) = \int_{- \infty }^x f(x)dxF(x)=∫
−∞
x

f(x)dx

– 期望：E ( X ) = k 1 E(X) = k_{1}E(X)=k
1

– 方差：D ( X ) = k 2 − k 1 2 D(X) = k_{2}-k_{1}^2D(X)=k
2

−k
1
2

– 特征函数：φ ( t ) = E ( e i t X ) \varphi(t) = E(e^{itX})φ(t)=E(e
itX
)

– 矩母函数：M ( t ) = E ( e t X ) M(t) = E(e^{tX})M(t)=E(e
tX
)

– 中心矩的关系：E ( X k ) = i − k φ ( k ) ( 0 ) = M ( k ) ( 0 ) E(X^k) = i^{-k}\varphi^{(k)}(0) = M^{(k)}(0)E(X
k
)=i
−k
φ
(k)
(0)=M
(k)
(0)

– 偏度：S k e w ( X ) = k 3 k 2 3 / 2 Skew(X) = \frac{k_{3}}{k_{2}^{3/2}}Skew(X)=
k
2
3/2

k
3

3

– 峰度：k u r t ( X ) = k 4 k 2 2 kurt(X) = \frac{k_{4}}{k_{2}^{2}}kurt(X)=
k
2
2

k
4

4

1、几何分布
– 密度函数:f ( x ) = ( 1 − p ) ( x − 1 ) p , f(x) = (1-p)^{(x-1)}p,f(x)=(1−p)
(x−1)
p, x = 1 , 2 , 3 , . . . . . . x = 1,2,3, ... ...x=1,2,3,......

– 分布函数：F ( x ) = ∑ k = 1 x f ( k ) = 1 − ( 1 − p ) x F(x) = \sum_{k=1}^x f(k) = 1 - (1-p)^xF(x)=∑
k=1
x

f(k)=1−(1−p)
x

– 期望：E ( X ) = ∑ k = 1 x k f ( k ) = 1 p E(X) = \sum_{k=1}^x kf(k) = \frac{1}{p}E(X)=∑
k=1
x

kf(k)=
p
1

– 方差：D ( X ) = ∑ k = 1 x k 2 f ( k ) − E ( X ) 2 = 1 − p p 2 D(X) = \sum_{k=1}^x k^2f(k) -E(X)^2= \frac{1-p}{p^2}D(X)=∑
k=1
x

k
2
f(k)−E(X)
2
=
p
2

1−p

– 矩母函数：M ( t ) = p e i t 1 − ( 1 − p ) e i t M(t) = \frac{pe^{it}}{1-(1-p)e^{it}}M(t)=
1−(1−p)e
it

pe
it

– 偏度：S k e w ( X ) = 2 ( 1 − p ) 1 / 2 Skew(X) = 2(1-p)^{1/2}Skew(X)=2(1−p)
1/2

– 峰度：k u r t ( X ) = 9 − 6 p kurt(X) = 9-6pkurt(X)=9−6p

函数功能
dgeom(x, prob, log = FALSE) 概率密度
pgeom(q, prob, lower.tail = TRUE, log.p = FALSE) 累计密度
qgeom(p, prob, lower.tail = TRUE, log.p = FALSE) 分位数
rgeom(n, prob) 随机数
几何分布的各中心距来自5：

2、负二项分布
– 矩母函数：M ( t ) = ( 1 − p ) r ( 1 − p e t ) − r M(t) = (1-p)^r(1-pe^t)^{-r}M(t)=(1−p)
r
(1−pe
t
)
−r

– 偏度：S k e w ( X ) = n 3 + 3 n 2 + 2 n − ( 3 n 2 + 3 n ) p + n p 2 ( n 2 + n ( 1 − p ) ) 3 / 2 Skew(X) = \frac{n^3+3n^2+2n-(3n^2+3n)p+np^2}{(n^2+n(1-p))^{3/2}}Skew(X)=
(n
2
+n(1−p))
3/2

n
3
+3n
2
+2n−(3n
2
+3n)p+np
2

– 峰度：k u r t ( X ) = 略 kurt(X) = 略kurt(X)=略（带入递推公式自行运算）

函数功能
dnbinom(x, size, prob, mu, log = FALSE) 概率密度
pnbinom(q, size, prob, mu, lower.tail = TRUE, log.p = FALSE) 累计密度
qnbinom (p, size, prob, mu, lower.tail = TRUE, log.p = FALSE) 分位数
rnbinom(n, size, prob, mu) 随机数
负二项分布的递推公式如下：6

3、帕斯卡分布
X XX服从r , θ r,\thetar,θ的负二项分布，Y YY服从r , θ r,\thetar,θ的帕斯卡分布。有Y = X + r Y = X +rY=X+r。即：在同样的实验中，帕斯卡分布是成功r rr次后实验（成功+失败）的次数，负二项分布是成功r rr次后失败的次数。负二项分布的数字特征推到见博文1。故帕斯卡分布的数字特征可以由负二项分布推出。
在R语言中我们仍可以使用下面负二项分布的函数做适当调整生成帕斯卡分布。
注：在百度百科7中还有另一种说法是:

帕斯卡分布，负二项分布的正整数形式，描述第n次成功发生在第x次的概率。

我们在课本中见的比较多的是负二项分布的整数形式，这时候课本会标注又称帕斯卡分布，类似于二项分布，负二项分布也可以推广到实数上去定义。大家感兴趣有的可以去自己找资料文献考察一下，这里就不过多的讨论其定义和数字特征的问题了。

函数功能
dnbinom(x, size, prob, mu, log = FALSE) 概率密度
pnbinom(q, size, prob, mu, lower.tail = TRUE, log.p = FALSE) 累计密度
qnbinom (p, size, prob, mu, lower.tail = TRUE, log.p = FALSE) 分位数
rnbinom(n, size, prob, mu) 随机数
4、泊松分布
– 矩母函数：M ( t ) = e λ ( e t − 1 ) M(t) = e^{\lambda(e^t-1)}M(t)=e
λ(e
t
−1)

– 偏度：S k e w ( X ) = λ 3 + 3 λ 2 + λ ( λ 2 + λ ) 3 / 2 Skew(X) = \frac{\lambda^3+3\lambda^2+\lambda}{(\lambda^2+\lambda)^{3/2}}Skew(X)=
(λ
2
+λ)
3/2

λ
3
+3λ
2
+λ

– 峰度：k u r t ( X ) = λ 3 + 6 λ 2 + 7 λ + 1 λ ( λ + 1 ) 2 kurt(X) = \frac{\lambda^3+6\lambda^2+7\lambda+1}{\lambda(\lambda+1)^2}kurt(X)=
λ(λ+1)
2

λ
3
+6λ
2
+7λ+1

函数功能
dpois(x, lambda, log = FALSE) 概率密度
ppois(q, lambda, lower.tail = TRUE, log.p = FALSE) 累计密度
qpois(p, lambda, lower.tail = TRUE, log.p = FALSE) 分位数
rpois(n, lambda) 随机数
中心矩的递推公式来自8：

5、参考链接
https://blog.csdn.net/weixin_46111814/article/details/105243202 ↩︎ ↩︎

https://blog.csdn.net/weixin_44602958/article/details/105261188 ↩︎

https://blog.csdn.net/STcyclone/article/details/84310450 ↩︎

https://baike.baidu.com/item/%E5%81%8F%E5%BA%A6 ↩︎

https://max.book118.com/html/2019/0412/6234220152002022.shtm ↩︎

朱成莲.关于负二项分布高阶矩的教学注记[J].高教学刊,2017,No.66,103-104+107. ↩︎

https://baike.baidu.com/item/%E5%B8%95%E6%96%AF%E5%8D%A1%E5%88%86%E5%B8%83/1188907 ↩︎

https://wenku.baidu.com/view/7f8328c10c22590102029d83.html ↩︎
————————————————
版权声明：本文为CSDN博主「统计学小王子」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
原文链接：https://blog.csdn.net/weixin_46111814/article/details/115499487

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)