官术网_书友最值得收藏!

2.11 連續隨機變量的變換

如果X是服從連續分布F的隨機變量,那么對任意函數gx),Y=gX)仍是隨機變量.Y的分布是什么?

首先考慮支撐.如果X的支撐是X且函數gXY,則Y的支撐為Y.例如,如果X的支撐是[0, 1]且gx)=1+2x,則Y=gX)的支撐是[1, 3].如果X的支撐是R+gx)=log(x),則Y=gX)的支撐為R.

Y的概率函數為FYy)=P[Yy]=P[gX)≤y].令By)是集合{x∈R:gx)≤y}.事件{gX)≤y}和{XBy)}是等價的.所以Y的分布函數為

FYy)=P[XBy)]

因此,Y的分布由X的概率函數確定.

gx)單調遞增時,gx)有反函數

hy)=g-1y

X=hY)且By)=(-∞, hy)].Y的分布函數為

FYy)=P[Xhy)]=FXhy))

它的密度函數是分布函數的導數.根據鏈式法則,可得

最后一個等號成立是因為hy)關于y的導數為正.

現考慮gx)單調遞減的情況,設其反函數為hy),則By)=[hy),∞),

FYy)=P[Xhy)]=1-FXhy))

它的密度函數是分布函數的導數,

最后一個等號成立是因為hy)關于y的導數為負.

gx)嚴格單調時,Y的密度為

fYy)=fXg-1y))Jy

其中

稱為變換的雅可比行列式(Jacobian).這個概念與微積分中類似.

由此得到下述定理.

定理2.4 XfXx),fx)是X上的連續函數.設gx)嚴格單調,g-1y)在Y上連續可微,則對所有的yY,都有

fYy)=fXg-1y))Jy

其中.

定理2.4給出了一個變換Y的密度函數的簡潔公式.下述四個例子具體說明定理2.4.

例13 fXx)=exp(-x), x≥0.設Y=λXλ>0,即gx)=λx.Y的支撐為Y=[0,∞).設函數gx)是單調遞增的,其反函數為hy)=y/λ.雅可比行列式是反函數的導數,

Y的密度函數為

其中y≥0.由于

其中第二個等號通過變量變換x=y/λ得到,故密度是合理的.

例14 fXx)=1,0≤x≤1.設Y=gX),其中gx)=-log(x).由于X的支撐為[0, 1],Y的支撐為Y=(0,∞).函數gx)單調遞減,其反函數為

hy)=g-1y)=exp(-y

hy)求導得到雅可比行列式:

注意,fXg-1y))=1,y≥0.計算Y的密度,

fYy)=fXg-1y))Jy)=exp(-y

其中y≥0.這是指數分布的密度函數.也就是說,如果X服從均勻分布,則Y=-log(X)服從指數分布.

例15 X有任意連續可逆(嚴格遞增)的累積分布函數FXx).定義隨機變量Y=FXx).Y的支撐為Y=[0, 1]:則在[0,1]上Y的累積分布函數為

求導得概率密度函數:

這是服從U[0, 1]隨機變量的密度函數.故YU[0, 1].

變換Y=FXx)被稱為概率積分變換(probability integral transformation).無論初始X的分布是什么,該公式都能將Y變換為均勻分布,這個事實是相當奇妙的.

例16fXx)表示圖2-7中時薪的密度函數.令Y=log(X).如果X的支撐為R+,則Y的支撐為R,其反函數為hy)=exp(y),雅可比行列式為exp(y).Y的密度為fYy)=fX(exp(y))exp(y).如圖2-8所示,該密度函數比時薪的密度更對稱,偏度更小.

圖2-8 時薪對數的密度函數

主站蜘蛛池模板: 容城县| 耒阳市| 安远县| 保靖县| 常宁市| 新宾| 芦溪县| 叙永县| 郎溪县| 澜沧| 冀州市| 金昌市| 岑巩县| 彭水| 舒兰市| 商水县| 于都县| 铜梁县| 墨江| 平顶山市| 申扎县| 右玉县| 大同市| 越西县| 渭源县| 孝义市| 绥江县| 泰顺县| 文安县| 司法| 图木舒克市| 招远市| 简阳市| 济源市| 新民市| 县级市| 玛多县| 宁远县| 新竹县| 交口县| 淮阳县|