跳到主要内容

大阪大学 基礎工学研究科 数理科学 (システム創成専攻) 2024年度 数理科学 [II-8]

Author

祭音Myyura (co-authored with GPT 5.6 SOL)

Description

nn は正の整数とする。実数値確率変数 X1,,XnX_1,\ldots,X_n は独立同分布で、密度は

f(x;μ)=11+μ2+1exp(μ1+μ2+1xx),xR,μRf(x;\mu)=\frac1{\sqrt{1+\mu^2}+1} \exp\left(\frac\mu{\sqrt{1+\mu^2}+1}x-|x|\right),\qquad x\in\mathbb R,\quad\mu\in\mathbb R

とする。Xˉn=n1iXi\bar X_n=n^{-1}\sum_iX_i について、次の問いに理由を付して答えよ。

(1) Xˉn\bar X_nμ\mu の不偏推定量か。

(2) Xˉn2\bar X_n^2μ2\mu^2 の不偏推定量か。

(3) nn\to\infty のとき、Xˉn3\bar X_n^3μ3\mu^3 の一致推定量か。

Kai

s=1+μ2s=\sqrt{1+\mu^2}, a=s+1a=s+1, b=μ/ab=\mu/a とおく。b<1|b|<1, 1b2=2/a1-b^2=2/a であり、密度は x0x\ge0a1e(1b)xa^{-1}e^{-(1-b)x}x<0x<0a1e(1+b)xa^{-1}e^{(1+b)x} と表される。

(1)

E[Xi]=1a{1(1b)21(1+b)2}=4ba(1b2)2=ab=μ.E[X_i]=\frac1a\left\{\frac1{(1-b)^2}-\frac1{(1+b)^2}\right\} =\frac{4b}{a(1-b^2)^2}=ab=\mu.

従って E[Xˉn]=μ\boxed{E[\bar X_n]=\mu} であり、不偏推定量である。

(2)

同様に

E[Xi2]=2a{1(1b)3+1(1+b)3}=2μ2+s+1,E[X_i^2]=\frac2a\left\{\frac1{(1-b)^3}+\frac1{(1+b)^3}\right\} =2\mu^2+s+1,

より Var(Xi)=μ2+s+1=s(s+1)>0\operatorname{Var}(X_i)=\mu^2+s+1=s(s+1)>0。したがって

E[Xˉn2]=μ2+s(s+1)n>μ2\boxed{E[\bar X_n^2]=\mu^2+\frac{s(s+1)}n>\mu^2}

であり、不偏ではない。

(3)

EXi<E|X_i|<\infty なので大数の法則から XˉnPμ\bar X_n\xrightarrow{P}\mu。関数 xx3x\mapsto x^3 は連続であるから Xˉn3Pμ3\boxed{\bar X_n^3\xrightarrow{P}\mu^3}。従って一致推定量である。