跳到主要内容

京都大学 情報学研究科 システム科学専攻 2013年8月実施 専門科目 確率統計

Author​

犬 (finalized by 祭音Myyura with assistance from GPT 6 Astra)

Description​

問題1​

既知の母平均 μ0\mu_0 と未知の母分散 v>0v>0 をもつ正規母集団 N(μ0,v)N(\mu_0,v) から無作為標本 X1,…,XnX_1,\ldots,X_n を得る。母集団の密度は

fG(x;v)=12πvexp⁡(−(x−μ0)22v).f_G(x;v)=\frac1{\sqrt{2\pi v}}\exp\left(-\frac{(x-\mu_0)^2}{2v}\right).

(1) モーメント母関数 MX(t)=E[etX]M_X(t)=E[e^{tX}] を求めよ。

(2) (1) を用いて、母平均周りの 33 次と 44 次のモーメントを求めよ。kk 次モーメントは E[(X−μ0)k]E[(X-\mu_0)^k] をいう。

(3) vv の最尤推定量 v^\widehat v を求めよ。

(4) E[v^]E[\widehat v] と V[v^]V[\widehat v] を求めよ。

(5) v^\widehat v が不偏推定量であることを示せ。

問題2​

(1) X,Y,ZX,Y,Z は互いに独立で、それぞれ (0,1)(0,1) 上の一様分布に従う。XYXY と Z2Z^2 の同時密度、および P(XY<Z2)P(XY<Z^2) を求めよ。

(2) (X,Y)(X,Y) の同時密度は

f(x,y)=12π1−ρ2exp⁡[−x2−2ρxy+y22(1−ρ2)],0≤ρ<1f(x,y)=\frac1{2\pi\sqrt{1-\rho^2}}\exp\left[-\frac{x^2-2\rho xy+y^2}{2(1-\rho^2)}\right],\qquad0\le\rho<1

である。XX と Z=(Y−ρX)/1−ρ2Z=(Y-\rho X)/\sqrt{1-\rho^2} は互いに独立で、それぞれ N(0,1)N(0,1) に従うことを示せ。

(3) (2) の結果を用いて、(2) の分布における P(X>0,Y>0)P(X>0,Y>0) を ρ\rho で表せ。

题目描述​

问题1 从均值 μ0\mu_0 已知、方差 v>0v>0 未知的正态总体 N(μ0,v)N(\mu_0,v) 抽取随机样本 X1,…,XnX_1,\ldots,X_n;其密度为 (2πv)−1/2exp⁡[−(x−μ0)2/(2v)](2\pi v)^{-1/2}\exp[-(x-\mu_0)^2/(2v)]。(1) 求矩母函数 MX(t)=E[etX]M_X(t)=E[e^{tX}]。(2) 利用 (1) 求三阶、四阶中心矩 E[(X−μ0)3]E[(X-\mu_0)^3]、E[(X−μ0)4]E[(X-\mu_0)^4]。(3) 求 vv 的最大似然估计量 v^\widehat v。(4) 求其期望和方差。(5) 证明其无偏。

问题2 (1) X,Y,ZX,Y,Z 独立且均服从 U(0,1)U(0,1),求 XYXY 与 Z2Z^2 的联合密度及 P(XY<Z2)P(XY<Z^2)。(2) 若 (X,Y)(X,Y) 的密度为

f(x,y)=exp⁡[−(x2−2ρxy+y2)/(2(1−ρ2))]2π1−ρ2,0≤ρ<1,f(x,y)=\frac{\exp[-(x^2-2\rho xy+y^2)/(2(1-\rho^2))]}{2\pi\sqrt{1-\rho^2}},\quad0\le\rho<1,

证明 XX 与 Z=(Y−ρX)/1−ρ2Z=(Y-\rho X)/\sqrt{1-\rho^2} 独立且均为标准正态变量。(3) 利用 (2) 将 P(X>0,Y>0)P(X>0,Y>0) 表为 ρ\rho 的函数。

Kai​

問題1​

(1) 指数を平方完成して積分すると MX(t)=exp⁡(μ0t+vt2/2)\boxed{M_X(t)=\exp(\mu_0t+vt^2/2)}。

(2) MX−μ0(t)=evt2/2=1+vt2/2+v2t4/8+⋯M_{X-\mu_0}(t)=e^{vt^2/2}=1+vt^2/2+v^2t^4/8+\cdots より、E[(X−μ0)3]=0\boxed{E[(X-\mu_0)^3]=0}、E[(X−μ0)4]=3v2\boxed{E[(X-\mu_0)^4]=3v^2}。

(3) 対数尤度を微分して

ℓ′(v)=−n2v+12v2∑i(Xi−μ0)2=0,v^=1n∑i(Xi−μ0)2.\ell'(v)=-\frac n{2v}+\frac1{2v^2}\sum_i(X_i-\mu_0)^2=0, \qquad\boxed{\widehat v=\frac1n\sum_i(X_i-\mu_0)^2}.

(4), (5) 独立性と (2) より

E[v^]=v,V[v^]=1n2∑i(3v2−v2)=2v2n.\boxed{E[\widehat v]=v},\qquad \boxed{V[\widehat v]=\frac1{n^2}\sum_i(3v^2-v^2)=\frac{2v^2}n}.

従って不偏である。

問題2​

(1)​

U=XY,V=Z2U=XY,V=Z^2 とおく。0<u,v<10<u,v<1 で

fU(u)=∫u1dxx=−log⁡u,fV(v)=12v.f_U(u)=\int_u^1\frac{dx}{x}=-\log u,\qquad f_V(v)=\frac1{2\sqrt v}.

独立性から、同時密度は fU,V(u,v)=−log⁡u/(2v)\boxed{f_{U,V}(u,v)=-\log u/(2\sqrt v)}(単位正方形内)、それ以外は零である。従って

P(U<V)=∫0112v∫0v(−log⁡u)du dv=12∫01v1/2(1−log⁡v)dv=59.P(U<V)=\int_0^1\frac1{2\sqrt v}\int_0^v(-\log u)du\,dv =\frac12\int_0^1v^{1/2}(1-\log v)dv =\boxed{\frac59}.

(2)​

y=ρx+1−ρ2zy=\rho x+\sqrt{1-\rho^2}z のヤコビアンは 1−ρ2\sqrt{1-\rho^2}。変換後は

fX,Z(x,z)=12πe−(x2+z2)/2=e−x2/22πe−z2/22π.f_{X,Z}(x,z)=\frac1{2\pi}e^{-(x^2+z^2)/2} =\frac{e^{-x^2/2}}{\sqrt{2\pi}}\frac{e^{-z^2/2}}{\sqrt{2\pi}}.

従って独立な標準正規変数である。

(3)​

X>0,Y>0X>0,Y>0 は x>0, z>−ρx/1−ρ2x>0,\ z>-\rho x/\sqrt{1-\rho^2} に対応する。(X,Z)(X,Z) の分布は回転対称であり、領域の角度は π/2+arcsin⁡ρ\pi/2+\arcsin\rho。従って

P(X>0,Y>0)=14+arcsin⁡ρ2π.\boxed{P(X>0,Y>0)=\frac14+\frac{\arcsin\rho}{2\pi}}.