跳到主要内容

大阪大学 基礎工学研究科 数理科学 (システム創成専攻) 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+1x−∣x∣),x∈R,μ∈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=n−1∑iXi\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) n→∞n\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, 1−b2=2/a1-b^2=2/a であり、密度は x≥0x\ge0 で a−1e−(1−b)xa^{-1}e^{-(1-b)x}、x<0x<0 で a−1e(1+b)xa^{-1}e^{(1+b)x} と表される。

(1)​

E[Xi]=1a{1(1−b)2−1(1+b)2}=4ba(1−b2)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(1−b)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)​

E∣Xi∣<∞E|X_i|<\infty なので大数の法則から Xˉn→Pμ\bar X_n\xrightarrow{P}\mu。関数 x↦x3x\mapsto x^3 は連続であるから Xˉn3→Pμ3\boxed{\bar X_n^3\xrightarrow{P}\mu^3}。従って一致推定量である。