跳到主要内容

北海道大学 情報科学院 情報科学専攻 情報理工学コース 2022年8月実施 専門科目1 問3 (確率・統計)

Author

祭音Myyura (with GPT 5.5)

Description

確率・統計に関する以下の問いに答えよ。 ただし、答えだけでなく、導出過程もわかるように解答すること。

[1]

連続型確率変数 XX の分布関数が

F(x)={0,(x<0),kx3,(0x<2),1,(2x)F(x)= \begin{cases} 0, & (x<0), \\ kx^3, & (0\leq x<2), \\ 1, & (2\leq x) \end{cases}

で与えられるとき、以下の問いに答えよ。

(1) 実数 kk の値を求めよ。

(2) 確率 P(0<X1)P(0<X\leq 1) を求めよ。

(3) 確率変数 XX の確率密度関数 f(x)f(x) を求めよ。

(4) 関数 ff のグラフの概形を示せ。

(5) 確率変数 XX の期待値 E[X]E[X] と分散 V[X]V[X] を求めよ。

[2]

母平均が μ\mu、母分散が σ2\sigma^2 の母集団からの大きさ nn の無作為標本を
X1,,XnX_1,\ldots,X_n とするとき、以下の問いに答えよ。
ただし、μ,σ2\mu,\sigma^2 はともに未知数であるとする。

(1) 標本平均を

Xˉ=1ni=1nXi\bar{X}=\frac{1}{n}\sum_{i=1}^{n}X_i

とするとき、Xˉ\bar{X}μ\mu の不偏推定量であること、すなわち、

E[Xˉ]=μE[\bar{X}]=\mu

が成り立つことを示せ。

(2) 等式

E[(Xˉ)2]=σ2n+μ2E[(\bar{X})^2]=\frac{\sigma^2}{n}+\mu^2

が成り立つことを示せ。

(3) 統計量

T1=1n1i=1n(XiXˉ)2T_1=\frac{1}{n-1}\sum_{i=1}^{n}(X_i-\bar{X})^2

σ2\sigma^2 の不偏推定量であることを示せ。

(4) 統計量

T2=(Xˉ)21n(n1)i=1n(XiXˉ)2T_2=(\bar{X})^2-\frac{1}{n(n-1)}\sum_{i=1}^{n}(X_i-\bar{X})^2

μ2\mu^2 の不偏推定量であることを示せ。

(5) 母集団が正規分布 N(μ,σ2)N(\mu,\sigma^2) であるとき、確率 P(T2<0)P(T_2<0) は正であること、すなわち、

P(T2<0)>0P(T_2<0)>0

が成り立つことを示せ。
ただし、X1,,XnX_1,\ldots,X_n が正規分布 N(μ,σ2)N(\mu,\sigma^2) からの無作為標本であるとき、
Xˉ\bar{X}T1T_1 は独立な確率変数となることを用いてよい。

Kai

[1]

(1)

x=2x=2 における F(x)F(x) の連続性より、

limx2F(x)=F(2)k23=1  k=18\begin{aligned} \lim_{x \to 2-} F(x) &= F(2) \\ k \cdot 2^3 &= 1 \\ \therefore \ \ k &= \frac{1}{8} \end{aligned}

がわかる。

(2)

P(0<X1)=F(1)F(0)=18 \begin{aligned} P \left( 0 \lt X \leq 1 \right) &= F(1) - F(0) \\ &= \frac{1}{8} \end{aligned}

(3)

f(x)=ddxF(x)={0(x<0)38x2(0x<2)0(22)\begin{aligned} f(x) &= \frac{d}{dx} F(x) \\ &= \begin{cases} 0 &(x \lt 0) \\ \frac{3}{8} x^2 &(0 \leq x \lt 2) \\ 0 &(2 \leq 2) \end{cases} \end{aligned}

(4)

omitted

(5)

E[X]=xf(x)dx=3802x3dx=38[x44]02=32E[X2]=x2f(x)dx        (X2 の期待値 )=3802x4dx=38[x55]02=125Var[X]=E[X2]E[X]2=320 \begin{aligned} \mathrm{E} \left[ X \right] &= \int_{-\infty}^\infty x f(x) dx \\ &= \frac{3}{8} \int_0^2 x^3 dx \\ &= \frac{3}{8} \left[ \frac{x^4}{4} \right]_0^2 \\ &= \frac{3}{2} \\ \mathrm{E} \left[ X^2 \right] &= \int_{-\infty}^\infty x^2 f(x) dx \ \ \ \ \ \ \ \ \left( X^2 \text{ の期待値 } \right) \\ &= \frac{3}{8} \int_0^2 x^4 dx \\ &= \frac{3}{8} \left[ \frac{x^5}{5} \right]_0^2 \\ &= \frac{12}{5} \\ \mathrm{Var} \left[ X \right] &= \mathrm{E} \left[ X^2 \right] - \mathrm{E} \left[ X \right]^2 \\ &= \frac{3}{20} \end{aligned}

[2]

X1,X2,,XnX_1,X_2,\ldots,X_n は、母平均 μ\mu、母分散 σ2\sigma^2 の母集団からの大きさ nn の無作為標本である。
したがって、

E[Xi]=μ,V[Xi]=σ2E[X_i]=\mu,\qquad V[X_i]=\sigma^2

であり、

E[Xi2]=V[Xi]+{E[Xi]}2=σ2+μ2E[X_i^2]=V[X_i]+\{E[X_i]\}^2=\sigma^2+\mu^2

が成り立つ。

(1)

E[Xˉ]=1ni=1nE[Xi]=μ \begin{aligned} \mathrm{E} \left[ \bar{X} \right] &= \frac{1}{n} \sum_{i=1}^n \mathrm{E} \left[ X_i \right] \\ &= \mu \end{aligned}

(2)

E[(Xˉ)2]=1n2i=1nj=1nE[XiXj]=1n2(i=1nE[Xi2]+i,j (ij)E[XiXj])=1n2(i=1nE[Xi2]+i,j (ij)E[Xi]E[Xj])=1n2(n(σ2+μ2)+(n2n)μ2)=σ2n+μ2 \begin{aligned} \mathrm{E} \left[ \left( \bar{X} \right)^2 \right] &= \frac{1}{n^2} \sum_{i=1}^n \sum_{j=1}^n \mathrm{E} \left[ X_i X_j \right] \\ &= \frac{1}{n^2} \left( \sum_{i=1}^n \mathrm{E} \left[ X_i^2 \right] + \sum_{i,j \ (i \ne j)} \mathrm{E} \left[ X_i X_j \right] \right) \\ &= \frac{1}{n^2} \left( \sum_{i=1}^n \mathrm{E} \left[ X_i^2 \right] + \sum_{i,j \ (i \ne j)} \mathrm{E} \left[ X_i \right] \mathrm{E} \left[ X_j \right] \right) \\ &= \frac{1}{n^2} \left( n \left( \sigma^2 + \mu^2 \right) + \left( n^2 - n \right) \mu^2 \right) \\ &= \frac{\sigma^2}{n} + \mu^2 \end{aligned}

(3)

T1=1n1i=1n(XiXˉ)2=1n1i=1n(Xi22XˉXi+(Xˉ)2)=1n1(i=1nXi2n(Xˉ)2)=1n1i=1nXi2nn1(Xˉ)2  E[T1]=1n1i=1nE[Xi2]nn1E[(Xˉ)2]=nn1(σ2+μ2)nn1(σ2n+μ2)=σ2 \begin{aligned} T_1 &= \frac{1}{n-1} \sum_{i=1}^n \left( X_i - \bar{X} \right)^2 \\ &= \frac{1}{n-1} \sum_{i=1}^n \left( X_i^2 - 2 \bar{X} X_i + \left( \bar{X} \right)^2 \right) \\ &= \frac{1}{n-1} \left( \sum_{i=1}^n X_i^2 - n \left( \bar{X} \right)^2 \right) \\ &= \frac{1}{n-1} \sum_{i=1}^n X_i^2 - \frac{n}{n-1} \left( \bar{X} \right)^2 \\ \therefore \ \ \mathrm{E} \left[ T_1 \right] &= \frac{1}{n-1} \sum_{i=1}^n \mathrm{E} \left[ X_i^2 \right] - \frac{n}{n-1} \mathrm{E} \left[ \left( \bar{X} \right)^2 \right] \\ &= \frac{n}{n-1} \left( \sigma^2 + \mu^2 \right) - \frac{n}{n-1} \left( \frac{\sigma^2}{n} + \mu^2 \right) \\ &= \sigma^2 \end{aligned}

(4)

T2=(Xˉ)21nT1  E[T2]=E[(Xˉ)2]1nE[T1]=σ2n+μ21nσ2=μ2 \begin{aligned} T_2 &= \left( \bar{X} \right)^2 - \frac{1}{n} T_1 \\ \therefore \ \ \mathrm{E} \left[ T_2 \right] &= \mathrm{E} \left[ \left( \bar{X} \right)^2 \right] - \frac{1}{n} \mathrm{E} \left[ T_1 \right] \\ &= \frac{\sigma^2}{n} + \mu^2 - \frac{1}{n} \sigma^2 \\ &= \mu^2 \end{aligned}

(5)

母集団が正規分布 N(μ,σ2)N(\mu,\sigma^2) であるとする。

(4) より、

T2=(Xˉ)21nT1T_2=(\bar{X})^2-\frac{1}{n}T_1

である。したがって、

T2<0T_2<0

となるためには、

(Xˉ)2<1nT1(\bar{X})^2<\frac{1}{n}T_1

が成り立てばよい。

任意の正の数 a>0a>0 をとる。
このとき、

Xˉ<a|\bar{X}|<a

かつ

T1>na2T_1>na^2

が成り立てば、

(Xˉ)2<a2(\bar{X})^2<a^2

かつ

1nT1>a2\frac{1}{n}T_1>a^2

となるので、

(Xˉ)2<1nT1(\bar{X})^2<\frac{1}{n}T_1

である。したがって、

T2<0T_2<0

が成り立つ。

よって、

{Xˉ<a, T1>na2}{T2<0}\{|\bar{X}|<a,\ T_1>na^2\} \subset \{T_2<0\}

であるから、

P(T2<0)P(Xˉ<a, T1>na2)P(T_2<0) \geq P(|\bar{X}|<a,\ T_1>na^2)

となる。

ここで、母集団が正規分布 N(μ,σ2)N(\mu,\sigma^2) であるから、

XˉN(μ,σ2n)\bar{X}\sim N\left(\mu,\frac{\sigma^2}{n}\right)

である。正規分布は任意の区間に正の確率をもつので、

P(Xˉ<a)>0P(|\bar{X}|<a)>0

である。

また、正規母集団からの標本について、

(n1)T1σ2χn12\frac{(n-1)T_1}{\sigma^2}\sim \chi^2_{n-1}

である。χ2\chi^2 分布は正の方向に広がりをもつので、

P(T1>na2)>0P(T_1>na^2)>0

である。

さらに、問題文より、Xˉ\bar{X}T1T_1 は独立である。したがって、

P(Xˉ<a, T1>na2)=P(Xˉ<a)P(T1>na2)>0\begin{aligned} P(|\bar{X}|<a,\ T_1>na^2) &= P(|\bar{X}|<a)P(T_1>na^2) \\ &>0 \end{aligned}

となる。

以上より、

P(T2<0)P(Xˉ<a, T1>na2)>0P(T_2<0)\geq P(|\bar{X}|<a,\ T_1>na^2)>0

であるから、

P(T2<0)>0P(T_2<0)>0

が示された。