跳到主要内容

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

Author

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

Description

問題1

X,YX,Y は独立な標準正規変数とする。一辺が a>0a>0 の正方形 [0,a]2[0,a]^2 と、原点を中心に (r,0)(r,0) から (0,r)(0,r) に至る半径 rr の四分円を考える。図のように、円内・正方形外を A,CA,C、正方形内・円外を BB とする。

等面積の正方形と四分円、差領域 A、B、C

(i) ACA\cup C の面積と BB の面積が等しいとする。このとき rraa で表せ。

(ii) (i) のもとで、同時密度 f(x,y)f(x,y) に対し Bf<ACf\iint_Bf<\iint_{A\cup C}f を示し、これから

P(0<X<a,0<Y<a)<14P(X2+Y2<r2)P(0<X<a,0<Y<a)<\frac14P(X^2+Y^2<r^2)

を導け。

(iii) (ii) の右辺を計算し、P(0<X<a)P(0<X<a) の上界を求めよ。

問題2

θ0\theta\ne0 の互いに独立な不偏推定量 X1,,XnX_1,\ldots,X_n が、正の分散 σ12,,σn2\sigma_1^2,\ldots,\sigma_n^2 をもつ。重み付き和 θ^=iλiXi\widehat\theta=\sum_i\lambda_iX_i を考え、iλi=c\sum_i\lambda_i=c とする。

(i) c=1c=1 なら θ^\widehat\theta が不偏であることを示し、その分散を最小化する重みと最小分散を求めよ。

(ii) (i) を用いて σ12,,σn2\sigma_1^2,\ldots,\sigma_n^2 の算術平均が調和平均以上であることを示せ。

(iii) 一般の cc に対し平均二乗誤差 J(c)=E[(θ^θ)2]J(c)=E[(\widehat\theta-\theta)^2] を最小化する重みと、その最小値 Jopt(c)J_{\rm opt}(c) を求めよ。

(iv) Jopt(c)J_{\rm opt}(c) を最小化する cc を求め、不偏の場合より平均二乗誤差が小さくなることを示せ。

题目描述

问题1 X,YX,Y 独立且均为标准正态变量。比较第一象限内的正方形 [0,a]2[0,a]^2a>0a>0)和半径 rr 的四分之一圆盘。图中圆内且正方形外的两部分记为 A,CA,C,正方形内且圆外的部分记为 BB。(i) 若 ACA\cup CBB 面积相等,求 rr。(ii) 证明联合密度满足 Bf<ACf\iint_Bf<\iint_{A\cup C}f,进而证明 P(0<X<a,0<Y<a)<P(X2+Y2<r2)/4P(0<X<a,0<Y<a)<P(X^2+Y^2<r^2)/4。(iii) 计算右边,求 P(0<X<a)P(0<X<a) 的上界。

问题2 参数 θ0\theta\ne0 的独立无偏估计量 XiX_i 有正方差 σi2\sigma_i^2。考虑 θ^=iλiXi\widehat\theta=\sum_i\lambda_iX_i,约束为 iλi=c\sum_i\lambda_i=c。(i) 当 c=1c=1 时证明无偏,求最小方差权重和方差。(ii) 由此证明这些方差的算术平均不小于调和平均。(iii) 对一般 cc 求最小化均方误差的权重和 Jopt(c)J_{\rm opt}(c)。(iv) 再优化 cc,证明所得均方误差小于无偏估计的最小均方误差。

Kai

問題1

(i) 二つの差領域の面積が等しいことは正方形と四分円の面積が等しいことと同値。従って a2=πr2/4a^2=\pi r^2/4 より r=2a/π\boxed{r=2a/\sqrt\pi}

(ii) f(x,y)=(2π)1e(x2+y2)/2f(x,y)=(2\pi)^{-1}e^{-(x^2+y^2)/2} は原点からの距離について厳密に減少する。BB では距離が rr より大きく、ACA\cup C では小さい。同面積なので

Bf<er2/22πB=er2/22πAC<ACf.\iint_Bf<\frac{e^{-r^2/2}}{2\pi}|B| =\frac{e^{-r^2/2}}{2\pi}|A\cup C|<\iint_{A\cup C}f.

共通領域の積分を加え、円対称性を用いると所定の確率不等式を得る。

(iii) 極座標で P(X2+Y2<r2)=0rses2/2ds=1er2/2P(X^2+Y^2<r^2)=\int_0^r se^{-s^2/2}ds=1-e^{-r^2/2}。独立性より左辺は P(0<X<a)2P(0<X<a)^2 なので

P(0<X<a)<121e2a2/π.\boxed{P(0<X<a)<\frac12\sqrt{1-e^{-2a^2/\pi}}}.

問題2

S=iσi2S=\sum_i\sigma_i^{-2} とおく。

(i) E[θ^]=θiλi=θE[\widehat\theta]=\theta\sum_i\lambda_i=\theta。コーシー・シュワルツ不等式から

1=(i(λiσi)σi1)2(iλi2σi2)S.1=\left(\sum_i(\lambda_i\sigma_i)\sigma_i^{-1}\right)^2 \le\left(\sum_i\lambda_i^2\sigma_i^2\right)S.

等号条件より λi=σi2/S\boxed{\lambda_i=\sigma_i^{-2}/S}Vmin=1/S\boxed{V_{\min}=1/S}

(ii) 等重み λi=1/n\lambda_i=1/n の分散は n2iσi21/Sn^{-2}\sum_i\sigma_i^2\ge1/S。従って n1iσi2n/S\boxed{n^{-1}\sum_i\sigma_i^2\ge n/S} となり、算術平均は調和平均以上である。

(iii) 偏りは (c1)θ(c-1)\theta であり、分散の最小化は (i) と同様なので

λi=cσi2S,Jopt(c)=c2S+(c1)2θ2.\boxed{\lambda_i=\frac{c\sigma_i^{-2}}S},\qquad \boxed{J_{\rm opt}(c)=\frac{c^2}S+(c-1)^2\theta^2}.

(iv) cc で微分すると 2c/S+2(c1)θ2=02c/S+2(c-1)\theta^2=0。従って

c=Sθ21+Sθ2,Jopt(c)=θ21+Sθ2<1S=Jopt(1).\boxed{c_* =\frac{S\theta^2}{1+S\theta^2}},\qquad \boxed{J_{\rm opt}(c_*)=\frac{\theta^2}{1+S\theta^2}<\frac1S=J_{\rm opt}(1)}.