概率与统计之五
\color{green}\Large\textbf{概率与统计之五}
概率与统计之五
文章为一些读书笔记,原书为机械工业出版社的《概率论导论》Joseph K.Blizstein和Jessica Hwang著,案例讲解和讲述证明是本书的一大特点,教给人们的不只是一些概率论的概念和公式,而是思考问题的思维方式和策略。不过就是书里的错印有点多,需要边看边改。
_(:з」∠)
文章目录
一. 变换
1. 变量变换公式
关于随机变量的光滑的可逆的变换,可以用变量变换公式处理
- 离 散 情 形 : P ( g ( X ) = y ) = ∑ x : g ( x ) = y P ( X = x ) ⇒ 如 果 y = g ( x ) 为 一 一 映 射 P ( g ( X ) = y ) = P ( X = g − 1 ( y ) ) \textcolor{blue}{离散情形:}\displaystyle P(g(X)=y)=\sum_{x:g(x)=y}P(X=x)\xRightarrow{如果y=g(x)为一一映射} P(g(X)=y)=P(X=g^{-1}(y)) 离散情形:P(g(X)=y)=x:g(x)=y∑P(X=x)如果y=g(x)为一一映射P(g(X)=y)=P(X=g−1(y))
- 连 续 情 形 : 对 于 连 续 情 形 而 言 , 现 在 只 考 虑 g 是 严 格 递 增 或 递 减 的 情 况 , 其 他 情 况 过 于 复 杂 \textcolor{blue}{连续情形:}对于连续情形而言,现在只考虑g是严格递增或递减的情况,其他情况过于复杂 连续情形:对于连续情形而言,现在只考虑g是严格递增或递减的情况,其他情况过于复杂
- 连 续 情 形 , 一 维 空 间 的 変 量 変 换 : 令 X 为 一 随 机 变 量 , 同 时 Y = g ( X ) 推 导 : F g ( x ) ( y ) = P ( g ( X ) ≤ y ) = P ( X ≤ g − 1 ( y ) ) = F X ( g − 1 ( y ) ) ⇒ 等 号 左 右 关 于 y 求 导 f Y ( y ) = f X ( x ) ( g − 1 ( y ) ) ′ = f X ( x ) 1 g ( x ) ′ = f X ( x ) ∣ d x d y ∣ \textcolor{blue}{连续情形,一维空间的変量変换:}令X为一随机变量,同时Y=g(X)\\ 推导:\displaystyle F_{g(x)}(y)=P(g(X)\leq y)=P(X\leq g^{-1}(y))=F_X(g^{-1}(y))\xRightarrow{等号左右关于y求导}f_Y(y)=f_X(x)(g^{-1}(y))^{'}=f_X(x)\frac{1}{g(x)^{'}} = f_X(x)|\frac{dx}{dy}| 连续情形,一维空间的変量変换:令X为一随机变量,同时Y=g(X)推导:Fg(x)(y)=P(g(X)≤y)=P(X≤g−1(y))=FX(g−1(y))等号左右关于y求导fY(y)=fX(x)(g−1(y))′=fX(x)g(x)′1=fX(x)∣dydx∣
- 连 续 情 形 , 多 维 空 间 的 変 量 変 换 : 令 X = ( X 1 , X 2 , … , X n ) 为 一 随 机 向 量 , 同 时 Y = g ( X ) 推 导 : F g ( x ) ( Y ) = F X ( g − 1 ( Y ) ) ⇒ 对 随 机 向 量 Y 求 导 f Y ( y 1 , y 2 , … ) = f X ( x 1 , x 2 , … ) ∣ ∂ ( x 1 , x 2 , … ) ∂ ( y 1 , y 2 , … ) ∣ 关 于 随 机 向 量 的 求 导 为 求 雅 克 比 行 列 式 的 绝 对 值 \textcolor{blue}{连续情形,多维空间的変量変换:}令X=(X_1,X_2,\ldots,X_n)为一随机向量,同时Y=g(X)\\ 推导:\displaystyle F_{g(x)}(Y)=F_X(g^{-1}(Y))\xRightarrow{对随机向量Y求导}f_Y(y_1,y_2,\ldots)=f_X(x_1,x_2,\ldots)|\frac{\partial(x_1,x_2,\ldots)}{\partial(y_1,y_2,\ldots)}|\\关于随机向量的求导为求雅克比行列式的绝对值 连续情形,多维空间的変量変换:令X=(X1,X2,…,Xn)为一随机向量,同时Y=g(X)推导:Fg(x)(Y)=FX(g−1(Y))对随机向量Y求导fY(y1,y2,…)=fX(x1,x2,…)∣∂(y1,y2,…)∂(x1,x2,…)∣关于随机向量的求导为求雅克比行列式的绝对值
- 不 难 看 出 无 论 是 一 维 空 间 还 是 多 维 空 间 的 变 换 , 对 于 连 续 情 形 而 言 , 求 概 率 便 是 求 空 间 中 的 面 积 或 体 积 , 而 变 量 变 换 求 概 率 其 实 就 是 利 用 换 元 法 求 定 积 分 或 多 元 积 分 , 对 积 分 形 式 求 导 变 得 到 了 变 换 后 的 概 率 密 度 函 数 , 另 外 关 于 求 导 或 者 求 雅 克 比 行 列 式 , 分 子 为 Y = g ( X ) 中 的 自 变 量 随 机 向 量 , 分 母 为 因 变 量 随 机 向 量 不难看出无论是一维空间还是多维空间的变换,对于连续情形而言,求概率便是求空间中的面积或体积,\\ 而变量变换求概率其实就是\textcolor{blue}{利用换元法求定积分或多元积分,对积分形式求导变得到了变换后的概率密度函数,}\\ 另外关于求导或者求雅克比行列式,\textcolor{red}{分子为Y=g(X)中的自变量随机向量,分母为因变量随机向量} 不难看出无论是一维空间还是多维空间的变换,对于连续情形而言,求概率便是求空间中的面积或体积,而变量变换求概率其实就是利用换元法求定积分或多元积分,对积分形式求导变得到了变换后的概率密度函数,另外关于求导或者求雅克比行列式,分子为Y=g(X)中的自变量随机向量,分母为因变量随机向量
1). Box-Muller算法
二维情形下的変量変换中有一个比较有名的例子便是
B
o
x
−
M
u
l
l
e
r
Box-Muller
Box−Muller算法,其中关键的便是类似于二元积分中的极坐标换元法
令
U
∼
U
n
i
f
(
0
,
2
π
)
,
T
∼
E
x
p
o
(
1
)
,
且
T
与
U
相
互
独
立
。
定
义
X
=
2
T
c
o
s
U
,
Y
=
2
T
s
i
n
U
。
求
(
X
,
Y
)
的
联
合
概
率
密
度
函
数
。
X
与
Y
是
相
互
独
立
吗
?
其
边
缘
分
布
是
什
么
?
\color{blue}令U\sim Unif(0,2\pi),T\sim Expo(1),且T与U相互独立。定义X=\sqrt{2T}cos\,U,Y=\sqrt{2T}sin\,U。\\求(X,Y)的联合概率密度函数。X与Y是相互独立吗?其边缘分布是什么?
令U∼Unif(0,2π),T∼Expo(1),且T与U相互独立。定义X=2TcosU,Y=2TsinU。求(X,Y)的联合概率密度函数。X与Y是相互独立吗?其边缘分布是什么?
由
于
U
,
T
相
互
独
立
,
所
以
其
联
合
概
率
密
度
函
数
为
f
U
,
T
(
u
,
t
)
=
1
2
π
e
−
t
u
∈
(
0
,
2
π
)
,
t
>
0
很
明
显
由
X
2
+
Y
2
=
2
T
(
c
o
s
2
U
+
s
i
n
2
U
)
=
2
T
可
知
(
2
T
,
U
)
是
(
X
,
Y
)
在
极
坐
标
下
的
表
示
方
法
。
而
且
又
因
为
可
以
从
(
X
,
Y
)
得
到
(
T
,
U
)
,
所
以
这
种
变
换
是
可
逆
的
,
可
以
应
用
便
变
换
公
式
。
首
先
求
雅
可
比
行
列
式
的
绝
对
值
∣
∂
(
u
,
t
)
∂
(
x
,
y
)
∣
=
1
/
∣
∂
(
x
,
y
)
∂
(
u
,
t
)
∣
=
1
/
∣
−
2
t
s
i
n
u
1
2
t
c
o
s
u
2
t
c
o
u
s
u
1
2
t
s
i
n
u
∣
=
1
/
1
=
1
接
下
来
列
出
变
量
变
换
式
f
X
,
Y
(
x
,
y
)
=
f
U
,
T
(
u
,
t
)
⋅
∣
∂
(
u
,
t
)
∂
(
x
,
y
)
∣
=
1
2
π
e
−
t
⋅
1
=
1
2
π
e
−
t
=
1
2
π
e
−
1
2
(
x
2
+
y
2
)
=
1
2
π
e
−
x
2
2
⋅
1
2
π
e
−
y
2
2
很
明
显
f
X
,
Y
可
以
分
解
成
g
(
x
)
h
(
y
)
一
个
x
的
函
数
乘
以
一
个
y
的
函
数
的
形
式
,
所
以
X
,
Y
是
相
互
独
立
的
,
又
∫
−
∞
+
∞
g
(
x
)
d
x
=
1
,
∫
−
∞
+
∞
h
(
y
)
d
y
=
1
,
皆
为
有
效
概
率
密
度
函
数
,
所
以
X
,
Y
的
边
缘
概
率
密
度
函
数
便
是
标
准
正
态
分
布
概
率
密
度
函
数
,
X
,
Y
皆
是
服
从
于
标
准
正
态
分
布
的
随
机
变
量
至
于
有
些
教
材
上
表
示
的
该
算
法
是
由
两
个
标
准
均
匀
分
布
随
机
变
量
U
1
,
U
2
构
成
标
准
正
态
分
布
的
说
法
其
实
也
很
好
理
解
,
无
外
乎
是
将
{
X
=
2
T
c
o
s
U
Y
=
2
T
s
i
n
U
⇒
变
为
{
X
=
−
2
∗
l
n
(
U
1
)
c
o
s
(
2
π
U
2
)
Y
=
−
2
∗
l
n
(
U
1
)
s
i
n
(
2
π
U
2
)
⇒
也
就
是
{
T
=
−
l
n
(
U
1
)
U
=
2
π
U
2
对
于
2
π
U
2
=
U
∼
U
n
i
f
(
0
,
2
π
)
这
很
好
理
解
,
主
要
是
对
于
T
=
−
l
n
(
U
1
)
的
理
解
。
这
是
因
为
根
据
均
匀
分
布
的
普
遍
性
,
将
一
个
随
机
变
量
代
入
其
本
身
的
累
积
分
布
函
数
公
式
,
将
得
到
一
个
标
准
均
匀
分
布
变
量
。
而
T
∼
E
x
p
o
(
1
)
,
服
从
指
数
分
布
,
于
是
一
个
标
准
均
匀
分
布
变
量
Z
=
1
−
e
−
T
⟹
T
=
−
l
n
(
1
−
Z
)
,
而
1
−
Z
同
样
是
一
个
标
准
均
匀
分
布
变
量
,
所
以
设
U
1
=
1
−
Z
,
便
可
以
证
明
由
两
个
标
准
均
匀
分
布
变
量
U
1
,
U
2
可
以
构
成
标
准
正
态
分
布
变
量
的
说
法
\begin{aligned} &由于U,T相互独立,所以其联合概率密度函数为f_{U,T}(u,t)=\frac{1}{2\pi}e^{-t}\quad u\in (0,2\pi),t>0\\ &很明显由X^2+Y^2=2T(cos^2U+sin^2U)=2T可知(\sqrt{2T},U)是(X,Y)在极坐标下的表示方法。\\ &而且又因为可以从(X,Y)得到(T,U),所以这种变换是可逆的,可以应用便变换公式。\\ &\displaystyle 首先求雅可比行列式的绝对值|\frac{\partial(u,t)}{\partial(x,y)}|=1/|\frac{\partial(x,y)}{\partial(u,t)}|=1 / \begin{vmatrix}-\sqrt{2t}sin\,u&\frac{1}{\sqrt{2t}cos\,u}\\\sqrt{2t}cous\,u&\frac{1}{\sqrt{2t}sin\,u}\end{vmatrix}=1/1=1\\ &接下来列出变量变换式\\ &f_{X,Y}(x,y)=f_{U,T}(u,t)\cdot |\frac{\partial(u,t)}{\partial(x,y)}|=\frac{1}{2\pi}e^{-t}\cdot 1=\frac{1}{2\pi}e^{-t}=\frac{1}{2\pi}e^{-\frac{1}{2}(x^2+y^2)}=\frac{1}{\sqrt{2\pi}}e^{\frac{-x^2}{2}}\cdot \frac{1}{\sqrt{2\pi}}e^{\frac{-y^2}{2}}\\ &很明显f_{X,Y}可以分解成g(x)h(y)一个x的函数乘以一个y的函数的形式,所以X,Y是相互独立的,\\ &又\int_{-\infty}^{+\infty}g(x)dx=1,\int_{-\infty}^{+\infty}h(y)dy=1,皆为有效概率密度函数,所以X,Y的边缘概率密度函数便是\\ &标准正态分布概率密度函数,X,Y皆是服从于标准正态分布的随机变量\\ &\color{red}至于有些教材上表示的该算法是由两个标准均匀分布随机变量U_1,U_2构成标准正态分布的说法其实也很好理解,\\ &\color{blue}无外乎是将\begin{cases}X=\sqrt{2T}cos\,U\\Y=\sqrt{2T}sin\,U\end{cases}\xRightarrow{变为}\begin{cases}X=\sqrt{-2 * ln(U_1)}cos(2\pi U_2)\\ Y=\sqrt{-2 * ln(U_1)}sin(2\pi U_2)\end{cases}\xRightarrow{也就是}\begin{cases}T=-ln(U_1)\\ U = 2\pi U_2\end{cases}\\ &对于2\pi U_2 = U \sim Unif(0,2\pi)这很好理解,主要是对于T=-ln(U_1)的理解。这是因为根据均匀分布的普遍性,\\ &将一个随机变量代入其本身的累积分布函数公式,将得到一个标准均匀分布变量。而T\sim Expo(1),服从指数分布,\\ &于是一个标准均匀分布变量Z=1-e^{-T}\implies T= -ln(1-Z),而1-Z同样是一个标准均匀分布变量,\\ &所以设U_1=1-Z,便可以证明由两个标准均匀分布变量U_1,U_2可以构成标准正态分布变量的说法 \end{aligned}
由于U,T相互独立,所以其联合概率密度函数为fU,T(u,t)=2π1e−tu∈(0,2π),t>0很明显由X2+Y2=2T(cos2U+sin2U)=2T可知(2T,U)是(X,Y)在极坐标下的表示方法。而且又因为可以从(X,Y)得到(T,U),所以这种变换是可逆的,可以应用便变换公式。首先求雅可比行列式的绝对值∣∂(x,y)∂(u,t)∣=1/∣∂(u,t)∂(x,y)∣=1/∣∣∣∣∣−2tsinu2tcousu2tcosu12tsinu1∣∣∣∣∣=1/1=1接下来列出变量变换式fX,Y(x,y)=fU,T(u,t)⋅∣∂(x,y)∂(u,t)∣=2π1e−t⋅1=2π1e−t=2π1e−21(x2+y2)=2π1e2−x2⋅2π1e2−y2很明显fX,Y可以分解成g(x)h(y)一个x的函数乘以一个y的函数的形式,所以X,Y是相互独立的,又∫−∞+∞g(x)dx=1,∫−∞+∞h(y)dy=1,皆为有效概率密度函数,所以X,Y的边缘概率密度函数便是标准正态分布概率密度函数,X,Y皆是服从于标准正态分布的随机变量至于有些教材上表示的该算法是由两个标准均匀分布随机变量U1,U2构成标准正态分布的说法其实也很好理解,无外乎是将{X=2TcosUY=2TsinU变为{X=−2∗ln(U1)cos(2πU2)Y=−2∗ln(U1)sin(2πU2)也就是{T=−ln(U1)U=2πU2对于2πU2=U∼Unif(0,2π)这很好理解,主要是对于T=−ln(U1)的理解。这是因为根据均匀分布的普遍性,将一个随机变量代入其本身的累积分布函数公式,将得到一个标准均匀分布变量。而T∼Expo(1),服从指数分布,于是一个标准均匀分布变量Z=1−e−T⟹T=−ln(1−Z),而1−Z同样是一个标准均匀分布变量,所以设U1=1−Z,便可以证明由两个标准均匀分布变量U1,U2可以构成标准正态分布变量的说法
2. 卷积
在这里卷积代表的是独立随机变量的和。也就是设两个相互独立的随机变量 X , Y X,Y X,Y来求他们的和变量 T = X + Y T=X+Y T=X+Y的分布情况。其实就是利用全概率公式。罗列出所有 X i + Y j = T X_i+Y_j=T Xi+Yj=T的事件,并把每一种事件的概率 P ( X i ) P ( Y j = T − X i ) P(X_i)P(Y_j=T-X_i) P(Xi)P(Yj=T−Xi)全部加起来即可,由于 X i + Y j X_i+Y_j Xi+Yj的值始终为 T T T所以才称之为卷积。只不过离散跟连续的情形不一样,一种是概率质量函数,一种是概率密度函数。
- 离 散 情 况 : P ( T = t ) = ∑ x P ( Y = t − x ) P ( X = x ) = ∑ y P ( X = t − y ) P ( Y = y ) \displaystyle\textcolor{blue}{离散情况:}P(T=t)=\sum_xP(Y=t-x)P(X=x)=\sum_yP(X=t-y)P(Y=y) 离散情况:P(T=t)=x∑P(Y=t−x)P(X=x)=y∑P(X=t−y)P(Y=y)
- 连 续 情 况 : F T ( t ) = P ( X + Y ≤ t ) = ∫ − ∞ + ∞ P ( Y ≤ t − x ) f X ( x ) d x = ∫ − ∞ + ∞ F Y ( t − x ) f X ( x ) d x 根 据 含 参 变 量 额 定 积 分 , 积 分 变 量 为 参 变 量 时 积 分 与 导 数 顺 序 可 以 互 换 的 原 则 F T ( t ) = ∫ − ∞ + ∞ F Y ( t − x ) f X ( x ) d x ⇒ 左 右 对 t 求 导 f T ( t ) = ∫ − ∞ + ∞ f Y ( t − x ) f X ( x ) d x \displaystyle\textcolor{blue}{连续情况:}F_T(t)=P(X+Y \leq t)=\int_{-\infty}^{+\infty}P(Y\leq t-x)f_X(x)dx=\int_{-\infty}^{+\infty}F_Y(t-x)f_X(x)dx\\ 根据含参变量额定积分,积分变量为参变量时积分与导数顺序可以互换的原则\\F_T(t)=\int_{-\infty}^{+\infty}F_Y(t-x)f_X(x)dx\xRightarrow{左右对t求导}f_T(t)=\int_{-\infty}^{+\infty}f_Y(t-x)f_X(x)dx 连续情况:FT(t)=P(X+Y≤t)=∫−∞+∞P(Y≤t−x)fX(x)dx=∫−∞+∞FY(t−x)fX(x)dx根据含参变量额定积分,积分变量为参变量时积分与导数顺序可以互换的原则FT(t)=∫−∞+∞FY(t−x)fX(x)dx左右对t求导fT(t)=∫−∞+∞fY(t−x)fX(x)dx
3. 贝塔分布
贝塔分布是定义在
(
0
,
1
)
(0,1)
(0,1)区间上的连续分布,他是均匀分布
U
n
i
f
(
0
,
1
)
Unif(0,1)
Unif(0,1)的推广形式,允许其概率密度函数在
(
0
,
1
)
(0,1)
(0,1)上非恒定。表达式为
X
∼
B
e
t
a
(
a
,
b
)
X\sim Beta(a,b)
X∼Beta(a,b)
相关性质:
- 概 率 密 度 函 数 : f ( x ) = 1 β ( a , b ) x a − 1 ( 1 − x ) b − 1 , 0 < x < 1 , a > 0 , b > 0 ∫ 0 1 1 β ( a , b ) x a − 1 ( 1 − x ) b − 1 d x = 1 \displaystyle \textcolor{blue}{概率密度函数:}f(x)=\frac{1}{\beta(a,b)}x^{a-1}(1-x)^{b-1},0<x<1,a>0,b>0\qquad \int_0^1\frac{1}{\beta(a,b)}x^{a-1}(1-x)^{b-1}dx=1 概率密度函数:f(x)=β(a,b)1xa−1(1−x)b−1,0<x<1,a>0,b>0∫01β(a,b)1xa−1(1−x)b−1dx=1
- 贝 塔 分 布 的 归 一 化 常 数 : { ∫ 0 1 ( n k ) x k ( 1 − x ) n − k d x = 1 n + 1 ∫ 0 1 x a − 1 ( 1 − x ) b − 1 d x = β ( a , b ) ⇒ 将 k = a − 1 , n = a + b − 2 代 入 上 式 , 并 联 立 方 程 组 可 求 得 β ( a , b ) = ( a − 1 ) ! ( b − 1 ) ! ( a + b − 1 ) ! 于 是 可 得 贝 塔 分 布 的 归 一 化 常 数 为 1 β ( a , b ) = ( a + b − 1 ) ! ( a − 1 ) ! ( b − 1 ) ! , 另 外 该 常 数 的 另 一 种 形 式 为 1 β ( a , b ) = Γ ( a + b ) Γ ( a ) Γ ( b ) , 这 是 涉 及 到 使 用 两 个 独 立 的 伽 马 变 量 来 表 示 一 个 贝 塔 变 量 时 使 用 \textcolor{blue}{贝塔分布的归一化常数:}\\ \begin{cases} \displaystyle \int_0^1{n \choose k}x^k(1-x)^{n-k}dx&\displaystyle=\frac{1}{n+1}\\ \displaystyle\int_0^1x^{a-1}(1-x)^{b-1}dx&\displaystyle=\beta(a,b)\end{cases} \xRightarrow{将k=a-1\ ,n=a+b-2代入上式,并联立方程组可求得}\displaystyle\beta(a,b)=\frac{(a-1)!(b-1)!}{(a+b-1)!}\\ 于是可得贝塔分布的归一化常数为\frac{1}{\beta(a,b)}=\frac{(a+b-1)!}{(a-1)!(b-1)!},另外该常数的另一种形式为\\ \frac{1}{\beta(a,b)}=\frac{\Gamma(a+b)}{\Gamma(a)\Gamma(b)},这是涉及到使用两个独立的伽马变量来表示一个贝塔变量时使用 贝塔分布的归一化常数:⎩⎪⎪⎨⎪⎪⎧∫01(kn)xk(1−x)n−kdx∫01xa−1(1−x)b−1dx=n+11=β(a,b)将k=a−1 ,n=a+b−2代入上式,并联立方程组可求得β(a,b)=(a+b−1)!(a−1)!(b−1)!于是可得贝塔分布的归一化常数为β(a,b)1=(a−1)!(b−1)!(a+b−1)!,另外该常数的另一种形式为β(a,b)1=Γ(a)Γ(b)Γ(a+b),这是涉及到使用两个独立的伽马变量来表示一个贝塔变量时使用
- 期 望 与 方 差 : E ( B ) = a a + b , E ( B 2 ) = a ( a + 1 ) ( a + b ) ( a + b + 1 ) , E ( B k ) = a ( a + 1 ) ⋯ ( a + k − 1 ) ( a + b ) ( a + b + 1 ) ⋯ ( a + b + k − 1 ) , V a r ( B ) = a b ( a + b ) 2 ( a + b + 1 ) \textcolor{blue}{期望与方差:}\displaystyle E(B)=\frac{a}{a+b},E(B^2)=\frac{a(a+1)}{(a+b)(a+b+1)},E(B^k)=\frac{a(a+1)\cdots(a+k-1)}{(a+b)(a+b+1)\cdots(a+b+k-1)},\\Var(B)=\frac{ab}{(a+b)^2(a+b+1)} 期望与方差:E(B)=a+ba,E(B2)=(a+b)(a+b+1)a(a+1),E(Bk)=(a+b)(a+b+1)⋯(a+b+k−1)a(a+1)⋯(a+k−1),Var(B)=(a+b)2(a+b+1)ab
-
概
率
密
度
曲
线
特
征
:
\textcolor{blue}{概率密度曲线特征:}
概率密度曲线特征:
a = 1 ∣ b = 1 , 分 布 曲 线 为 斜 率 不 为 零 的 直 线 , 若 a = b = 1 则 B e t a ( 1 , 1 ) = U n i f ( 0 , 1 ) a < 1 , b < 1 , 曲 线 为 U 型 开 口 朝 上 ; a > 1 , b > 1 曲 线 开 口 朝 下 a = b ≠ 1 , 则 曲 线 关 于 1 / 2 对 称 , a > b , 向 大 于 1 / 2 方 向 偏 移 , a < b , 向 小 于 1 / 2 方 向 偏 移 \begin{aligned} &a=1|b=1,分布曲线为斜率不为零的直线,若a=b=1则Beta(1,1)=Unif(0,1)\\ &a<1,b<1,曲线为U型开口朝上;a>1,b>1曲线开口朝下\\ &a=b\neq1,则曲线关于1/2对称,a > b,向大于1/2方向偏移,a< b,向小于1/2方向偏移 \end{aligned} a=1∣b=1,分布曲线为斜率不为零的直线,若a=b=1则Beta(1,1)=Unif(0,1)a<1,b<1,曲线为U型开口朝上;a>1,b>1曲线开口朝下a=b=1,则曲线关于1/2对称,a>b,向大于1/2方向偏移,a<b,向小于1/2方向偏移 - 案 例 : 贝 塔 分 布 的 案 例 有 很 多 , 一 个 常 用 的 案 例 就 是 用 一 个 贝 塔 分 布 变 量 来 表 示 未 知 概 率 , 也 就 是 说 , 可 以 用 贝 塔 分 布 给 未 知 概 率 定 义 概 率 。 由 此 引 出 贝 塔 分 布 与 二 项 分 布 的 共 轭 性 \textcolor{blue}{案例:}贝塔分布的案例有很多,一个常用的案例就是用一个贝塔分布变量来表示未知概率,也就是说,\\可以用贝塔分布给未知概率定义概率。由此\textcolor{blue}{引出贝塔分布与二项分布的共轭性} 案例:贝塔分布的案例有很多,一个常用的案例就是用一个贝塔分布变量来表示未知概率,也就是说,可以用贝塔分布给未知概率定义概率。由此引出贝塔分布与二项分布的共轭性
- 贝 塔 分 布 与 二 项 分 布 的 共 轭 性 : 基 本 逻 辑 假 设 要 推 断 一 个 试 验 每 次 测 试 的 成 功 率 p 。 当 然 最 简 单 的 方 法 就 是 做 大 量 测 试 n , 统 计 其 成 功 次 数 k 和 失 败 次 数 n − k , 便 可 以 大 体 估 算 出 该 试 验 每 一 次 的 成 功 率 。 n 越 大 , 则 估 计 的 p 就 越 接 近 真 实 。 那 么 这 个 过 程 的 能 否 进 行 数 学 抽 象 化 的 实 现 , 一 般 常 用 的 办 法 是 利 用 贝 叶 斯 准 则 的 贝 叶 斯 推 断 。 把 所 有 未 知 量 都 视 作 随 机 变 量 , 再 用 贝 叶 斯 准 则 推 出 p 的 随 机 分 布 。 既 然 使 用 到 了 贝 叶 斯 准 则 , 那 么 就 需 要 给 每 一 个 变 量 赋 予 一 个 随 机 分 布 而 使 之 成 为 一 个 随 机 变 量 , 对 于 二 项 分 布 试 验 的 成 功 次 数 k , 很 显 然 可 以 定 义 它 是 一 个 服 从 二 项 分 布 的 随 机 变 量 X ∼ B i n ( n , p ) , 现 在 的 问 题 在 于 试 验 成 功 概 率 p 是 一 个 服 从 什 么 分 布 的 随 机 变 量 。 根 据 已 经 论 证 的 结 果 p 服 从 的 是 贝 塔 分 布 。 分 布 说 明 一 假 设 一 个 概 率 值 p 为 一 随 机 变 量 , 其 先 验 分 布 为 贝 塔 分 布 , 即 p ∼ B e t a ( a , b ) , 二 以 p 为 一 次 伯 努 利 试 验 的 实 际 成 功 率 , n 为 总 试 验 次 数 , 于 是 在 该 条 件 下 成 功 次 数 的 随 机 变 量 为 X ∣ p ∼ B i n ( n , p ) , 服 从 二 项 分 布 。 三 现 在 已 经 知 道 了 一 个 边 缘 分 布 p ∼ B e t a ( a , b ) 和 一 个 条 件 分 布 X ∣ p ∼ B i n ( n , p ) , 那 么 便 可 以 构 造 一 个 联 合 分 布 : f X , p ( x , p ) = P ( X = k ∣ p ) f ( p ) , 这 是 一 个 混 合 分 布 并 可 以 求 出 X 的 边 缘 概 率 质 量 函 数 P ( X = k ) = ∫ 0 1 P ( X = k ∣ p ) f ( p ) d p 以 及 在 n 次 试 验 后 , 成 功 次 数 为 k 的 概 率 P ( X = k ∣ p ) = ( n k ) p k ( 1 − p ) n − k 四 在 弄 明 白 了 所 有 条 件 后 , 便 可 利 用 贝 叶 斯 准 则 求 出 q 的 条 件 概 率 密 度 函 数 , f ( p ∣ X = k ) = P ( X = k ∣ p ) f ( p ) P ( X = k ) ⇒ f ( p ∣ X = k ) ∝ p a + k − 1 ( 1 − p ) b + n − k − 1 ⇒ p ∣ X = k ∼ B e t a ( a + k , b + n − k ) 也 就 是 说 p 的 后 验 分 布 同 样 服 从 于 贝 塔 分 布 。 由 于 p 的 先 验 后 验 分 布 全 部 服 从 贝 塔 分 布 , 且 在 给 定 p 的 条 件 下 数 据 服 从 二 项 分 布 , 于 是 称 贝 塔 分 布 为 二 项 分 布 的 共 轭 先 验 分 布 五 关 于 B e t a ( a , b ) 和 B e t a ( a + k , b + n − k ) 中 参 数 的 解 释 , a 为 事 先 观 察 的 成 功 试 验 次 数 , b 为 事 先 观 察 的 失 败 试 验 次 数 , 而 在 进 行 了 n 次 试 验 后 , 参 数 变 为 a + k , 即 试 验 的 累 积 成 功 次 数 , 以 及 b + n − k , 即 试 验 的 累 积 失 败 次 数 。 六 使 用 贝 塔 分 布 的 意 义 。 一 般 而 言 随 着 我 们 试 验 次 数 增 多 , 收 集 到 的 数 据 积 累 , 会 慢 慢 更 新 信 息 : 即 从 先 验 分 布 计 算 得 到 后 验 分 布 , 然 后 用 后 验 分 布 替 代 旧 的 先 验 分 布 从 而 得 到 新 的 先 验 分 布 , 再 用 新 的 先 验 分 布 计 算 得 到 新 的 后 验 分 布 , 如 此 循 环 往 复 。 通 过 上 面 的 证 明 我 们 可 以 得 知 , 对 于 一 个 二 项 分 布 的 成 功 概 率 p , 它 的 这 些 更 新 信 息 的 循 环 过 程 是 被 包 含 在 贝 塔 分 布 族 中 进 行 的 。 我 们 只 需 要 通 过 试 验 不 断 更 新 贝 塔 分 布 中 的 两 个 参 数 , 即 试 验 的 历 史 累 积 成 功 次 数 和 试 验 的 历 史 累 积 失 败 次 数 , 便 可 以 从 贝 塔 分 布 中 生 成 一 个 符 合 试 验 结 果 的 成 功 概 率 p \textcolor{blue}{贝塔分布与二项分布的共轭性}:\\ \textcolor{blue}{基本逻辑}\\ 假设要推断一个试验每次测试的成功率p。当然最简单的方法就是做大量测试n,统计其成功次数k和\\ 失败次数n-k,便可以大体估算出该试验每一次的成功率。n越大,则估计的p就越接近真实。\\ 那么这个过程的能否进行数学抽象化的实现,一般常用的办法是利用贝叶斯准则的贝叶斯推断。\\ 把所有未知量都视作随机变量,再用贝叶斯准则推出p的随机分布。既然使用到了贝叶斯准则,\\ 那么就需要给每一个变量赋予一个随机分布而使之成为一个随机变量,\\ 对于二项分布试验的成功次数k,很显然可以定义它是一个服从二项分布的随机变量X\sim Bin(n,p),\\ 现在的问题在于试验成功概率p是一个服从什么分布的随机变量。根据已经论证的结果p服从的是贝塔分布。\\ \textcolor{blue}{分布说明}\\ \textcolor{blue}{一}\quad假设一个概率值p为一随机变量,其先验分布为贝塔分布,即p\sim Beta(a,b),\\ \textcolor{blue}{二}\quad以p为一次伯努利试验的实际成功率,n为总试验次数,于是在该条件下成功次数的随机变量为\\ \qquad X|p\sim Bin(n,p),服从二项分布。\\ \textcolor{blue}{三}\quad现在已经知道了一个边缘分布p\sim Beta(a,b)和一个条件分布X|p\sim Bin(n,p),\\ \qquad\displaystyle那么便可以构造一个联合分布:f_X,p(x,p)=P(X=k|p)f(p),这是一个混合分布\\ \qquad\displaystyle并可以求出X的边缘概率质量函数P(X=k)=\int_0^1P(X=k|p)f(p)dp\\ \qquad\displaystyle以及在n次试验后,成功次数为k的概率P(X=k|p)={n \choose k}p^k(1-p)^{n-k}\\ \textcolor{blue}{四}\quad\displaystyle 在弄明白了所有条件后,便可利用贝叶斯准则求出q的条件概率密度函数,\\ \quad \displaystyle f(p|X=k)=\frac{P(X=k|p)f(p)}{P(X=k)}\rArr f(p|X=k) \propto p^{a+k-1}(1-p)^{b+n-k-1}\rArr p|X=k \sim Beta(a+k,b+n-k)\\ \qquad\displaystyle 也就是说p的后验分布同样服从于贝塔分布。由于p的先验后验分布全部服从贝塔分布,\\ \qquad\displaystyle 且在给定p的条件下数据服从二项分布,于是称贝塔分布为二项分布的共轭先验分布\\ \textcolor{blue}{五}\quad关于Beta(a,b)和Beta(a+k,b+n-k)中参数的解释,\\ \qquad\displaystyle a为事先观察的成功试验次数,b为事先观察的失败试验次数,而在进行了n次试验后,\\ \qquad\displaystyle 参数变为a+k,即试验的累积成功次数,以及b+n-k,即试验的累积失败次数。\\ \textcolor{blue}{六}\quad使用贝塔分布的意义。一般而言随着我们试验次数增多,收集到的数据积累,会慢慢更新信息:\\ \qquad\displaystyle 即\textcolor{blue}{从先验分布计算得到后验分布,然后用后验分布替代旧的先验分布从而得到新的先验分布,\\ \qquad\displaystyle 再用新的先验分布计算得到新的后验分布,如此循环往复}。\\ \qquad\displaystyle 通过上面的证明我们可以得知,对于一个二项分布的成功概率p,它的这些更新信息的循环过程\\ \qquad\displaystyle是被包含在贝塔分布族中进行的。我们只需要通过试验不断更新贝塔分布中的两个参数,\\ \qquad\displaystyle 即试验的历史累积成功次数和试验的历史累积失败次数,便可以从贝塔分布中生成一个符合试验结果的成功概率p 贝塔分布与二项分布的共轭性:基本逻辑假设要推断一个试验每次测试的成功率p。当然最简单的方法就是做大量测试n,统计其成功次数k和失败次数n−k,便可以大体估算出该试验每一次的成功率。n越大,则估计的p就越接近真实。那么这个过程的能否进行数学抽象化的实现,一般常用的办法是利用贝叶斯准则的贝叶斯推断。把所有未知量都视作随机变量,再用贝叶斯准则推出p的随机分布。既然使用到了贝叶斯准则,那么就需要给每一个变量赋予一个随机分布而使之成为一个随机变量,对于二项分布试验的成功次数k,很显然可以定义它是一个服从二项分布的随机变量X∼Bin(n,p),现在的问题在于试验成功概率p是一个服从什么分布的随机变量。根据已经论证的结果p服从的是贝塔分布。分布说明一假设一个概率值p为一随机变量,其先验分布为贝塔分布,即p∼Beta(a,b),二以p为一次伯努利试验的实际成功率,n为总试验次数,于是在该条件下成功次数的随机变量为X∣p∼Bin(n,p),服从二项分布。三现在已经知道了一个边缘分布p∼Beta(a,b)和一个条件分布X∣p∼Bin(n,p),那么便可以构造一个联合分布:fX,p(x,p)=P(X=k∣p)f(p),这是一个混合分布并可以求出X的边缘概率质量函数P(X=k)=∫01P(X=k∣p)f(p)dp以及在n次试验后,成功次数为k的概率P(X=k∣p)=(kn)pk(1−p)n−k四在弄明白了所有条件后,便可利用贝叶斯准则求出q的条件概率密度函数,f(p∣X=k)=P(X=k)P(X=k∣p)f(p)⇒f(p∣X=k)∝pa+k−1(1−p)b+n−k−1⇒p∣X=k∼Beta(a+k,b+n−k)也就是说p的后验分布同样服从于贝塔分布。由于p的先验后验分布全部服从贝塔分布,且在给定p的条件下数据服从二项分布,于是称贝塔分布为二项分布的共轭先验分布五关于Beta(a,b)和Beta(a+k,b+n−k)中参数的解释,a为事先观察的成功试验次数,b为事先观察的失败试验次数,而在进行了n次试验后,参数变为a+k,即试验的累积成功次数,以及b+n−k,即试验的累积失败次数。六使用贝塔分布的意义。一般而言随着我们试验次数增多,收集到的数据积累,会慢慢更新信息:即从先验分布计算得到后验分布,然后用后验分布替代旧的先验分布从而得到新的先验分布,再用新的先验分布计算得到新的后验分布,如此循环往复。通过上面的证明我们可以得知,对于一个二项分布的成功概率p,它的这些更新信息的循环过程是被包含在贝塔分布族中进行的。我们只需要通过试验不断更新贝塔分布中的两个参数,即试验的历史累积成功次数和试验的历史累积失败次数,便可以从贝塔分布中生成一个符合试验结果的成功概率p
4. 伽马分布
伽马分布是定义在正实数轴上的一个连续分布;它是指数分布的推广。若一个指数随机变量表示在无记忆条件下首次成功的等待时间,则伽马分布的随机变量表示多次成功的总等待时间。也就是说它同样是可以通过泊松过程来得到。至于伽马分布的概率密度表达式部分是基于伽马函数
Γ
(
a
)
=
∫
0
+
∞
x
a
e
−
x
d
x
x
,
a
>
0
\displaystyle \Gamma(a)=\int_0^{+\infty}x^ae^{-x}\frac{dx}{x},a>0
Γ(a)=∫0+∞xae−xxdx,a>0而构成的。表达式为
X
∼
G
a
m
m
a
(
a
,
λ
)
X\sim Gamma(a,\lambda)
X∼Gamma(a,λ)
相关性质为
- 概 率 密 度 函 数 : f ( y ) = 1 Γ ( a ) ( λ y ) a e − λ y 1 y y > 0 , a > 0 , λ > 0 ∫ 0 + ∞ 1 Γ ( a ) ( λ y ) a e − λ y 1 y d y = 1 a 可 以 看 做 是 等 待 成 功 的 次 数 , λ 则 与 泊 松 分 布 , 指 数 分 布 中 的 定 义 一 致 , 为 成 功 速 率 。 当 然 这 不 是 绝 对 , 在 共 轭 案 例 中 又 有 别 的 定 义 \displaystyle \textcolor{blue}{概率密度函数:} f(y)=\frac{1}{\Gamma(a)}(\lambda y)^ae^{-\lambda y}\frac{1}{y}\quad y>0,a>0, \lambda > 0\quad \int_0^{+\infty}\frac{1}{\Gamma(a)}(\lambda y)^ae^{-\lambda y}\frac{1}{y}dy=1 \\a可以看做是等待成功的次数,\lambda则与泊松分布,指数分布中的定义一致,为成功速率。\\ 当然这不是绝对,在共轭案例中又有别的定义 概率密度函数:f(y)=Γ(a)1(λy)ae−λyy1y>0,a>0,λ>0∫0+∞Γ(a)1(λy)ae−λyy1dy=1a可以看做是等待成功的次数,λ则与泊松分布,指数分布中的定义一致,为成功速率。当然这不是绝对,在共轭案例中又有别的定义
- 伽 马 函 数 的 几 个 重 要 性 质 : 一 Γ ( s ) = ∫ 0 + ∞ e − x x s − 1 d x , s > 0 二 Γ ( a + 1 ) = a Γ ( a ) . 对 所 有 a > 0 成 立 三 Γ ( n ) = ( n − 1 ) ! , n 为 正 整 数 四 Γ ( 1 2 ) = π Γ ( 1 ) = 1 \displaystyle \textcolor{blue}{伽马函数的几个重要性质:}\\ 一\quad\Gamma(s)=\int_0^{+\infty}e^{-x}x^{s-1}dx,s>0\\ 二\quad\Gamma(a+1)=a\Gamma(a).对所有a>0成立\\ 三\quad\Gamma(n)=(n-1)!,n为正整数\\ 四\quad\Gamma(\frac{1}{2})=\sqrt{\pi}\quad \Gamma(1)=1 伽马函数的几个重要性质:一Γ(s)=∫0+∞e−xxs−1dx,s>0二Γ(a+1)=aΓ(a).对所有a>0成立三Γ(n)=(n−1)!,n为正整数四Γ(21)=πΓ(1)=1
- 期 望 与 方 差 : 若 X ∼ G a m m a ( a , 1 ) E ( X ) = Γ ( a + 1 ) Γ ( a ) = a , E ( X 2 ) = Γ ( a + 2 ) Γ ( a ) = ( a + 1 ) a , E ( X n ) = Γ ( a + n ) Γ ( a ) = a ( a + 1 ) ⋯ ( a + n − 1 ) V a r ( X ) = a 若 Y ∼ G a m m a ( a , λ ) , E ( Y ) = a λ , E ( Y n ) = 1 λ n ⋅ Γ ( a + n ) Γ ( a ) , V a r ( Y ) = a λ 2 \textcolor{blue}{期望与方差:}若X\sim Gamma(a,1)\\ \displaystyle E(X)=\frac{\Gamma(a+1)}{\Gamma(a)}=a,E(X^2)=\frac{\Gamma(a+2)}{\Gamma(a)}=(a+1)a,E(X^n)=\frac{\Gamma(a+n)}{\Gamma(a)}=a(a+1)\cdots(a+n-1)\\ Var(X)=a\\ 若Y\sim Gamma(a,\lambda),E(Y)=\frac{a}{\lambda},E(Y^n)=\frac{1}{\lambda^n}\cdot\frac{\Gamma(a+n)}{\Gamma(a)},Var(Y)=\frac{a}{\lambda^2} 期望与方差:若X∼Gamma(a,1)E(X)=Γ(a)Γ(a+1)=a,E(X2)=Γ(a)Γ(a+2)=(a+1)a,E(Xn)=Γ(a)Γ(a+n)=a(a+1)⋯(a+n−1)Var(X)=a若Y∼Gamma(a,λ),E(Y)=λa,E(Yn)=λn1⋅Γ(a)Γ(a+n),Var(Y)=λ2a
- 指 数 分 布 与 伽 马 分 布 : 令 X 1 , … , X n 独 立 同 分 布 与 指 数 分 布 E x p o ( λ ) , 则 X 1 + ⋯ + X n ∼ G a m m a ( n , λ ) 简 单 点 说 就 是 几 个 独 立 同 分 布 的 指 数 变 量 的 和 服 从 伽 马 分 布 。 如 果 n = 1 的 话 , G a m m a ( 1 , λ ) 即 为 指 数 分 布 。 直 观 来 讲 , 伽 马 分 布 G a m m a ( n , λ ) 为 第 n 次 试 验 成 功 时 所 需 要 等 待 的 时 间 变 量 的 分 布 , 而 两 次 成 功 试 验 的 间 隔 时 间 变 量 就 是 一 个 指 数 分 布 , 由 于 指 数 分 布 的 无 记 忆 性 , 这 些 间 隔 变 量 都 是 独 立 同 分 布 的 为 E x p o ( λ ) , 很 明 显 把 这 些 时 间 间 隔 都 加 起 来 就 是 总 共 等 待 的 时 间 , 所 以 独 立 同 分 布 的 指 数 变 量 的 和 , 服 从 伽 马 分 布 \textcolor{blue}{指数分布与伽马分布:}令X_1,\ldots,X_n独立同分布与指数分布Expo(\lambda),则X_1+\cdots+X_n\sim Gamma(n,\lambda)\\ 简单点说就是几个独立同分布的指数变量的和服从伽马分布。如果n=1的话,Gamma(1,\lambda)即为指数分布。\\ 直观来讲,伽马分布Gamma(n,\lambda)为第n次试验成功时所需要等待的时间变量的分布,而两次成功试验的间隔\\ 时间变量就是一个指数分布,由于指数分布的无记忆性,这些间隔变量都是独立同分布的为Expo(\lambda),\\ 很明显把这些时间间隔都加起来就是总共等待的时间,所以独立同分布的指数变量的和,服从伽马分布 指数分布与伽马分布:令X1,…,Xn独立同分布与指数分布Expo(λ),则X1+⋯+Xn∼Gamma(n,λ)简单点说就是几个独立同分布的指数变量的和服从伽马分布。如果n=1的话,Gamma(1,λ)即为指数分布。直观来讲,伽马分布Gamma(n,λ)为第n次试验成功时所需要等待的时间变量的分布,而两次成功试验的间隔时间变量就是一个指数分布,由于指数分布的无记忆性,这些间隔变量都是独立同分布的为Expo(λ),很明显把这些时间间隔都加起来就是总共等待的时间,所以独立同分布的指数变量的和,服从伽马分布
- 伽 马 变 量 的 和 : 若 X ∼ G a m m a ( a , λ ) , Y ∼ G a m m a ( b , λ ) 相 互 独 立 , 其 中 a 和 b 是 整 数 , 则 X + Y ∼ G a m m a ( a + b , λ ) 。 这 个 由 第 四 条 很 容 易 证 明 , X , Y 可 以 表 示 为 一 系 列 相 同 独 立 同 分 布 的 指 数 变 量 的 和 , 而 所 有 变 量 的 λ 相 同 , X 为 a 个 E x p o ( λ ) 相 加 , Y 为 b 个 E x p o ( λ ) 相 加 , 所 以 两 拨 指 数 变 量 的 和 仍 然 服 从 伽 马 分 布 , 并 且 参 数 为 a + b 和 λ \textcolor{blue}{伽马变量的和:}\\ 若X\sim Gamma(a,\lambda),Y\sim Gamma(b,\lambda)相互独立,其中a和b是整数,则X+Y\sim Gamma(a+b,\lambda)。\\ 这个由第四条很容易证明,X,Y可以表示为一系列相同独立同分布的指数变量的和,而所有变量的\lambda相同,\\ X为a个Expo(\lambda)相加,Y为b个Expo(\lambda)相加,所以两拨指数变量的和仍然服从伽马分布,并且参数为a+b和\lambda 伽马变量的和:若X∼Gamma(a,λ),Y∼Gamma(b,λ)相互独立,其中a和b是整数,则X+Y∼Gamma(a+b,λ)。这个由第四条很容易证明,X,Y可以表示为一系列相同独立同分布的指数变量的和,而所有变量的λ相同,X为a个Expo(λ)相加,Y为b个Expo(λ)相加,所以两拨指数变量的和仍然服从伽马分布,并且参数为a+b和λ
- 伽 马 分 布 与 负 二 项 分 布 : 类 似 于 指 数 分 布 和 几 何 分 布 。 伽 马 分 布 变 量 是 完 成 n 次 成 功 试 验 时 所 等 待 的 时 间 , 而 负 二 项 分 布 变 量 是 完 成 n 次 成 功 试 验 时 所 经 历 的 所 有 试 验 次 数 , 如 果 每 次 试 验 测 试 的 足 够 快 两 次 试 验 间 的 间 隔 足 够 短 , 那 么 就 可 以 用 试 验 次 数 近 似 表 示 试 验 时 间 , 也 就 是 说 极 限 负 二 项 分 布 近 似 于 伽 马 分 布 \textcolor{blue}{伽马分布与负二项分布:}类似于指数分布和几何分布。伽马分布变量是完成n次成功试验时所等待的时间,\\ 而负二项分布变量是完成n次成功试验时所经历的所有试验次数,如果每次试验测试的足够快两次试验间的间隔足够短,\\ 那么就可以用试验次数近似表示试验时间,也就是说极限负二项分布近似于伽马分布 伽马分布与负二项分布:类似于指数分布和几何分布。伽马分布变量是完成n次成功试验时所等待的时间,而负二项分布变量是完成n次成功试验时所经历的所有试验次数,如果每次试验测试的足够快两次试验间的间隔足够短,那么就可以用试验次数近似表示试验时间,也就是说极限负二项分布近似于伽马分布
- 伽 马 分 布 与 泊 松 分 布 的 共 轭 性 : 与 之 前 贝 塔 分 布 是 二 项 分 布 概 率 未 知 的 共 轭 性 先 验 一 样 , 伽 马 分 布 是 一 个 泊 松 分 布 速 率 未 知 的 共 轭 先 验 。 对 于 一 个 未 知 速 率 λ 使 其 先 验 分 布 服 从 伽 马 分 布 且 λ ∼ G a m m a ( r 0 , b 0 ) 对 试 验 观 察 t 时 间 , 并 且 记 录 这 期 间 的 成 功 次 数 Y , 则 Y 的 条 件 分 布 服 从 泊 松 分 布 , 即 Y = y ∣ λ ∼ P o i s ( λ t ) 由 Y 的 条 件 分 布 和 λ 的 边 缘 分 布 构 成 一 个 联 合 分 布 f Y , λ ( y , λ ) = P ( Y = y ∣ λ ) f ( λ ) , 并 求 出 Y 的 边 缘 概 率 质 量 函 数 P ( Y = y ) = ∫ 0 + ∞ P ( Y = y ∣ λ ) f ( λ ) d λ 。 于 是 基 于 上 面 的 条 件 我 们 可 以 利 用 贝 叶 斯 准 则 求 出 λ 的 后 验 概 率 密 度 函 数 和 后 验 分 布 f ( λ ∣ Y = y ) = P ( Y = y ∣ λ ) f ( λ ) P ( Y = y ) ⇒ f ( λ ∣ Y = y ) ∝ e − ( b 0 + t ) λ λ r 0 + y 1 λ ⇒ λ ∣ Y = y ∼ G a m m a ( r 0 + y , b 0 + t ) 即 λ 的 后 验 分 布 依 然 服 从 伽 马 分 布 , 因 此 伽 马 分 布 是 泊 松 分 布 的 共 轭 先 验 解 释 一 下 伽 马 分 布 中 的 参 数 , r 0 为 事 先 观 察 到 的 在 b 0 时 间 内 成 功 的 次 数 , 随 后 在 经 过 t 时 间 后 有 观 察 到 的 成 功 次 数 为 y 则 累 计 的 成 功 次 数 为 r 0 + y , 累 积 的 观 察 时 间 为 b 0 + t 。 由 此 我 们 可 以 知 道 使 用 伽 马 分 布 的 意 义 , 通 过 泊 松 分 布 试 验 收 集 数 据 并 通 过 贝 叶 斯 准 则 不 断 更 新 未 知 速 率 λ 的 先 验 后 验 的 循 环 过 程 是 被 包 含 在 伽 马 分 布 族 中 进 行 的 , 也 就 是 说 我 们 只 需 要 更 新 伽 马 分 布 的 两 个 参 数 , 历 史 累 积 成 功 次 数 和 历 史 累 积 观 察 时 间 , 就 可 以 从 更 新 后 的 伽 马 分 布 中 生 成 一 个 符 合 试 验 结 果 的 成 功 速 率 λ 并 加 以 使 用 \textcolor{blue}{伽马分布与泊松分布的共轭性:}与之前贝塔分布是二项分布概率未知的共轭性先验一样,伽马分布是一个泊松分布速率未知的共轭先验。\\ \displaystyle对于一个未知速率\lambda使其先验分布服从伽马分布且\lambda \sim Gamma(r_0,b_0)\\ \displaystyle对试验观察t时间,并且记录这期间的成功次数Y,则Y的条件分布服从泊松分布,即Y=y|\lambda\sim Pois(\lambda t)\\ \displaystyle由Y的条件分布和\lambda的边缘分布构成一个联合分布f_{Y,\lambda}(y,\lambda)=P(Y=y|\lambda)f(\lambda),\\ \displaystyle并求出Y的边缘概率质量函数P(Y=y)=\int_0^{+\infty}P(Y=y|\lambda)f(\lambda)d\lambda。\\ 于是基于上面的条件我们可以利用贝叶斯准则求出\lambda的后验概率密度函数和后验分布\\ f(\lambda|Y=y)=\frac{P(Y=y|\lambda)f(\lambda)}{P(Y=y)}\rArr f(\lambda|Y=y) \propto e^{-(b_0+t)\lambda}\lambda^{r_0+y}\frac{1}{\lambda}\rArr \lambda|Y=y \sim Gamma(r_0+y,b_0+t)\\ 即\lambda的后验分布依然服从伽马分布,因此伽马分布是泊松分布的共轭先验\\ \textcolor{blue}{解释一下伽马分布中的参数},r_0为事先观察到的在b_0时间内成功的次数,随后在经过t时间后有观察到的成功次数为y\\ 则累计的成功次数为r_0+y,累积的观察时间为b_0+t。\\ 由此我们可以知道\textcolor{blue}{使用伽马分布的意义},通过泊松分布试验收集数据并通过贝叶斯准则不断更新未知速率\lambda的先验后验的循环过程\\ 是被包含在伽马分布族中进行的,也就是说我们只需要更新伽马分布的两个参数,历史累积成功次数和历史累积观察时间,\\ 就可以从更新后的伽马分布中生成一个符合试验结果的成功速率\lambda并加以使用 伽马分布与泊松分布的共轭性:与之前贝塔分布是二项分布概率未知的共轭性先验一样,伽马分布是一个泊松分布速率未知的共轭先验。对于一个未知速率λ使其先验分布服从伽马分布且λ∼Gamma(r0,b0)对试验观察t时间,并且记录这期间的成功次数Y,则Y的条件分布服从泊松分布,即Y=y∣λ∼Pois(λt)由Y的条件分布和λ的边缘分布构成一个联合分布fY,λ(y,λ)=P(Y=y∣λ)f(λ),并求出Y的边缘概率质量函数P(Y=y)=∫0+∞P(Y=y∣λ)f(λ)dλ。于是基于上面的条件我们可以利用贝叶斯准则求出λ的后验概率密度函数和后验分布f(λ∣Y=y)=P(Y=y)P(Y=y∣λ)f(λ)⇒f(λ∣Y=y)∝e−(b0+t)λλr0+yλ1⇒λ∣Y=y∼Gamma(r0+y,b0+t)即λ的后验分布依然服从伽马分布,因此伽马分布是泊松分布的共轭先验解释一下伽马分布中的参数,r0为事先观察到的在b0时间内成功的次数,随后在经过t时间后有观察到的成功次数为y则累计的成功次数为r0+y,累积的观察时间为b0+t。由此我们可以知道使用伽马分布的意义,通过泊松分布试验收集数据并通过贝叶斯准则不断更新未知速率λ的先验后验的循环过程是被包含在伽马分布族中进行的,也就是说我们只需要更新伽马分布的两个参数,历史累积成功次数和历史累积观察时间,就可以从更新后的伽马分布中生成一个符合试验结果的成功速率λ并加以使用
- 伽 马 分 布 与 泊 松 分 布 的 共 轭 性 的 补 充 说 明 : 可 以 看 出 在 泊 松 分 布 的 共 轭 先 验 分 布 中 , 伽 马 分 布 的 第 二 个 参 数 已 不 再 是 成 功 速 率 λ 而 是 历 史 累 积 观 察 时 间 , 这 是 需 要 注 意 的 另 外 E ( λ ∣ Y = y ) = r 0 + y b 0 + t , V a r ( λ ∣ Y = y ) = r 0 + y ( b 0 + t ) 2 \textcolor{blue}{伽马分布与泊松分布的共轭性的补充说明:} 可以看出在泊松分布的共轭先验分布中,伽马分布的第二个参数已不再是成功速率\lambda而是历史累积观察时间,这是需要注意的\\ \displaystyle 另外E(\lambda|Y=y)=\frac{r_0+y}{b_0+t},\quad Var(\lambda|Y=y)=\frac{r_0+y}{(b_0+t)^2} 伽马分布与泊松分布的共轭性的补充说明:可以看出在泊松分布的共轭先验分布中,伽马分布的第二个参数已不再是成功速率λ而是历史累积观察时间,这是需要注意的另外E(λ∣Y=y)=b0+tr0+y,Var(λ∣Y=y)=(b0+t)2r0+y
5. 贝塔分布与伽马分布的关系
令 X ∼ G a m m a ( a , λ ) , Y ∼ G a m m a ( b , λ ) , X 与 Y 相 互 独 立 , T = X + Y , W = X X + Y 则 T ∼ G a m m a ( a + b , λ ) , W ∼ B e t a ( a , b ) , 且 T 与 W 相 互 独 立 , E ( W ) = a a + b \displaystyle 令X\sim Gamma(a,\lambda),Y\sim Gamma(b,\lambda),X与Y相互独立,T=X+Y,W= \frac{X}{X+Y}\\ 则T \sim Gamma(a+b,\lambda),W\sim Beta(a,b),且T与W相互独立,E(W)=\frac{a}{a+b} 令X∼Gamma(a,λ),Y∼Gamma(b,λ),X与Y相互独立,T=X+Y,W=X+YX则T∼Gamma(a+b,λ),W∼Beta(a,b),且T与W相互独立,E(W)=a+ba
6. 顺序统计量
将 n 个 随 机 变 量 X 1 , ⋯ , X n 由 小 到 大 进 行 排 序 , 得 到 变 换 后 的 随 机 变 量 序 列 X ( 1 ) , ⋯ , X ( n ) , X ( j ) 表 示 变 量 在 序 列 中 排 第 j 位 。 则 变 换 后 的 随 机 变 量 序 列 被 称 为 顺 序 统 计 量 将n个随机变量X_1,\cdots,X_n由小到大进行排序,得到变换后的随机变量序列X_{(1)},\cdots,X_{(n)},X_{(j)}表示变量在序列中排第j位。\\则变换后的随机变量序列被称为顺序统计量 将n个随机变量X1,⋯,Xn由小到大进行排序,得到变换后的随机变量序列X(1),⋯,X(n),X(j)表示变量在序列中排第j位。则变换后的随机变量序列被称为顺序统计量
- 顺 序 统 计 量 的 累 积 分 布 函 数 : 令 X 1 , ⋯ , X n 是 独 立 同 步 分 布 的 连 续 性 随 机 变 量 序 列 , 其 累 积 分 布 函 数 为 F , 则 第 j 个 顺 序 统 计 量 X ( j ) 的 累 积 分 布 函 数 为 P ( X ( j ) ≤ x ) = ∑ k = j n ( n k ) F ( x ) k ( 1 − F ( x ) ) n − k \textcolor{blue}{顺序统计量的累积分布函数:}令X_1,\cdots,X_n是独立同步分布的连续性随机变量序列,其累积分布函数为F,\\则第j个\displaystyle顺序统计量X_{(j)}的累积分布函数为P(X_{(j)}\leq x)=\sum_{k=j}^n{n \choose k}F(x)^k(1-F(x))^{n-k} 顺序统计量的累积分布函数:令X1,⋯,Xn是独立同步分布的连续性随机变量序列,其累积分布函数为F,则第j个顺序统计量X(j)的累积分布函数为P(X(j)≤x)=k=j∑n(kn)F(x)k(1−F(x))n−k
- 顺 序 统 计 量 的 概 率 密 度 函 数 : 令 X 1 , ⋯ , X n 是 独 立 同 步 分 布 的 连 续 性 随 机 变 量 序 列 , 其 累 积 分 布 函 数 为 F , 概 率 密 度 函 数 为 f 则 第 j 个 顺 序 统 计 量 X ( j ) 的 概 率 密 度 函 数 为 f X ( j ) ( x ) = n ( n − 1 j − 1 ) f ( x ) F ( x ) j − 1 ( 1 − F ( x ) ) n − j 从 式 子 结 构 可 以 看 出 , 一 般 的 顺 序 统 计 量 不 会 服 从 一 个 一 致 的 分 布 , 但 服 从 标 准 均 匀 分 布 的 顺 序 统 计 量 的 联 合 分 布 是 个 例 外 \textcolor{blue}{顺序统计量的概率密度函数:}令X_1,\cdots,X_n是独立同步分布的连续性随机变量序列,其累积分布函数为F,概率密度函数为f\\ \displaystyle则第j个顺序统计量X_{(j)}的概率密度函数为f_{X_{(j)}}(x)=n{n-1 \choose j-1}f(x)F(x)^{j-1}(1-F(x))^{n-j}\\ 从式子结构可以看出,一般的顺序统计量不会服从一个一致的分布,但服从标准均匀分布的顺序统计量的联合分布是个例外 顺序统计量的概率密度函数:令X1,⋯,Xn是独立同步分布的连续性随机变量序列,其累积分布函数为F,概率密度函数为f则第j个顺序统计量X(j)的概率密度函数为fX(j)(x)=n(j−1n−1)f(x)F(x)j−1(1−F(x))n−j从式子结构可以看出,一般的顺序统计量不会服从一个一致的分布,但服从标准均匀分布的顺序统计量的联合分布是个例外
- 均 匀 分 布 中 的 次 序 统 计 量 : 令 U 1 , ⋯ , U n 是 独 立 同 步 分 布 于 U ( 0 , 1 ) 的 随 机 变 量 。 那 么 当 0 ≤ x ≤ 1 时 , f ( x ) = 1 , 且 F ( x ) = x , 因 此 U ( j ) 的 概 率 密 度 函 数 为 f U ( j ) ( x ) = n ( n − 1 j − 1 ) x j − 1 ( 1 − F ( x ) ) n − j 这 正 是 B e t a ( j , n − j + 1 ) 的 概 率 密 度 函 数 。 因 此 U ( j ) ∼ B e t a ( j , n − j + 1 ) , 且 E ( U ( j ) ) = j n + 1 \textcolor{blue}{均匀分布中的次序统计量:}令U_1,\cdots,U_n是独立同步分布于U(0,1)的随机变量。那么当0\leq x\leq 1时,\\ \displaystyle f(x)=1,且F(x)=x,因此U_{(j)}的概率密度函数为f_{U_{(j)}}(x)=n{n-1 \choose j-1}x^{j-1}(1-F(x))^{n-j}\\ \displaystyle 这正是Beta(j,n-j+1)的概率密度函数。因此U(j)\sim Beta(j,n-j+1),且E(U_{(j)})=\frac{j}{n+1} 均匀分布中的次序统计量:令U1,⋯,Un是独立同步分布于U(0,1)的随机变量。那么当0≤x≤1时,f(x)=1,且F(x)=x,因此U(j)的概率密度函数为fU(j)(x)=n(j−1n−1)xj−1(1−F(x))n−j这正是Beta(j,n−j+1)的概率密度函数。因此U(j)∼Beta(j,n−j+1),且E(U(j))=n+1j
二. 条件期望
1. 两种条件期望
1). 基于给定事件的条件期望
设A是发生概率大于0的事件,则在给定A的条件下
{
E
(
Y
∣
A
)
=
∑
y
y
P
(
Y
=
y
∣
A
)
Y
是
离
散
型
变
量
E
(
Y
∣
A
)
=
∫
−
∞
+
∞
y
f
(
Y
=
y
∣
A
)
d
y
Y
是
连
续
型
变
量
f
(
y
∣
A
)
=
P
(
A
∣
Y
=
y
)
f
(
x
)
P
(
A
)
Y
是
连
续
型
变
量
\begin{cases} \displaystyle E(Y|A)&=\displaystyle\sum_yyP(Y=y|A)&Y是离散型变量\\ \displaystyle E(Y|A)&=\displaystyle\int_{-\infty}^{+\infty}yf(Y=y|A)dy&Y是连续型变量\\ \displaystyle f(y|A)&=\displaystyle\frac{P(A|Y=y)f(x)}{P(A)}&Y是连续型变量\\ \end{cases}
⎩⎪⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎪⎪⎧E(Y∣A)E(Y∣A)f(y∣A)=y∑yP(Y=y∣A)=∫−∞+∞yf(Y=y∣A)dy=P(A)P(A∣Y=y)f(x)Y是离散型变量Y是连续型变量Y是连续型变量
简单来说就是
A
A
A的状态如何,是否发生已经给定了,即
A
A
A变量的值已知,所以
E
(
Y
∣
A
)
E(Y|A)
E(Y∣A)是一个常数
样本法
E
(
Y
)
≈
1
n
∑
j
=
1
n
y
j
E
(
Y
∣
A
)
≈
∑
j
=
1
n
y
j
I
j
∑
j
=
1
n
I
j
做
n
次
试
验
,
y
j
表
示
每
一
次
事
件
Y
的
状
态
,
I
j
为
表
示
事
件
A
是
否
发
生
的
示
性
随
机
变
量
,
1
为
发
生
,
0
为
没
发
生
\displaystyle E(Y)\approx\frac{1}{n}\sum_{j=1}^ny_j\qquad E(Y|A)\approx\frac{\displaystyle\sum_{j=1}^ny_jI_j}{\displaystyle\sum_{j=1}^nI_j}\\ 做n次试验,y_j表示每一次事件Y的状态,I_j为表示事件A是否发生的示性随机变量,1为发生,0为没发生
E(Y)≈n1j=1∑nyjE(Y∣A)≈j=1∑nIjj=1∑nyjIj做n次试验,yj表示每一次事件Y的状态,Ij为表示事件A是否发生的示性随机变量,1为发生,0为没发生
2). 基于给定随机变量的条件期望
- 对 于 随 机 变 量 X 而 言 , E ( Y ∣ X ) 实 质 为 X 的 函 数 g ( X ) 对于随机变量X而言,E(Y|X)实质为X的函数g(X) 对于随机变量X而言,E(Y∣X)实质为X的函数g(X)
- 在 已 知 X 的 情 况 下 , 某 种 意 义 上 E ( Y ∣ X ) 是 Y 的 最 佳 预 测 在已知X的情况下,某种意义上E(Y|X)是Y的最佳预测 在已知X的情况下,某种意义上E(Y∣X)是Y的最佳预测
- 对 于 2 的 应 用 最 多 的 就 是 线 性 回 归 和 曲 线 拟 合 , 拟 合 出 来 的 直 线 或 曲 线 其 实 就 是 一 个 关 于 X 的 函 数 , 是 Y 关 于 X 的 条 件 期 望 对于2的应用最多的就是线性回归和曲线拟合,拟合出来的直线或曲线其实就是一个关于X的函数,是Y关于X的条件期望 对于2的应用最多的就是线性回归和曲线拟合,拟合出来的直线或曲线其实就是一个关于X的函数,是Y关于X的条件期望
- { E ( Y ∣ X = x ) = ∑ y y P ( Y = y ∣ X = x ) Y 是 离 散 型 变 量 E ( Y ∣ X = x ) = ∫ − ∞ + ∞ y f Y ∣ X ( y ∣ x ) d y Y 是 连 续 型 变 量 \begin{cases} E(Y|X=x)&\displaystyle =\sum_yyP(Y=y|X=x)&Y是离散型变量\\ E(Y|X=x)&\displaystyle =\int_{-\infty}^{+\infty}yf_{Y|X}(y|x)dy&Y是连续型变量\\ \end{cases} ⎩⎪⎪⎨⎪⎪⎧E(Y∣X=x)E(Y∣X=x)=y∑yP(Y=y∣X=x)=∫−∞+∞yfY∣X(y∣x)dyY是离散型变量Y是连续型变量
2. 条件期望的性质
- 如 果 X 与 Y 相 互 独 立 则 E ( Y ∣ X ) = E ( Y ) 如果X与Y相互独立则E(Y|X)=E(Y) 如果X与Y相互独立则E(Y∣X)=E(Y)
- 对 于 X 的 任 意 函 数 h , E ( h ( X ) Y ∣ X ) = h ( X ) E ( Y ∣ X ) 对于X的任意函数h,E(h(X)\,Y|X)=h(X)\,E(Y|X) 对于X的任意函数h,E(h(X)Y∣X)=h(X)E(Y∣X)
- E ( Y 1 + Y 2 ∣ X ) = E ( Y 1 ∣ X ) + E ( Y 2 ∣ X ) , 对 于 常 数 c , E ( c Y ∣ X ) = c E ( Y ∣ X ) E(Y_1+Y_2|X)=E(Y_1|X)+E(Y_2|X),对于常数c,E(cY|X)=cE(Y|X) E(Y1+Y2∣X)=E(Y1∣X)+E(Y2∣X),对于常数c,E(cY∣X)=cE(Y∣X)
- 亚 当 定 律 : E ( E ( Y ∣ X ) ) = E ( Y ) , 这 是 连 接 条 件 期 望 与 非 条 件 期 望 的 关 键 , 另 外 E ( E ( Y ∣ Z , X ) ∣ Z ) = E ( Y ∣ X ) \color{blue}亚当定律 :E(E(Y|X))=E(Y),这是连接条件期望与非条件期望的关键,另外E(E(Y|Z,X)|Z)=E(Y|X) 亚当定律:E(E(Y∣X))=E(Y),这是连接条件期望与非条件期望的关键,另外E(E(Y∣Z,X)∣Z)=E(Y∣X)
- 投 影 解 释 : 随 机 变 量 Y − E ( Y ∣ X ) 是 用 X 预 测 Y 的 残 差 。 他 与 任 意 函 数 h ( X ) 是 不 相 关 的 。 用 一 个 三 维 空 间 的 直 线 投 影 到 一 个 平 面 上 的 例 子 去 类 比 就 好 。 也 就 是 说 E ( Y ∣ X ) 是 在 h ( X ) 平 面 内 最 接 近 Y 的 函 数 了 投影解释:随机变量Y-E(Y|X)是用X预测Y的残差。他与任意函数h(X)是不相关的。\\ 用一个三维空间的直线投影到一个平面上的例子去类比就好。也就是说E(Y|X)是在h(X)平面内最接近Y的函数了 投影解释:随机变量Y−E(Y∣X)是用X预测Y的残差。他与任意函数h(X)是不相关的。用一个三维空间的直线投影到一个平面上的例子去类比就好。也就是说E(Y∣X)是在h(X)平面内最接近Y的函数了
3. 条件方差
在给定
X
X
X的条件下,
Y
Y
Y的条件方差为
V
a
r
(
Y
∣
X
)
=
E
(
(
Y
−
E
(
Y
∣
X
)
)
2
∣
X
)
⟹
V
a
r
(
Y
∣
X
)
=
E
(
Y
2
∣
X
)
−
(
E
(
Y
∣
X
)
)
2
Var(Y|X)=E((Y-E(Y|X))^2|X)\implies Var(Y|X)=E(Y^2|X)-(E(Y|X))^2
Var(Y∣X)=E((Y−E(Y∣X))2∣X)⟹Var(Y∣X)=E(Y2∣X)−(E(Y∣X))2
夏娃定律
V
a
r
(
Y
)
=
E
(
V
a
r
(
Y
∣
X
)
)
+
V
a
r
(
E
(
Y
∣
X
)
)
,
很
明
显
夏
娃
定
律
是
将
非
条
件
方
差
与
条
件
方
差
和
条
件
期
望
结
合
起
来
的
定
律
其
中
我
们
称
E
(
V
a
r
(
Y
∣
X
)
)
为
组
内
变
动
,
V
a
r
(
E
(
Y
∣
X
)
)
为
组
间
变
动
,
也
就
是
说
把
一
个
X
的
值
看
做
分
组
的
标
准
,
V
a
r
(
Y
∣
X
=
x
j
)
很
显
然
是
一
个
组
内
的
Y
值
波
动
的
表
示
,
其
后
求
均
值
E
(
V
a
r
(
Y
∣
X
)
)
只
是
这
个
波
动
的
均
数
而
V
a
r
(
E
(
Y
∣
X
)
)
先
求
的
是
组
内
的
期
望
均
值
,
然
后
看
这
个
值
在
各
个
组
间
的
波
动
,
即
求
V
a
r
(
E
(
Y
∣
X
)
)
Var(Y)=E(Var(Y|X))+Var(E(Y|X)),很明显夏娃定律是将非条件方差与条件方差和条件期望结合起来\\的定律 其中我们称E(Var(Y|X))为\textcolor{blue}{组内变动},Var(E(Y|X))为\textcolor{blue}{组间变动},也就是说把一个X的值看做分组的标准,\\ Var(Y|X=x_j)很显然是一个组内的Y值波动的表示,其后求均值E(Var(Y|X))只是这个波动的均数\\ 而Var(E(Y|X))先求的是组内的期望均值,然后看这个值在各个组间的波动,即求Var(E(Y|X))
Var(Y)=E(Var(Y∣X))+Var(E(Y∣X)),很明显夏娃定律是将非条件方差与条件方差和条件期望结合起来的定律其中我们称E(Var(Y∣X))为组内变动,Var(E(Y∣X))为组间变动,也就是说把一个X的值看做分组的标准,Var(Y∣X=xj)很显然是一个组内的Y值波动的表示,其后求均值E(Var(Y∣X))只是这个波动的均数而Var(E(Y∣X))先求的是组内的期望均值,然后看这个值在各个组间的波动,即求Var(E(Y∣X))
4. 线性回归
线性回归模型使用单个解释变量
X
X
X来预测响应变量
Y
Y
Y,并且假设
Y
Y
Y的条件期望与
X
X
X呈线性关系即
E
(
Y
∣
X
)
=
a
+
b
X
E(Y|X)=a+bX
E(Y∣X)=a+bX,或者说
E
(
Y
∣
X
)
=
a
+
b
X
+
ε
,
其
中
ε
是
一
个
随
机
变
量
称
为
误
差
,
E
(
ε
∣
X
)
=
0
E(Y|X)=a+bX+\varepsilon,其中\varepsilon是一个随机变量称为误差,E(\varepsilon|X)=0
E(Y∣X)=a+bX+ε,其中ε是一个随机变量称为误差,E(ε∣X)=0,显然我们想要求出
E
(
Y
∣
X
)
E(Y|X)
E(Y∣X)
就需要求出
a
a
a斜率和
b
b
b截距,
b
=
C
o
v
(
X
,
Y
)
V
a
r
(
X
)
,
a
=
E
(
Y
)
−
b
E
(
X
)
=
E
(
Y
)
−
C
o
v
(
X
,
Y
)
V
a
r
(
X
)
⋅
E
(
X
)
\displaystyle b=\frac{Cov(X,Y)}{Var(X)},a=E(Y)-bE(X)=E(Y)-\frac{Cov(X,Y)}{Var(X)}\cdot E(X)
b=Var(X)Cov(X,Y),a=E(Y)−bE(X)=E(Y)−Var(X)Cov(X,Y)⋅E(X)