跳到主要内容

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

Author

AKIRA (小红书:94184092292)

Description

問題1

確率変数 XX は確率密度関数

f(x;θ)={θeθx,x>00,x0f(x; \theta) = \begin{cases} \theta e^{-\theta x}, & x > 0 \\ 0, & x \le 0 \end{cases}

の指数分布にしたがう。ただし θ>0\theta > 0 はパラメータである。 ある定数 θ0>0\theta_0 > 0 に対して、帰無仮説 H0:θ=θ0H_0 : \theta = \theta_0、対立仮説 H1:θ<θ0H_1 : \theta < \theta_0 の仮説検定を有意水準 α (0<α<1)\alpha \ (0 < \alpha < 1) で行いたい。 以下の設問に答えなさい。その導出過程も示すこと。

(1) 定数 b>0b > 0 を定めておき、X>bX > b のとき H0H_0 を棄却する。 この仮説検定の有意水準が α\alpha となるような定数 bb を求めよ。

(2) 定数 c>0c > 0 を定めておき、パラメータ θ\theta の信頼区間を

S(x)={θ0<θcx}S(x) = \left\{ \theta \mid 0 < \theta \le \frac{c}{x} \right\}

とする。P(θS(X))=1αP(\theta \in S(X)) = 1 - \alpha を満たすような定数 cc を求めよ。

(3) 与えられた定数 d>0d > 0 に対して、事象 {X>d}\{ X > d \} を条件とする XX の条件付き分布にしたがう確率変数 YY を定義する。 すなわち、任意の y>dy > d に対して P(Y>y)=P(X>yX>d)P(Y > y) = P(X > y \mid X > d) である。 定数 b>db' > d を定めておき、Y>bY > b' のとき H0H_0 を棄却する。この仮説検定の有意水準が α\alpha となるような定数 bb' を求めよ。

(4) 設問(3)の確率変数 YY からパラメータ θ\theta の信頼区間をつくりたい。 ある関数 h(y)h(y) を用いて、

T(y)={θ0<θh(y)}T(y) = \{ \theta \mid 0 < \theta \le h(y) \}

とする。P(θT(Y))=1αP(\theta \in T(Y)) = 1 - \alpha を満たすような関数 h(y)h(y) を求めよ。

問題2

以下の設問に答えなさい。 ただし、N(μ,σ2)N(\mu, \sigma^2) は平均 μ\mu、分散 σ2\sigma^2 の正規分布、E[]E[\cdot] は期待値を表す。

(1) 独立ではないが無相関であるような実確率変数の組 (X,Y)(X, Y) の例を1つ挙げよ。 またそれが独立ではないこと、無相関であることの証明も示せ。

(2) XXYY を独立な N(0,1)N(0, 1) に従う確率変数とする。 Z=X/YZ = X / Y が従う確率分布の確率密度関数を求めよ。

以下の設問では、次のように定義される関数 Γ(a)\Gamma(a)

Γ(a)=0xa1exdx\Gamma(a) = \int_0^\infty x^{a-1} e^{-x} dx

および a,b>0a, b > 0 なるパラメータを持つ確率密度関数 g(x;a,b)g(x; a, b)

g(x;a,b)={(baΓ(a))1xa1exb,x>00,x0g(x; a, b) = \begin{cases} (b^a \Gamma(a))^{-1} x^{a-1} e^{-\frac{x}{b}}, &x > 0 \\ 0, &x \le 0 \end{cases}

を用いる。 確率密度関数 g(x;a,b)g(x; a, b) を持つ確率分布を G(a,b)G(a, b) とする。 また必要であれば a>0a > 0 に対して Γ(a+1)=aΓ(a)\Gamma(a+1) = a\Gamma(a) を用いてよい。

(3) G(a,b)G(a, b) に従う確率変数 XX について、そのモーメント母関数

MX(r)=E[erX]M_X(r) = E[e^{rX}]

が有限であるような実数 rr の条件を示し、そのときの MX(r)M_X(r) の値を求めよ。

(4) 確率変数 X1,X2,,XnX_1, X_2, \ldots, X_n は独立で N(0,v)N(0, v) に従うものとする。 このとき Y=i=1nXi2Y = \sum_{i=1}^n X_i^2 は、あるパラメータ a1,b1a_1, b_1 を持つ G(a1,b1)G(a_1, b_1) に従うことを示すとともに、a1,b1a_1, b_1v,nv, n を用いて表せ。

(5) G(a,b)G(a, b) に従う確率変数 XX の期待値 μ\mu と分散 σ2\sigma^2a,ba, b を用いて表せ。

(6) x1,x2,,xnx_1, x_2, \ldots, x_n をパラメータ a,ba, b が未知である G(a,b)G(a, b) からの無作為標本とする。 設問(5)の結果とモーメント法を用いて、a,ba, b に対する推定値 a^,b^\hat{a}, \hat{b} を標本平均

xˉ=1ni=1nxi,s2=1ni=1n(xixˉ)2\bar{x} = \frac{1}{n} \sum_{i=1}^n x_i, \quad s^2 = \frac{1}{n} \sum_{i=1}^n (x_i - \bar{x})^2

の関数として表せ。ただしモーメント法とは、パラメータを KK 個持つ確率密度関数 f(x;θ1,θ2,,θK)f(x; \theta_1, \theta_2, \ldots, \theta_K) のモーメント

mk(θ1,θ2,,θK)=xkf(x;θ1,θ2,,θK)dxm_k(\theta_1, \theta_2, \ldots, \theta_K) = \int_{-\infty}^\infty x^k f(x; \theta_1, \theta_2, \ldots, \theta_K) dx

を標本モーメント

m^k=1ni=1nxik\hat{m}_k = \frac{1}{n} \sum_{i=1}^n x_i^k

と等しいと置き、m^k=mk(θ1,θ2,,θK)\hat{m}_k = m_k(\theta_1, \theta_2, \ldots, \theta_K) (k=1,,K)(k = 1, \ldots, K) なる KK 個の連立方程式を θ1,θ2,,θK\theta_1, \theta_2, \ldots, \theta_K について解くことで、推定値 θ1,θ2,,θK\theta_1, \theta_2, \ldots, \theta_K を得る方法である。

Kai

問題1

問題2