間違いしかありません.コメントにてご指摘いただければ幸いです(気が付いた点を特に断りなく頻繁に書き直していますのでご注意ください).

カイ二乗分布の導出

カイ二乗χ2分布の導出

k=1のXの確率密度凾数

X=Z12Z1:F1(X=t)=0tf1(x)dxk=1XF1,f1F1(X=t)=t+tg(z1)dz1Z1=±X,x:0t,z1:0±tg=20tg(z1)dz1(y(x0)),Z10=20tg(x)12x12dxZ1=X(Z10)dz1dx=ddxx=ddxx12=12x12,dz1=12x12dxz1:0t,x:0t=0tg(x)x12dx0tf1(x)dx=0tg(x)x12dx(0t(>0))f1(x)=g(x)x120t(>0)0tt=12πex22x12:g(x)=12πex22=12πex2x12

k=2のXの確率密度凾数

k=2 X=Z12+Z22Z1,Z2:Y=Z12XY=Z22f2(x)=0f1(y)f1(xy)dyk=2Xf2,Z1Z2f1(Z1)f1(Z2)=0x{12πey2y12}{12πe(xy)2(xy)12}dy=12π12π0xey2e(xy)2y12(xy)12dy=12π0xey2(xy)2y12(xy)12dy=12πex20xy12(xy)12dy=12πex2π0a1xaxdx=π=12ex2

k=3のXの確率密度凾数

X=Z12+Z22+Z32Z1,Z2,Z3:Y=Z12+Z22XY=Z32f3(x)=0f2(y)f1(xy)dyk=3Xf3,YZ3f2(Y)f1(Z3)=0x{12ey2}{12πe(xy)2(xy)12}dy=1212π0xey2e(xy)2(xy)12dy=122π0xey2(xy)2(xy)12dy=122πex20x(xy)12dy=122πex22x120a1axdx=2a12=12πex2x12

k=4のXの確率密度凾数

X=Z12+Z22+Z32+Z42Z1,Z2,Z3,Z4:Y=Z12+Z22+Z32XY=Z42f4(x)=0f3(y)f1(xy)dyk=4Xf4,YZ4f3(Y)f1(Z4)=0x{12πey2y12}{12πe(xy)2(xy)12}dy=12π12π0xey2e(xy)2y12(xy)12dy=12π0xey2(xy)2y12(xy)12dy=12πex20xy12(xy)12dy=12πex2xπ20axaxdx=aπ2=14ex2x

k=5のXの確率密度凾数

X=Z12+Z22+Z32+Z42+Z52Z1,Z2,Z3,Z4,Z5:Y=Z12+Z22+Z32+Z42XY=Z52f5(x)=0f4(y)f1(xy)dyk=5Xf5,YZ4f4(Y)f1(Z5)=0x{14ey2y}{12πe(xy)2(xy)12}dy=1412π0xey2e(xy)2y(xy)12dy=142π0xey2(xy)2y(xy)12dy=142πex20xy(xy)12dy=142πex243x320axaxdx=43a32=132πex2x32

k=nXの確率密度凾数の類推

k=1:f1(x)=12πex2x12k=2:f2(x)=12ex2x0k=3:f3(x)=12πex2x12k=4:f4(x)=14ex2x1k=5:f5(x)=132πex2x32k=n:fn(x)=αnex2xn21

係数αnを求める

Fn(t)=0tfn(x)dxnXFn,fnFn(t)=0tαnex2xn21dx=αn0tex2xn21dxαnx,cf(x)dx=cf(x)dx=αn0tet(2u)n212duu=x2,x=2u,dxdu=2,x:0,u:0=αn0tet2n21tn212dt=αn0tet2n2tn21dt=αn2n20tettn21dtFn(t=)=0fn(x)dx=αn2n20ettn21dt=αn2n2Γ(n2)Γ(s)=0etts1dt=αn2n2Γ(n2)=11αn=Fn(t=)2n2Γ(n2)=12n2Γ(n2)

係数αkの確認

α1=1212Γ(12)=1212π=12πα2=1222Γ(22)=1211=12α3=1232Γ(32)=1232π2=12πα4=1242Γ(42)=1221=14α5=1252Γ(52)=12523π4=132π

累積分布凾数の等式よりk=nXの確率密度凾数を求める

0tfn(x)dx=0t12n2Γ(n2)ex2xn21dx(0t(>0))fn(x)=12n2Γ(n2)ex2xn210t(>0)0tt

数学的帰納法による証明

f1(x),fn(x)を認めた上で,fn+1(x)fn(x)nn+1とした式になるかを確認する. f1(x)=12πex2x12fn(x)=12n2Γ(n2)ex2xn21fn+1(x)=12n+12Γ(n+12)ex2xn+121fn(x)nn+1,2 X=Z12++Zn2+Zn+12Z1,,Zn,Zn+1:Y=Z12++ZnXY=Zn+12fn+1(x)=0fn(y)f1(xy)dy=0{12n2Γ(n2)ey2yn21}{12πexy2(xy)12}dy=12n2Γ(n2)12π0ey2exy2yn21(xy)12dycf(x)dx=cf(x)dx=12n2Γ(n2)12π0ey2xy2yn21(xy)12dyABAC=AB+C=12n2Γ(n2)1212π0ex2yn21(xy)12dy=12n2212Γ(n2)1πex20yn21(xy)12dycf(x)dx=cf(x)dx=12n+12Γ(n2)1πex201(xu)n21(xxu)12xduy=xu,dydu=x,x:0,u:01=12n+12Γ(n2)1πex201(xu)n21{x(1u)}12xdu=12n+12Γ(n2)1πex201xn21un21x12(1u)12xdu(AB)C=ACBC=12n+12Γ(n2)1πex2xn21x12x01un21(1u)12ducf(x)dx=cf(x)dx=12n+12Γ(n2)1πex2xn2112+101un21(1u)12duABAC=AB+C=12n+12Γ(n2)1πex2xn+12101un21(1u)12du=12n+12Γ(n2)1πex2xn+121B(n21,12)B(p,q)=01xp(1x)qdx=12n+12Γ(n2)1πex2xn+121(n21)!(12)!(n2112+1)!B(p,q)=p!q!(p+q+1)!=12n+12Γ(n2)1πex2xn+121(n21)!(12+11)!(n2+121)!=12n+12Γ(n2)1πex2xn+121(n21)!(121)!(n+121)!=12n+12Γ(n2)1πex2xn+121Γ(n2)Γ(12)Γ(n+12)Γ(n)=(n1)!=12n+12Γ(n2)1πex2xn+121Γ(n2)Γ(12)Γ(n+12)=12n+121πex2xn+121Γ(12)Γ(n+12)=12n+121πex2xn+121πΓ(n+12)Γ(12)=π=12n+12Γ(n+12)ex2xn+121 任意のn(n1)において上記式が成り立つことを証明できた.
以上によりχ2分布の確率密度凾数が導出された.

0 件のコメント:

コメントを投稿