复旦大学经济学院 ECON130001
先研究联合、边缘、条件分布及独立性,再求由 $X,Y$ 构成的新变量的分布,例如和、积、最大值与最小值。
求和的分布时,将 $X+Y\le z$ 写成平面区域,再用联合分布计算该区域的概率;独立的连续型变量可用卷积求和的密度。
| 分布 | 两个变量各自的分布 | 和的分布 |
|---|---|---|
| 二项 | $X\sim B(n_1,p)$,$Y\sim B(n_2,p)$ | $X+Y\sim B(n_1+n_2,p)$ |
| 泊松 | $X\sim\mathrm{Poisson}(\lambda_1)$,$Y\sim\mathrm{Poisson}(\lambda_2)$ | $X+Y\sim\mathrm{Poisson}(\lambda_1+\lambda_2)$ |
| 正态 | $X\sim N(\mu_1,\sigma_1^2)$,$Y\sim N(\mu_2,\sigma_2^2)$ | $X+Y\sim N(\mu_1+\mu_2,\sigma_1^2+\sigma_2^2)$ |
二项分布还要求成功概率相同;此时合并的是两组具有相同成功概率的独立试验。
知道分布后,怎样用几个数描述平均水平、波动大小,以及两个变量的共同变化?
设航班甲每次延误10分钟;航班乙有一半概率准点,另一半概率延误20分钟。
甲的延误时间固定,乙的延误时间会变化。怎样分别描述平均水平和波动大小?
投资者 A、B 的收益分布不同,谁更好?AI 模型 A、B 的误差分布不同,哪个更优?答案还取决于评价标准。本讲用几个数概括分布的部分特征,并说明它们的适用边界。
某航班过去四天分别延误 $0,10,20,30$ 分钟,平均延误多久?——$15$ 分钟。
换个问法:过去四天里两天延误 10 分钟、一天 20 分钟、一天没延误呢? $$\frac{10+10+20+0}{4}=10\times\tfrac12+20\times\tfrac14+0\times\tfrac14=10$$ 右式是以观察频率为权重的平均数。用分布中的概率代替这些频率,就得到期望的计算形式。
离散型:若级数 $\sum_k x_kf(x_k)$ 绝对收敛,则 $$EX=\sum_{k=1}^{\infty}x_kf(x_k)$$
连续型:若积分 $\int_{-\infty}^{\infty}xf(x)\mathrm{d}x$ 绝对收敛,则 $$EX=\int_{-\infty}^{\infty}xf(x)\,\mathrm{d}x$$
把定义中的取值 $x$ 换成 $|x|$,再按概率求和或积分,结果仍须是有限数。
只有有限个可能取值时,这一条件自动满足;有无限多个取值或取值无界时,需要检查。本讲讨论满足这一条件的有限期望。
$X$ 服从 0-1 分布,$P(X=1)=p$。
$$E(X)=0\times(1-p)+1\times p=p$$伯努利分布的期望就是 $p$ 本身——「平均成功次数」等于「成功概率」。它也可用于计算多个试验中成功总次数的期望。
$X$ 有密度 $f(x)=2x$($0<x<1$)。
$$E(X)=\int_0^1 x\cdot 2x\,\mathrm{d}x=\frac{2}{3}$$直觉核对:密度随 $x$ 增大而增大,取值更偏向区间右端,因此期望应大于中点 $0.5$;$\frac23$ 与此一致。
盒中红球比例 $p$,有放回取 $n$ 个,其中 $X$ 个红球。求 $EX$。
关键一步:$k\,C_n^k=k\cdot\dfrac{n!}{k!(n-k)!}=n\cdot\dfrac{(n-1)!}{(k-1)!(n-k)!}=n\,C_{n-1}^{k-1}$,于是
$$EX=np\sum_{k=1}^{n}C_{n-1}^{k-1}p^{k-1}(1-p)^{(n-1)-(k-1)}=np\cdot 1=np$$末尾那个和式是 $B(n-1,p)$ 的全部概率,等于 1。
第15页将利用期望的可加性给出另一种计算方法。
$$EX=\sum_{k=1}^{\infty}k(1-p)^{k-1}p$$ 两边乘 $(1-p)$ 并把指标平移: $$(1-p)EX=\sum_{k=1}^{\infty}k(1-p)^{k}p=\sum_{k=1}^{\infty}(k-1)(1-p)^{k-1}p$$ 上下两式相减,$k-(k-1)=1$: $$pEX=\sum_{k=1}^{\infty}(1-p)^{k-1}p=1 \ \Longrightarrow\ EX=\frac1p$$
结果的含义:成功率 $p$,平均要试 $1/p$ 次。$p=0.1$ 就平均试 10 次。
以上级数在 $0<p<1$ 时绝对收敛,因此可以逐项相减;$p=1$ 时 $X=1$,结论同样成立。
$x=0$ 那项为 0,从 1 开始;提出一个 $\lambda$ 后剩下的又是全部概率。
分部积分。速率 $\lambda$ 的倒数就是平均等待时间。例如,泊松到达模型中速率为每小时4次,平均等待时间为15分钟。
正态分布有有限期望,由密度关于 $\mu$ 对称可得。
| 分布 | $EX$ | 怎么记 |
|---|---|---|
| 伯努利 $(p)$ | $p$ | 平均成功次数 = 成功概率 |
| 二项 $B(n,p)$ | $np$ | 做 $n$ 次,每次期望 $p$ |
| 几何 $(p)$ | $1/p$ | 成功率的倒数 |
| 泊松 $(\lambda)$ | $\lambda$ | 由泊松概率函数计算得到 |
| 指数 $(\lambda)$ | $1/\lambda$ | 速率的倒数 |
| 正态 $N(\mu,\sigma^2)$ | $\mu$ | 对称中心 |
表中结果可直接用于后续的期望计算;使用前先确认分布及参数含义。
令 $Y=r(X)$,若 $E|r(X)|<\infty$,则 $$E[r(X)]=\sum_{\text{所有 }x}r(x)f(x)\qquad\text{或}\qquad E[r(X)]=\int_{-\infty}^{\infty}r(x)f(x)\,\mathrm{d}x$$
第5讲介绍了求 $Y=r(X)$ 分布的方法。 但如果只要期望,用 $X$ 的密度直接积就行,不必求 $Y$ 的分布。
$f(x)=2x$($0<x<1$),求 $E(X^{1/2})$。
$$\int_0^1 x^{1/2}\cdot2x\,\mathrm{d}x=\frac45$$$f(x)=3x^2$($0<x<1$),求 $E(1/X)$。
$$\int_0^1 \frac1x\cdot3x^2\,\mathrm{d}x=\frac32$$$f(x,y)=1$($0\le x,y\le1$),求 $E(X^2+Y^2)$。
$$\int_0^1\!\!\int_0^1(x^2+y^2)\,\mathrm{d}x\,\mathrm{d}y=\frac13+\frac13=\frac23$$$f(x,y)=12y^2$($0\le y\le x\le1$),求 $E(XY)$。
$$\int_0^1\!\!\int_0^x xy\cdot12y^2\,\mathrm{d}y\,\mathrm{d}x =\int_0^1 3x^5\,\mathrm{d}x=\frac12$$积分限由三角形取值区域确定:$y$ 从 0 到 $x$。
若 $E|r(X,Y)|<\infty$,则 $E[r(X,Y)]=\displaystyle\iint r(x,y)f(x,y)\,\mathrm{d}x\,\mathrm{d}y$——同样不必先求 $r(X,Y)$ 的分布。
以下假设 $E|X|$、$E|Y|$ 均有限。
只要期望存在,$E(X+Y)=E(X)+E(Y)$ 就成立,不要求独立。
独立可以推出 $E(XY)=E(X)E(Y)$,反过来不一定;后面会给出反例。
$X\sim B(n,p)$ 可写成 $n$ 个独立伯努利之和 $X=\sum_{i=1}^n X_i$, 每个 $E(X_i)=p$。由可加性 $$EX=\sum_{i=1}^{n}E(X_i)=np$$
直接利用可加性,可省去第8页的组合恒等式推导。
$n$ 封信随机装进 $n$ 个信封,$X$ 为装对的封数。求 $EX$。
令 $X_i=1$ 表示第 $i$ 封装对,则 $P(X_i=1)=\frac1n$,$E(X_i)=\frac1n$: $$EX=\sum_{i=1}^{n}E(X_i)=n\cdot\frac1n=1$$
无论 $n$ 多大,平均装对的封数都是 1。 当 $n\ge2$ 时,这些 $X_i$ 并不独立 ——可加性不需要独立,这里就用上了。
$X_1,X_2,X_3$ 独立,$E(X_i)=0$、$E(X_i^2)=1$。求 $E[X_1^2(X_2-4X_3)^2]$。
$X_1^2$ 与 $(X_2-4X_3)^2$ 独立,故可乘: $$=E(X_1^2)\cdot E[(X_2-4X_3)^2]$$ $$=1\cdot E(X_2^2-8X_2X_3+16X_3^2)=1-0+16=17$$ 中间项 $E(X_2X_3)=E(X_2)E(X_3)=0$。
第一次滤去 $X_1$ 比例杂质,第二次滤去剩余杂质中的 $X_2$ 比例。 $X_1,X_2$ 独立同分布,$f(x)=4x^3$($0<x<1$)。 求两次后剩余比例 $Y$ 的期望。
$Y=(1-X_1)(1-X_2)$,由独立性 $$E(Y)=\big[1-E(X_1)\big]^2=\Big(1-\tfrac45\Big)^2=0.04$$ ($E(X_1)=\int_0^1 x\cdot4x^3=\frac45$。)
两级过滤后的剩余杂质比例的期望为4%;单级过滤后的相应期望为20%。
组合总收益率 $R_p=\sum_i w_iR_i$($w_i$ 是各资产权重)。
组合的预期收益就是各资产预期收益的加权平均——不需要知道资产之间怎么联动。
这个公式回答了组合平均能获得多少收益。即使平均收益相同,实际收益的波动也可能不同;下一部分再学习怎样描述这种波动。
给定 $X=x$,且相应条件分布存在,且在该条件下 $|Y|$ 的期望有限时,$Y$ 的条件期望,就是条件分布的期望: $$E[Y\mid x]=\sum_{\text{所有 }y}y\,g_2(y\mid x) \qquad\text{或}\qquad E[Y\mid x]=\int_{-\infty}^{\infty}y\,g_2(y\mid x)\,\mathrm{d}y$$
固定 $x$ 后,$E[Y\mid x]$ 是一个数;随 $x$ 变化,这些数构成一个函数。 第6讲比较的是给定评分后的高分概率;这里进一步考察条件平均值:给定用户对甲片打 5 分,他对乙片的平均打分。
甲的延误始终为10分钟;乙以各1/2的概率延误0分钟或20分钟。两者的期望相同,但甲没有波动,乙会偏离平均值。
能否用一个非负的数,概括随机变量偏离自身期望的程度?
$E(X^2)$ 是先将 $X$ 平方,再求期望,称为二阶矩。它通常不同于 $(EX)^2$。以下假设 $E(X^2)<\infty$。
称 $X-EX$ 为离差。离差平方的期望称为方差: $$D(X)=E\big[(X-EX)^2\big]$$ 离散型 $\sum_x(x-EX)^2f(x)$,连续型 $\int(x-EX)^2f(x)\mathrm{d}x$。 其平方根 $\sqrt{D(X)}$ 称为标准差。
离差的期望为0,正负偏差会抵消;平方后不再抵消,也便于推导和的方差。绝对偏差也是可用的离散程度指标。
方差的单位是原单位的平方;标准差与原变量单位相同。
股价 $A\sim U[25,35]$,$B\sim U[15,45]$。求各自的期望与方差。
$E(A)=E(B)=30$——期望完全相同。
$$D(A)=\int_{25}^{35}(a-30)^2\cdot\tfrac1{10}\,\mathrm{d}a=\frac{25}{3}\approx8.33$$ $$D(B)=\int_{15}^{45}(b-30)^2\cdot\tfrac1{30}\,\mathrm{d}b=75$$标准差分别约为 $2.89$ 与 $8.66$,后者约为前者的3倍。
一般地,$U[a,b]$ 的方差是 $\dfrac{(b-a)^2}{12}$——区间宽度变为原来的3倍,方差变为原来的9倍。本例均值相同,方差与标准差反映了取值范围不同所带来的波动差异。
以下假设 $X,Y$ 均有有限二阶矩。
$E(X+Y)=E(X)+E(Y)$ 不需要独立; $D(X+Y)=D(X)+D(Y)$ 在独立时成立;第28页将给出一般公式:协方差为0时,方差也可直接相加 。
而且 $D$ 里的系数是 $a^2$ 不是 $a$——所以标准差才是 $|a|$ 倍。
$$(1-p)E(X^2)=\sum_{k\ge1}(k-1)^2(1-p)^{k-1}p$$ 相减,$k^2-(k-1)^2=2k-1$: $$pE(X^2)=E(2X-1)=\frac2p-1$$ $$D(X)=E(X^2)-(EX)^2=\frac{2-p}{p^2}-\frac1{p^2}=\frac{1-p}{p^2}$$
先计算 $E[X(X-1)]$: $$E(X^2)=E[X(X-1)]+EX=n(n-1)p^2+np$$ $$D(X)=n^2p^2-np^2+np-n^2p^2=np(1-p)$$
利用独立性:$X=\sum X_i$ 独立, 每个 $D(X_i)=p(1-p)$,由可加性 $$D(X)=np(1-p)$$
注意右边这次用到了独立——而 第15页 求期望时没用到。 同一个分解,方差比期望多一个前提。
期望与方差都等于 $\lambda$,这是该分布的一项性质,但不能仅凭期望与方差相等判定分布类型。
标准差为 $1/\lambda$,与期望相等;但仅凭这一等式不能判定变量服从指数分布。
记号里那个 $\sigma^2$ 本来就是方差;标准差为 $\sigma$。
设物体的真实长度为常数 $c$,测量结果为随机变量 $Y$。误差为 $Y-c$,平方误差为 $(Y-c)^2$。平方误差的期望称为均方误差。
若 $E(Y^2)<\infty$,则
$$E[(Y-c)^2]=D(Y)+(EY-c)^2$$方差描述测量值围绕自身均值的波动;后一项描述平均测量值与真实值的偏离。
真实长度 $c=100$ 毫米,$EY=100.1$ 毫米,$D(Y)=0.04$ 平方毫米,则
$$E[(Y-c)^2]=0.04+(100.1-100)^2=0.05\ \text{平方毫米}$$| 分布 | $EX$ | $D(X)$ |
|---|---|---|
| 伯努利 $(p)$ | $p$ | $p(1-p)$ |
| 二项 $B(n,p)$ | $np$ | $np(1-p)$ |
| 几何 $(p)$ | $1/p$ | $(1-p)/p^2$ |
| 泊松 $(\lambda)$ | $\lambda$ | $\lambda$ |
| 均匀 $U[a,b]$ | $(a+b)/2$ | $(b-a)^2/12$ |
| 指数 $(\lambda)$ | $1/\lambda$ | $1/\lambda^2$ |
| 正态 $N(\mu,\sigma^2)$ | $\mu$ | $\sigma^2$ |
比较各分布的参数与数字特征,计算时确认变量的取值范围及参数约定。
怎样描述两个变量偏离各自平均值时的共同变化? 例如,数学成绩较高时,语文成绩是否也倾向于较高?
设 $E(X^2)$ 与 $E(Y^2)$ 均有限;这一条件也称为具有有限二阶矩,见第20页。
$$\mathrm{Cov}(X,Y)=E\big[(X-EX)(Y-EY)\big]$$展开后得到常用的计算式: $$\mathrm{Cov}(X,Y)=E(XY)-E(X)E(Y)$$
协方差等于 $E(XY)$ 与 $E(X)E(Y)$ 的差。独立时差为0;不独立时,差也可能为0。
同号离差产生正的乘积,异号离差产生负的乘积;协方差的符号取决于这些乘积按概率加权后的结果。
第4条给出一般答案:除了两个方差之和,还要计入两倍的协方差。当协方差为0时,两个方差仍可直接相加。
$f(x,y)=2xy+0.5$($0\le x,y\le1$),其他位置为0。求 $\mathrm{Cov}(X,Y)$。
协方差为正。其大小受计量单位影响,后面用相关系数比较线性关联的强弱。
设 $R_p=w_AR_A+w_BR_B$,权重固定且 $w_A,w_B\ge0$、$w_A+w_B=1$。
$$D(R_p)=w_A^2D(R_A)+w_B^2D(R_B)+2w_Aw_B\mathrm{Cov}(R_A,R_B)$$资产的预期收益、方差与权重给定后,协方差降低不会增加组合方差;两项权重均为正时,组合方差严格降低。组合预期收益保持不变。
负协方差使交叉项非正,组合方差不超过两个加权方差项之和。这既不保证组合方差小于任一资产的方差,也不保证组合实际收益优于任一单资产。
两个资产的期望收益与波动率可调,关键是相关系数滑块。
先将相关系数理解为两资产线性关联的标准化指标,定义见下一页;此处权重在0至1之间。
交互演示:拖动相关系数,观察风险—收益曲线如何弯折
$\mathrm{Cov}(2X,Y)=2\mathrm{Cov}(X,Y)$——换个单位,数就变了。 用元还是用万元,算出来差一万倍,无法横向比较。
要求 $0<D(X),D(Y)<\infty$。
$$\rho(X,Y)=\frac{\mathrm{Cov}(X,Y)}{\sqrt{D(X)}\sqrt{D(Y)}}$$可以证明 $|\rho|\le1$。特别地:
以万元记录 $X$,改为以元记录 $X^*=10000X$;保持 $Y$ 的单位不变。
$$\mathrm{Cov}(X^*,Y)=10000\mathrm{Cov}(X,Y)$$正的单位换算系数不会改变相关系数。协方差的数值变大,并不意味着两个变量的关联变强。
$X$ 等可能取 $-1,0,1$,令 $Y=X^2$。
不独立,例如 $P(Y=0\mid X=0)=1\ne P(Y=0)=1/3$——知道 $X$ 就完全知道 $Y$。
但算协方差:$E(X)=0$, $$\mathrm{Cov}(X,Y)=E(X^3)-E(X)E(X^2)=0-0=0$$ 不相关。
$\rho$ 只度量线性关系。$Y$ 由 $X$ 完全确定,但 $Y=X^2$ 是非线性关系;正、负取值在协方差计算中的贡献恰好抵消。
$\rho=0$ 只表示没有线性相关,不表示「没有关系」。
若 $Y=aX+b$,$a\ne0$ 且 $0<D(X)<\infty$,则 $\mathrm{Cov}=aD(X)$、$D(Y)=a^2D(X)$, $$\rho=\frac{aD(X)}{|a|D(X)}=\mathrm{sgn}(a)$$ 此时 $\rho=\pm1$,符号由斜率决定。
这里讨论 $\sigma_1,\sigma_2>0$、$|\rho|<1$ 的非退化二元正态分布,其联合密度为 $$f(x,y)=\frac{1}{2\pi\sigma_1\sigma_2\sqrt{1-\rho^2}} \exp\left\{-\frac{1}{2(1-\rho^2)}\left[ \Big(\tfrac{x-\mu_1}{\sigma_1}\Big)^2 -2\rho\Big(\tfrac{x-\mu_1}{\sigma_1}\Big)\Big(\tfrac{y-\mu_2}{\sigma_2}\Big) +\Big(\tfrac{y-\mu_2}{\sigma_2}\Big)^2\right]\right\}$$
在联合分布为二元正态这一条件下,不相关与独立等价。仅知道两个边缘分布都为正态,不能应用这一结论。
以下密度公式要求协方差矩阵正定;退化正态分布可能没有这种密度。
记 $\boldsymbol{X},\boldsymbol{\mu}$ 为 $n$ 维列向量,$\boldsymbol{C}$ 为协方差矩阵 ($c_{ij}=\mathrm{Cov}(X_i,X_j)$),以 $\boldsymbol{x}$ 表示取值,则 $$f(x_1,\dots,x_n)=(2\pi)^{-n/2}|\boldsymbol{C}|^{-1/2} \exp\left\{-\tfrac12(\boldsymbol{x}-\boldsymbol{\mu})^{T}\boldsymbol{C}^{-1}(\boldsymbol{x}-\boldsymbol{\mu})\right\}$$
注意第一条的「且相互独立」不能去掉: 各分量都正态,不足以保证整体是多维正态。 第四条是二维那条性质的推广,同样依赖联合正态这一前提。
若 $g$ 是凹函数,$EX$ 与 $E[g(X)]$ 均有限,则 $$E[g(X)]\le g(E[X])$$ 若 $g$ 严格凹,等号成立当且仅当 $X$ 以概率1为常数。对一般凹函数,若 $g$ 在包含 $X$ 取值的区间上为线性函数,也会取等号。
当密度为 $f(x)=2x$($0<x<1$)时,已算得 $E\sqrt X=4/5$、$EX=2/3$。凹函数 $g(x)=\sqrt x$ 给出 $4/5\le\sqrt{2/3}$。
假设 $E(X^2)$、$E(Y^2)$ 均有限。
$$\big[E(XY)\big]^2\le E(X^2)\,E(Y^2)$$把 $X,Y$ 换成离差 $X-EX$、$Y-EY$,立刻得到 $$\mathrm{Cov}(X,Y)^2\le D(X)D(Y)$$ 当两个方差均非零时,即得 $|\rho|\le1$。
期望概括平均水平,方差概括偏离自身期望的程度。相同的期望或方差可以对应不同的分布。
协方差参与和的方差计算;相关系数消除了单位影响。独立推出不相关,反向结论需要额外条件。
使用有限期望、方差和相关系数前,分别确认可积性、有限二阶矩及非零方差。数字特征不能替代完整分布,也不能单独确定决策。
当独立同分布变量的个数增加时,样本平均为何趋近总体均值?它在均值附近的随机波动又呈现什么分布?
大数定律讨论前一个问题;中心极限定理讨论适当标准化后的波动。两者是不同的结论。