113 年 國立臺北大學統計學系碩士班《數理統計(含機率概論)》

📄 試題原卷 免費註冊後即可對照原始考卷 PDF免費註冊

第 I.1 題15 分

Let the joint pdf of random variables X and Y be given by

f(x,y)={cxif 0<x<y<10elsewhere.f(x,y) = \begin{cases} cx & \text{if } 0 < x < y < 1 \\ 0 & \text{elsewhere.} \end{cases}

Let Z = X + Y. Derive the CDF of Z.

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

本題考查:

  1. 聯合機率密度函數的正規化。

  2. 由聯合密度求函數 Z=X+YZ=X+Y 的累積分布函數:

    FZ(z)=P(Z≤z)=P(X+Y≤z)F_Z(z)=P(Z\le z)=P(X+Y\le z)

  3. 依照直線 x+y=zx+y=z 與原始定義域的交集,分段計算二重積分。

原始定義域為

0<x<y<10<x<y<1

這是在單位正方形中,由 y=xy=x、y=1y=1 與 x=0x=0 所圍成的三角形。


第一步:求常數 cc

由聯合機率密度函數的總積分必須等於 11:

∫01∫x1cx dy dx=1\int_0^1\int_x^1 cx\,dy\,dx=1

因此

c∫01x(1−x) dx=1c\int_0^1 x(1-x)\,dx=1 c(12−13)=1c\left(\frac12-\frac13\right)=1 c6=1\frac{c}{6}=1

所以

c=6c=6

故聯合密度為

f(x,y)=6x,0<x<y<1f(x,y)=6x,\qquad 0<x<y<1

第二步:建立 FZ(z)F_Z(z) 的積分區域

由定義,

FZ(z)=P(X+Y≤z)F_Z(z)=P(X+Y\le z)

因此需在區域

0<x<y<1,x+y≤z0<x<y<1,\qquad x+y\le z

上積分:

FZ(z)=∬0<x<y<1x+y≤z6x dy dxF_Z(z)=\iint_{\substack{0<x<y<1\\x+y\le z}}6x\,dy\,dx

由於 X>0X>0 且 Y>XY>X,所以 Z=X+Y>0Z=X+Y>0;又因為 X<1X<1、Y<1Y<1,所以 Z<2Z<2。因此只需討論 0<z<20<z<2。

直線 x+y=zx+y=z 與 y=xy=x 的交點滿足

x=z−xx=z-x

故

x=z2x=\frac z2

這會造成 z=1z=1 前後的積分區域不同,因此必須分成兩段。


情形一:0<z<10<z<1

此時直線 x+y=zx+y=z 完全位於 y=1y=1 以下。

對固定的 xx 而言,yy 的範圍為

x<y<z−xx<y<z-x

且必須滿足

x<z−xx<z-x

因此

0<x<z20<x<\frac z2

所以

FZ(z)=∫0z/2∫xz−x6x dy dxF_Z(z) =\int_0^{z/2}\int_x^{z-x}6x\,dy\,dx

先對 yy 積分:

FZ(z)=6∫0z/2x[(z−x)−x] dxF_Z(z) =6\int_0^{z/2}x\bigl[(z-x)-x\bigr]\,dx =6∫0z/2x(z−2x) dx=6\int_0^{z/2}x(z-2x)\,dx

計算得

FZ(z)=6[zx22−2x33]0z/2F_Z(z) =6\left[\frac{zx^2}{2}-\frac{2x^3}{3}\right]_0^{z/2} =z34=\frac{z^3}{4}

情形二:1≤z<21\le z<2

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 I.2 題15 分

Let the joint pdf of random variables X and Y be given by

f(x,y)={6xif 0<x<y<10elsewhere.f(x,y) = \begin{cases} 6x & \text{if } 0 < x < y < 1 \\ 0 & \text{elsewhere.} \end{cases}

Let Z = X + Y. Compute E(X|Z).

登入後即可作答並保存紀錄。

這一題的完整詳解

本題要求計算條件期望值 E(X∣Z)E(X|Z)。我們需要先求出 X 在給定 Z 下的條件機率密度函數 (conditional pdf),然後計算其期望值。

步驟一:計算 X 的邊際機率密度函數 (Marginal pdf of X)
fX(x)=∫−∞∞f(x,y) dyf_X(x) = \int_{-\infty}^{\infty} f(x,y) \, dy
對於 0<x<10 < x < 1,積分範圍為 x<y<1x < y < 1。
fX(x)=∫x16x dy=6x[y]x1=6x(1−x)for 0<x<1f_X(x) = \int_{x}^{1} 6x \, dy = 6x [y]_{x}^{1} = 6x(1-x) \quad \text{for } 0 < x < 1
fX(x)=0elsewhere.f_X(x) = 0 \quad \text{elsewhere.}

步驟二:計算 Z 的機率密度函數 (pdf of Z)
由上一題的 CDF FZ(z)F_Z(z),我們可以求得 Z 的 pdf fZ(z)f_Z(z)。
對於 0<z≤10 < z \le 1: fZ(z)=FZ′(z)=ddz(z34)=3z24f_Z(z) = F_Z'(z) = \frac{d}{dz} \left( \frac{z^3}{4} \right) = \frac{3z^2}{4}。
對於 1<z<21 < z < 2: fZ(z)=FZ′(z)=ddz(z34−(z−1)3)=3z24−3(z−1)2f_Z(z) = F_Z'(z) = \frac{d}{dz} \left( \frac{z^3}{4} - (z-1)^3 \right) = \frac{3z^2}{4} - 3(z-1)^2。

fZ(z)={3z24if 0<z≤13z24−3(z−1)2if 1<z<20elsewhere.f_Z(z) = \begin{cases} \frac{3z^2}{4} & \text{if } 0 < z \le 1 \\ \frac{3z^2}{4} - 3(z-1)^2 & \text{if } 1 < z < 2 \\ 0 & \text{elsewhere.} \end{cases}

步驟三:計算條件機率密度函數 fX∣Z(x∣z)f_{X|Z}(x|z)
fX∣Z(x∣z)=f(x,z) for X=x,Z=x+Y=zfZ(z)f_{X|Z}(x|z) = \frac{f(x,z) \text{ for } X=x, Z=x+Y=z}{f_Z(z)}
這裡的 f(x,z)f(x,z) 指的是 X 和 Z 的聯合機率密度函數。然而,我們已知 X 和 Y 的聯合 pdf f(x,y)f(x,y),且 Z=X+YZ = X+Y。
更直接的定義是:
fX∣Z(x∣z)=f(x,y=z−x)fZ(z)f_{X|Z}(x|z) = \frac{f(x, y=z-x)}{f_Z(z)}
我們需要確定 X 的取值範圍,在給定 z 的條件下。
我們知道 0<x<y<10 < x < y < 1 且 x+y=zx+y = z。
將 y=z−xy = z-x 代入 x<y<1x < y < 1:
x<z−x  ⟹  2x<z  ⟹  x<z/2x < z-x \implies 2x < z \implies x < z/2。
z−x<1  ⟹  z−1<xz-x < 1 \implies z-1 < x。
所以,對於給定的 z,x 的範圍是 max⁡(0,z−1)<x<z/2\max(0, z-1) < x < z/2。

因此,
fX∣Z(x∣z)=6xfZ(z)for max⁡(0,z−1)<x<z/2f_{X|Z}(x|z) = \frac{6x}{f_Z(z)} \quad \text{for } \max(0, z-1) < x < z/2

我們需要分情況討論 z 的範圍:

  • 情況 1:0<z≤10 < z \le 1
    此時,z−1≤0z-1 \le 0,所以 max⁡(0,z−1)=0\max(0, z-1) = 0。
    x 的範圍是 0<x<z/20 < x < z/2。
    fZ(z)=3z24f_Z(z) = \frac{3z^2}{4}。
    fX∣Z(x∣z)=6x3z2/4=24x3z2=8xz2for 0<x<z/2f_{X|Z}(x|z) = \frac{6x}{3z^2/4} = \frac{24x}{3z^2} = \frac{8x}{z^2} \quad \text{for } 0 < x < z/2
    驗證積分為 1:
    ∫0z/28xz2 dx=8z2[x22]0z/2=8z2((z/2)22)=8z2z2/42=8z2z28=1\int_{0}^{z/2} \frac{8x}{z^2} \, dx = \frac{8}{z^2} \left[ \frac{x^2}{2} \right]_{0}^{z/2} = \frac{8}{z^2} \left( \frac{(z/2)^2}{2} \right) = \frac{8}{z^2} \frac{z^2/4}{2} = \frac{8}{z^2} \frac{z^2}{8} = 1
    正確。

  • 情況 2:1<z<21 < z < 2
    此時,z−1>0z-1 > 0,所以 max⁡(0,z−1)=z−1\max(0, z-1) = z-1。
    x 的範圍是 z−1<x<z/2z-1 < x < z/2。
    fZ(z)=3z24−3(z−1)2=34(z2−4(z−1)2)=34(z2−4(z2−2z+1))=34(−3z2+8z−4)f_Z(z) = \frac{3z^2}{4} - 3(z-1)^2 = \frac{3}{4} (z^2 - 4(z-1)^2) = \frac{3}{4} (z^2 - 4(z^2 - 2z + 1)) = \frac{3}{4} (-3z^2 + 8z - 4)。

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 I.3 題10 分

Let the joint pdf of random variables X and Y be given by

f(x,y)={6xif 0<x<y<10elsewhere.f(x,y) = \begin{cases} 6x & \text{if } 0 < x < y < 1 \\ 0 & \text{elsewhere.} \end{cases}

Compute E[E(X|Z)] based on your answer in 2.

登入後即可作答並保存紀錄。

這一題的完整詳解

本題要求計算 E[E(X∣Z)]E[E(X|Z)]。根據全期望定理 (Law of Total Expectation),我們知道 E[E(X∣Z)]=E(X)E[E(X|Z)] = E(X)。因此,我們可以利用上一題的結果 E(X∣Z=z)E(X|Z=z) 來計算 E[E(X∣Z)]E[E(X|Z)],或者直接計算 E(X)E(X)。直接計算 E(X)E(X) 通常比較簡便。

方法一:利用全期望定理 E[E(X∣Z)]=E(X)E[E(X|Z)] = E(X)

首先,我們需要計算 X 的邊際機率密度函數 (marginal pdf of X)。
fX(x)=∫−∞∞f(x,y) dyf_X(x) = \int_{-\infty}^{\infty} f(x,y) \, dy
對於 0<x<10 < x < 1,積分範圍為 x<y<1x < y < 1。
fX(x)=∫x16x dy=6x[y]x1=6x(1−x)for 0<x<1f_X(x) = \int_{x}^{1} 6x \, dy = 6x [y]_{x}^{1} = 6x(1-x) \quad \text{for } 0 < x < 1
fX(x)=0elsewhere.f_X(x) = 0 \quad \text{elsewhere.}
現在計算 E(X)E(X):
E(X)=∫−∞∞xfX(x) dx=∫01x⋅6x(1−x) dx=∫01(6x2−6x3) dxE(X) = \int_{-\infty}^{\infty} x f_X(x) \, dx = \int_{0}^{1} x \cdot 6x(1-x) \, dx = \int_{0}^{1} (6x^2 - 6x^3) \, dx
E(X)=[6x33−6x44]01=[2x3−32x4]01E(X) = \left[ \frac{6x^3}{3} - \frac{6x^4}{4} \right]_{0}^{1} = \left[ 2x^3 - \frac{3}{2}x^4 \right]_{0}^{1}
E(X)=(2(1)3−32(1)4)−(0)=2−32=42−32=12E(X) = (2(1)^3 - \frac{3}{2}(1)^4) - (0) = 2 - \frac{3}{2} = \frac{4}{2} - \frac{3}{2} = \frac{1}{2}
所以,E[E(X∣Z)]=E(X)=12E[E(X|Z)] = E(X) = \frac{1}{2}。

方法二:利用上一題的 E(X∣Z=z)E(X|Z=z) 結果
根據定義,E[E(X∣Z)]=∫−∞∞E(X∣Z=z)fZ(z) dzE[E(X|Z)] = \int_{-\infty}^{\infty} E(X|Z=z) f_Z(z) \, dz。
我們需要分段積分:
E[E(X∣Z)]=∫01E(X∣Z=z)fZ(z) dz+∫12E(X∣Z=z)fZ(z) dzE[E(X|Z)] = \int_{0}^{1} E(X|Z=z) f_Z(z) \, dz + \int_{1}^{2} E(X|Z=z) f_Z(z) \, dz
其中,E(X∣Z=z)E(X|Z=z) 和 fZ(z)f_Z(z) 的表達式來自前兩題。

第一部分:∫01(z3)(3z24) dz\int_{0}^{1} \left(\frac{z}{3}\right) \left(\frac{3z^2}{4}\right) \, dz
∫01z3⋅3z24 dz=∫01z34 dz=[z416]01=116\int_{0}^{1} \frac{z}{3} \cdot \frac{3z^2}{4} \, dz = \int_{0}^{1} \frac{z^3}{4} \, dz = \left[ \frac{z^4}{16} \right]_{0}^{1} = \frac{1}{16}

第二部分:

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 I.4 題5 分

Let the joint pdf of random variables X and Y be given by

f(x,y)={6xif 0<x<y<10elsewhere.f(x,y) = \begin{cases} 6x & \text{if } 0 < x < y < 1 \\ 0 & \text{elsewhere.} \end{cases}

Compute E(X) based on the pdf of X.

登入後即可作答並保存紀錄。

這一題的完整詳解

本題要求根據 X 的機率密度函數 (pdf) 來計算 X 的期望值 E(X)E(X)。我們需要先求出 X 的邊際機率密度函數,然後再計算其期望值。

步驟一:計算 X 的邊際機率密度函數 fX(x)f_X(x)
聯合機率密度函數為 f(x,y)=6xf(x,y) = 6x 在 0<x<y<10 < x < y < 1 的區域。
fX(x)=∫−∞∞f(x,y) dyf_X(x) = \int_{-\infty}^{\infty} f(x,y) \, dy
對於 0<x<10 < x < 1,積分範圍是 x<y<1x < y < 1。
fX(x)=∫x16x dyf_X(x) = \int_{x}^{1} 6x \, dy
fX(x)=6x[y]x1f_X(x) = 6x [y]_{x}^{1}
fX(x)=6x(1−x)f_X(x) = 6x (1-x)
因此,X 的邊際機率密度函數為:

fX(x)={6x(1−x)if 0<x<10elsewhere.f_X(x) = \begin{cases} 6x(1-x) & \text{if } 0 < x < 1 \\ 0 & \text{elsewhere.} \end{cases}
🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 I.5 題5 分

Suppose you are dealing with an unknown joint pdf f(x, y). Suppose you only have two random samples of size n: one is X1, X2, ..., Xn from the distribution of X, and the other is W1, W2, ..., Wn from the distribution of W, where W = E(X|Z). Can you determine which random sample would allow you to derive a more precise estimator of E(X)? Provide an explanation for your choice.

登入後即可作答並保存紀錄。

這一題的完整詳解

本題旨在探討,在不知道聯合機率密度函數 f(x,y)f(x,y) 的情況下,我們有兩組獨立同分佈 (i.i.d.) 的隨機樣本:一組是來自隨機變數 X 的樣本 X1,…,XnX_1, \dots, X_n,另一組是來自隨機變數 W 的樣本 W1,…,WnW_1, \dots, W_n,其中 W=E(X∣Z)W = E(X|Z)。我們要判斷哪一個樣本集(X 或 W)能提供對 E(X)E(X) 更精確的估計,並給出解釋。

核心觀念:
估計的精確度通常與估計量的變異數 (variance) 有關。變異數越小,估計越精確。對於來自同一分佈的樣本,其算術平均數是該分佈期望值的無偏估計量,而其變異數的大小決定了估計的精確度。

分析:
我們有兩個樣本集:

  1. 樣本集 1:{X1,X2,…,Xn}\{X_1, X_2, \dots, X_n\},來自 X 的分佈。我們想估計 E(X)E(X)。一個自然的估計量是樣本平均數 Xˉ=1n∑i=1nXi\bar{X} = \frac{1}{n} \sum_{i=1}^n X_i。
  2. 樣本集 2:{W1,W2,…,Wn}\{W_1, W_2, \dots, W_n\},來自 W 的分佈,其中 W=E(X∣Z)W = E(X|Z)。我們想估計 E(X)E(X)。一個自然的估計量是樣本平均數 Wˉ=1n∑i=1nWi\bar{W} = \frac{1}{n} \sum_{i=1}^n W_i。

由於 W=E(X∣Z)W = E(X|Z),這意味著 W 是 X 的一個「條件期望值」。通常情況下,條件期望值會比原始隨機變數「更穩定」或「更集中」。

數學解釋:
我們需要比較 Xˉ\bar{X} 和 Wˉ\bar{W} 作為 E(X)E(X) 的估計量的精確度。精確度由估計量的變異數來衡量。
Var(Xˉ)=Var(1n∑Xi)=1n2∑Var(Xi)=1nVar(X)Var(\bar{X}) = Var(\frac{1}{n} \sum X_i) = \frac{1}{n^2} \sum Var(X_i) = \frac{1}{n} Var(X) (假設 Var(X)Var(X) 存在)。
Var(Wˉ)=Var(1n∑Wi)=1n2∑Var(Wi)=1nVar(W)Var(\bar{W}) = Var(\frac{1}{n} \sum W_i) = \frac{1}{n^2} \sum Var(W_i) = \frac{1}{n} Var(W) (假設 Var(W)Var(W) 存在)。

我們知道,對於任何隨機變數 X 和 Z,全期望定理告訴我們 E(X)=E[E(X∣Z)]E(X) = E[E(X|Z)]。
由於 W=E(X∣Z)W = E(X|Z),所以 E(W)=E[E(X∣Z)]=E(X)E(W) = E[E(X|Z)] = E(X)。這表示 Wˉ\bar{W} 和 Xˉ\bar{X} 都是 E(X)E(X) 的無偏估計量。

現在我們比較 Var(X)Var(X) 和 Var(W)Var(W)。
我們有 Var(X)=E[Var(X∣Z)]+Var[E(X∣Z)]Var(X) = E[Var(X|Z)] + Var[E(X|Z)]。
由於 W=E(X∣Z)W = E(X|Z),所以 Var(X)=E[Var(X∣Z)]+Var(W)Var(X) = E[Var(X|Z)] + Var(W)。
因為 Var(X∣Z)Var(X|Z) 是條件變異數,它是一個非負的隨機變數。所以 E[Var(X∣Z)]≥0E[Var(X|Z)] \ge 0。
因此, Var(X)≥Var(W)Var(X) \ge Var(W)。

如果 Var(X∣Z)Var(X|Z) 在 Z 的分佈上不是恆為零(也就是說,對於某些 Z 的值,X 的條件變異數不為零),那麼 E[Var(X∣Z)]>0E[Var(X|Z)] > 0,這將導致 Var(X)>Var(W)Var(X) > Var(W)。

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 II.1.a 題5 分

Let X1,X2,…,XnX_1, X_2, \dots, X_n be a random sample from Bernoulli(p). Consider two estimators, p^U=Xˉ\hat{p}_U = \bar{X} and p^B=Xˉ+12n+1\hat{p}_B = \frac{\bar{X} + \frac{1}{2}}{\sqrt{n}+1}. Find the limiting distribution of n(p^U−p)\sqrt{n}(\hat{p}_U - p) and n(p^B−p)\sqrt{n}(\hat{p}_B - p), respectively.

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

本題考查:

  1. 樣本平均數的中央極限定理。
  2. 估計量的偏誤與一致性。
  3. 判斷標準化後的隨機變數是否具有有限的極限分布。

對於 Xi∼Bernoulli⁡(p)X_i\sim\operatorname{Bernoulli}(p),

E(Xi)=p,Var⁡(Xi)=p(1−p).E(X_i)=p,\qquad \operatorname{Var}(X_i)=p(1-p).

當 0<p<10<p<1 時,由中央極限定理,

n(Xˉ−p)→dN(0,p(1−p)).\sqrt{n}(\bar X-p)\xrightarrow{d}N\bigl(0,p(1-p)\bigr).

一、估計量 p^U=Xˉ\hat p_U=\bar X

直接套用中央極限定理:

n(p^U−p)=n(Xˉ−p).\sqrt{n}(\hat p_U-p) = \sqrt{n}(\bar X-p).

因此,

n(p^U−p)→dN(0,p(1−p))\boxed{ \sqrt{n}(\hat p_U-p) \xrightarrow{d} N\bigl(0,p(1-p)\bigr) }

二、估計量 p^B\hat p_B

給定

p^B=Xˉ+12n+1.\hat p_B=\frac{\bar X+\frac12}{\sqrt n+1}.

考慮其極限:

p^B=Xˉ+12n+1.\hat p_B = \frac{\bar X+\frac12}{\sqrt n+1}.

由大數法則,Xˉ→pp\bar X\xrightarrow{p}p,而分母 n+1→∞\sqrt n+1\to\infty,所以

p^B→p0.\hat p_B\xrightarrow{p}0.

因此,當 p>0p>0 時,p^B\hat p_B 並不一致估計 pp。進一步計算標準化形式:

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 II.1.b 題5 分

Let X1,X2,…,XnX_1, X_2, \dots, X_n be a random sample from Bernoulli(p). Consider two estimators, p^U=Xˉ\hat{p}_U = \bar{X} and p^B=Xˉ+12n+1\hat{p}_B = \frac{\bar{X} + \frac{1}{2}}{\sqrt{n}+1}. Show that p^U\hat{p}_U and p^B\hat{p}_B are the consistent estimators of p.

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

估計量 p^n\hat p_n 為參數 pp 的一致估計量,定義為

p^n→Pp,\hat p_n \xrightarrow{P} p,

亦即對任意 ε>0\varepsilon>0,

lim⁡n→∞P(∣p^n−p∣>ε)=0.\lim_{n\to\infty} P\left(\left|\hat p_n-p\right|>\varepsilon\right)=0.

本題使用大數法則:

Xˉ=1n∑i=1nXi→PE(Xi)=p,\bar X=\frac{1}{n}\sum_{i=1}^n X_i \xrightarrow{P} E(X_i)=p,

因為 Xi∼Bernoulli⁡(p)X_i\sim\operatorname{Bernoulli}(p)。


一、證明 p^U=Xˉ\hat p_U=\bar X 為一致估計量

由於

E(Xi)=p,Var⁡(Xi)=p(1−p),E(X_i)=p,\qquad \operatorname{Var}(X_i)=p(1-p),

可得

E(Xˉ)=p,Var⁡(Xˉ)=p(1−p)n.E(\bar X)=p, \qquad \operatorname{Var}(\bar X) =\frac{p(1-p)}{n}.

由切比雪夫不等式,對任意 ε>0\varepsilon>0,

P(∣Xˉ−p∣≥ε)≤Var⁡(Xˉ)ε2=p(1−p)nε2.P\left(\left|\bar X-p\right|\ge \varepsilon\right) \le \frac{\operatorname{Var}(\bar X)}{\varepsilon^2} = \frac{p(1-p)}{n\varepsilon^2}.

當 n→∞n\to\infty 時,

p(1−p)nε2→0.\frac{p(1-p)}{n\varepsilon^2}\to 0.

因此

P(∣Xˉ−p∣≥ε)→0,P\left(\left|\bar X-p\right|\ge \varepsilon\right)\to 0,

所以

p^U=Xˉ→Pp.\boxed{\hat p_U=\bar X\xrightarrow{P}p}.

故 p^U\hat p_U 是 pp 的一致估計量。


二、檢查 p^B\hat p_B 是否一致

題目給定

p^B=Xˉ+12n+1.\hat p_B = \frac{\bar X+\frac12}{\sqrt n+1}.

由於 0≤Xˉ≤10\le \bar X\le 1,因此

0≤p^B≤32n+1.0\le \hat p_B \le \frac{\frac32}{\sqrt n+1}.

右側隨 n→∞n\to\infty 趨近於 00,故由夾擠定理,

p^B→P0.\hat p_B\xrightarrow{P}0.

也可利用 Xˉ→Pp\bar X\xrightarrow{P}p 寫成

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 II.1.c 題5 分

Let X1,X2,…,XnX_1, X_2, \dots, X_n be a random sample from Bernoulli(p). Consider two estimators, p^U=Xˉ\hat{p}_U = \bar{X} and p^B=Xˉ+12n+1\hat{p}_B = \frac{\bar{X} + \frac{1}{2}}{\sqrt{n}+1}. Find a consistent estimator of 1p\frac{1}{p}.

登入後即可作答並保存紀錄。

這一題的完整詳解

本題要求我們尋找一個參數 1p\frac{1}{p} 的一致估計量。我們已經知道 Xˉ\bar{X} 是 pp 的一個一致估計量,且 E(Xˉ)=pE(\bar{X}) = p,Var(Xˉ)=p(1−p)nVar(\bar{X}) = \frac{p(1-p)}{n}。

核心觀念:
如果 θ^n\hat{\theta}_n 是 θ\theta 的一致估計量,且 gg 是一個在 θ\theta 的鄰域內連續的函數,那麼 g(θ^n)g(\hat{\theta}_n) 通常是 g(θ)g(\theta) 的一致估計量。

尋找 1p\frac{1}{p} 的估計量:
由於 Xˉ\bar{X} 是一致估計量,且 lim⁡n→∞E(Xˉ)=p\lim_{n \to \infty} E(\bar{X}) = p 且 lim⁡n→∞Var(Xˉ)=0\lim_{n \to \infty} Var(\bar{X}) = 0。
考慮函數 g(x)=1xg(x) = \frac{1}{x}。我們想找到 g(p)=1pg(p) = \frac{1}{p} 的一致估計量。
一個自然的嘗試是使用 g(Xˉ)=1Xˉg(\bar{X}) = \frac{1}{\bar{X}}。

證明 1Xˉ\frac{1}{\bar{X}} 是一致估計量:
我們需要證明 lim⁡n→∞P(∣1Xˉ−1p∣<ϵ)=1\lim_{n \to \infty} P\left(\left|\frac{1}{\bar{X}} - \frac{1}{p}\right| < \epsilon\right) = 1。
這等價於證明 1Xˉ→P1p\frac{1}{\bar{X}} \xrightarrow{P} \frac{1}{p}。
根據 Slutsky 定理,如果 Xn→PcX_n \xrightarrow{P} c,其中 cc 是一個常數,且 gg 在 cc 點連續,則 g(Xn)→Pg(c)g(X_n) \xrightarrow{P} g(c)。
我們知道 Xˉ→Pp\bar{X} \xrightarrow{P} p (根據大數法則)。
如果 p∈(0,1)p \in (0, 1),則函數 g(x)=1xg(x) = \frac{1}{x} 在 x=px=p 點是連續的。
因此,1Xˉ→P1p\frac{1}{\bar{X}} \xrightarrow{P} \frac{1}{p}。
這意味著 1Xˉ\frac{1}{\bar{X}} 是 1p\frac{1}{p} 的一致估計量,只要 p∈(0,1)p \in (0, 1)。

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 II.1.d 題5 分

Let X1,X2,…,XnX_1, X_2, \dots, X_n be a random sample from Bernoulli(p). Consider two estimators, p^U=Xˉ\hat{p}_U = \bar{X} and p^B=Xˉ+12n+1\hat{p}_B = \frac{\bar{X} + \frac{1}{2}}{\sqrt{n}+1}. Use n=1n = 1 to show that no unbiased estimator of 1p\frac{1}{p} exists.

登入後即可作答並保存紀錄。

這一題的完整詳解

本題要求使用 n=1n=1 的情況,證明不存在 1p\frac{1}{p} 的無偏估計量 (unbiased estimator)。

無偏估計量的定義:
一個估計量 θ^\hat{\theta} 被稱為參數 θ\theta 的無偏估計量,如果 E(θ^)=θE(\hat{\theta}) = \theta。

考慮 n=1n=1 的情況:
當 n=1n=1,我們只有一個樣本 X1X_1。
X1∼Bernoulli(p)X_1 \sim \text{Bernoulli}(p),所以 P(X1=1)=pP(X_1=1) = p 且 P(X1=0)=1−pP(X_1=0) = 1-p。
E(X1)=1⋅p+0⋅(1−p)=pE(X_1) = 1 \cdot p + 0 \cdot (1-p) = p。
Var(X1)=p(1−p)Var(X_1) = p(1-p)。

我們需要證明不存在任何函數 T(X1)T(X_1) 使得 E[T(X1)]=1pE[T(X_1)] = \frac{1}{p} 對於所有的 p∈(0,1)p \in (0, 1) 都成立。

考慮我們已知的估計量:
在上一題中,我們找到了 1Xˉ\frac{1}{\bar{X}} 作為 1p\frac{1}{p} 的一致估計量(當 p∈(0,1]p \in (0, 1])。
當 n=1n=1,Xˉ=X1\bar{X} = X_1。所以,我們考慮估計量 1X1\frac{1}{X_1}。

計算 E[1X1]E[\frac{1}{X_1}]:
1X1\frac{1}{X_1} 的定義域為 X1∈{0,1}X_1 \in \{0, 1\}。
當 X1=1X_1 = 1,1X1=1\frac{1}{X_1} = 1。
當 X1=0X_1 = 0,1X1\frac{1}{X_1} 是無窮大,其期望值是未定義的。
即使我們忽略 X1=0X_1=0 的情況,只考慮 X1=1X_1=1 的情況,其期望值是:
E[1X1∣X1=1]=1E[\frac{1}{X_1} | X_1=1] = 1。
然而,我們需要 E[1X1]=1pE[\frac{1}{X_1}] = \frac{1}{p} 對於所有的 p∈(0,1)p \in (0, 1)。
如果 p=1p=1,則 X1X_1 必定是 1,所以 E[1X1]=1=11E[\frac{1}{X_1}] = 1 = \frac{1}{1}。這似乎是滿足的。
但是,如果 p<1p < 1,則 X1X_1 有可能為 0。
P(X1=0)=1−p>0P(X_1 = 0) = 1-p > 0。
在這種情況下,E[1X1]E[\frac{1}{X_1}] 的計算涉及到除以零,其期望值是無窮大,或者在嚴格的機率論定義下,這個期望值是未定義的。

嚴格證明:
假設存在一個 n=1n=1 的無偏估計量 T(X1)T(X_1),使得 E[T(X1)]=1pE[T(X_1)] = \frac{1}{p} 對於所有的 p∈(0,1)p \in (0, 1) 成立。
E[T(X1)]=T(1)P(X1=1)+T(0)P(X1=0)E[T(X_1)] = T(1) P(X_1=1) + T(0) P(X_1=0)
E[T(X1)]=T(1)⋅p+T(0)⋅(1−p)E[T(X_1)] = T(1) \cdot p + T(0) \cdot (1-p)

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 II.2.a 題5 分

Let X1,X2,…,XnX_1, X_2, \dots, X_n be a random sample from U(θ,θ+1)U(\theta, \theta+1). Find the shortest length 100(1−α)%100(1 - \alpha)\% confidence interval for θ\theta in the class C(X):=[X(1)a,X(n)b]C(X) := \left[\frac{X_{(1)}}{a}, \frac{X_{(n)}}{b}\right].

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

本題給定的信賴區間類別為

C(X)=[X(1)a,X(n)b],C(X)=\left[\frac{X_{(1)}}{a},\frac{X_{(n)}}{b}\right],

其中 a,ba,b 應為與樣本無關的常數。要成為 θ\theta 的信賴區間,必須使覆蓋率

Pθ(X(1)a≤θ≤X(n)b)P_\theta\left(\frac{X_{(1)}}{a}\leq \theta\leq \frac{X_{(n)}}{b}\right)

等於或至少達到 1−α1-\alpha,且應與未知參數 θ\theta 無關。

然而,本題所給的「除以 a,ba,b」形式會使覆蓋率依賴 θ\theta,因此無法形成固定信賴水準的信賴區間。

覆蓋率分析

令

Yi=Xi−θ,i=1,…,n.Y_i=X_i-\theta,\qquad i=1,\dots,n.

則 Y1,…,YnY_1,\dots,Y_n 為來自 U(0,1)U(0,1) 的隨機樣本,且

X(1)=θ+Y(1),X(n)=θ+Y(n).X_{(1)}=\theta+Y_{(1)},\qquad X_{(n)}=\theta+Y_{(n)}.

覆蓋事件為

θ+Y(1)a≤θ且θ+Y(n)b≥θ.\frac{\theta+Y_{(1)}}{a}\leq\theta \quad\text{且}\quad \frac{\theta+Y_{(n)}}{b}\geq\theta.

若假設 a,b>0a,b>0,上式等價於

Y(1)≤(a−1)θ,Y(n)≥(b−1)θ.Y_{(1)}\leq (a-1)\theta, \qquad Y_{(n)}\geq (b-1)\theta.

可見覆蓋事件直接含有 θ\theta,所以其機率一般會隨 θ\theta 改變。

例如,當 a>1a>1 時,

Pθ(X(1)a≤θ)=P(Y(1)≤(a−1)θ).P_\theta\left(\frac{X_{(1)}}{a}\leq\theta\right) = P\left(Y_{(1)}\leq(a-1)\theta\right).

在 0<(a−1)θ<10<(a-1)\theta<1 時,由於

P(Y(1)>t)=(1−t)n,P(Y_{(1)}>t)=(1-t)^n,

因此

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 II.2.b 題5 分

Let X1,X2,…,XnX_1, X_2, \dots, X_n be a random sample from U(θ,θ+1)U(\theta, \theta+1). Prove that P(θ′∈C(X))<1−αP(\theta' \in C(X)) < 1 - \alpha for θ′≠θ\theta' \neq \theta, where C(X):=[X(1)a,X(n)b]C(X) := \left[\frac{X_{(1)}}{a}, \frac{X_{(n)}}{b}\right] is the confidence interval found in part (a).

登入後即可作答並保存紀錄。

這一題的完整詳解

本題要求證明,對於一個真實值 θ′≠θ\theta' \neq \theta,我們在 part (a) 中找到的信賴區間 C(X):=[X(1)a,X(n)b]C(X) := \left[\frac{X_{(1)}}{a}, \frac{X_{(n)}}{b}\right],包含 θ′\theta' 的機率 P(θ′∈C(X))P(\theta' \in C(X)) 小於 1−α1-\alpha。
這實際上是要求證明信賴區間的覆蓋機率。對於一個置信水準為 1−α1-\alpha 的信賴區間 [L(X),U(X)][L(X), U(X)],我們要求 P(θ∈[L(X),U(X)])=1−αP(\theta \in [L(X), U(X)]) = 1-\alpha 對於所有 θ\theta 成立。
題目要求證明 P(θ′∈C(X))<1−αP(\theta' \in C(X)) < 1 - \alpha for θ′≠θ\theta' \neq \theta。這句話的表述有些不尋常。通常,信賴區間的定義是 P(θ∈C(X))=1−αP(\theta \in C(X)) = 1-\alpha 對於所有 θ\theta 成立。
如果 C(X)C(X) 是一個 100(1−α)%100(1-\alpha)\% 的信賴區間,那麼對於任何 θ0\theta_0,我們有 P(θ0∈C(X))=1−αP(\theta_0 \in C(X)) = 1-\alpha。
這裡 θ′\theta' 是某個特定的值,而不是隨機變數。

重新理解題目:
這句話的意思可能是:如果我們固定了 aa 和 bb,使得對於某個 θ0\theta_0,有 P(θ0∈C(X))=1−αP(\theta_0 \in C(X)) = 1-\alpha,那麼對於另一個 θ′≠θ0\theta' \neq \theta_0,我們有 P(θ′∈C(X))<1−αP(\theta' \in C(X)) < 1-\alpha。
這是在問信賴區間的覆蓋機率是否「精確」為 1−α1-\alpha 或「大於」1−α1-\alpha(稱為過度覆蓋,over-coverage)。

假設 part (a) 找到了這樣一組 a,ba, b 使得 P(θ∈C(X))=1−αP(\theta \in C(X)) = 1-\alpha 對於所有 θ\theta 成立。
那麼,對於任何一個特定的值 θ′\theta',我們將 C(X)C(X) 中的 θ\theta 替換為 θ′\theta'。
P(θ′∈C(X))=P(X(1)a≤θ′≤X(n)b)P(\theta' \in C(X)) = P\left(\frac{X_{(1)}}{a} \le \theta' \le \frac{X_{(n)}}{b}\right)。
由於 X(1)X_{(1)} 和 X(n)X_{(n)} 是基於 θ\theta 的隨機變數,它們的分佈取決於 θ\theta。
P(θ′∈C(X))=P(X(1)a≤θ′ and X(n)b≥θ′)P(\theta' \in C(X)) = P\left(\frac{X_{(1)}}{a} \le \theta' \text{ and } \frac{X_{(n)}}{b} \ge \theta'\right)。
=P(X(1)≤aθ′ and X(n)≥bθ′)= P\left(X_{(1)} \le a\theta' \text{ and } X_{(n)} \ge b\theta'\right)。

令 Y(1)=X(1)−θY_{(1)} = X_{(1)} - \theta 和 Y(n)=X(n)−θY_{(n)} = X_{(n)} - \theta。
X(1)=Y(1)+θX_{(1)} = Y_{(1)} + \theta
X(n)=Y(n)+θX_{(n)} = Y_{(n)} + \theta
P(Y(1)+θ≤aθ′ and Y(n)+θ≥bθ′)P(Y_{(1)} + \theta \le a\theta' \text{ and } Y_{(n)} + \theta \ge b\theta')
=P(Y(1)≤aθ′−θ and Y(n)≥bθ′−θ)= P(Y_{(1)} \le a\theta' - \theta \text{ and } Y_{(n)} \ge b\theta' - \theta)。

如果 part (a) 找到了 a,ba, b 使得 P(θ∈C(X))=1−αP(\theta \in C(X)) = 1-\alpha 對所有 θ\theta 成立,這意味著 aa 和 bb 必須是常數。
如果 aa 和 bb 是常數,那麼 P(θ′∈C(X))P(\theta' \in C(X)) 的計算將會依賴於 θ′\theta'。
P(θ′∈C(X))=P(X(1)a≤θ′≤X(n)b)P(\theta' \in C(X)) = P\left(\frac{X_{(1)}}{a} \le \theta' \le \frac{X_{(n)}}{b}\right)。

核心問題:
信賴區間的覆蓋機率 P(θ∈C(X))P(\theta \in C(X)) 應該是 1−α1-\alpha 對於所有 θ\theta。
如果 part (a) 確實找到了這樣的 C(X)C(X),那麼 P(θ′∈C(X))P(\theta' \in C(X)) 對於任何 θ′\theta' 的值,其意義何在?

一種可能的解釋:
「Confidence interval」的定義是 P(θ∈C(X))=1−αP(\theta \in C(\mathbf{X})) = 1-\alpha。
如果 part (a) 找到了這樣的 C(X)C(X),那麼對於任何 θ′\theta',這個機率定義就是 P(θ′∈C(X))=1−αP(\theta' \in C(X)) = 1-\alpha。
這與題目要求證明 P(θ′∈C(X))<1−αP(\theta' \in C(X)) < 1 - \alpha 相矛盾。

另一種解釋:
也許 part (a) 找到的區間是「最短長度」的,但可能不是「精確」的 1−α1-\alpha 覆蓋。
例如,對於某些分佈,最短長度的 CI 可能會導致過度覆蓋 (over-coverage),即 P(θ∈C(X))≥1−αP(\theta \in C(X)) \ge 1-\alpha,且對於某些 θ\theta,嚴格大於 1−α1-\alpha。

如果 part (a) 確實找到了 1−α1-\alpha 的 CI:
那麼 P(θ∈C(X))=1−αP(\theta \in C(X)) = 1-\alpha。
題目要求證明 P(θ′∈C(X))<1−αP(\theta' \in C(X)) < 1 - \alpha for θ′≠θ\theta' \neq \theta。
這句話的結構是:對於一個固定的 θ′\theta',其機率小於 1−α1-\alpha。
這意味著,對於某個 θ′\theta',我們計算 P(C(X) contains θ′)P(C(X) \text{ contains } \theta'),這個機率小於 1−α1-\alpha。

讓我們假設 part (a) 找到了 a,ba, b 使得 P(θ∈C(X))=1−αP(\theta \in C(X)) = 1-\alpha 對於所有 θ\theta 成立。
那麼,對於任何一個特定的值 θ′\theta', 我們計算 P(θ′∈C(X))P(\theta' \in C(X))。
P(θ′∈C(X))=P(X(1)a≤θ′≤X(n)b)P(\theta' \in C(X)) = P\left(\frac{X_{(1)}}{a} \le \theta' \le \frac{X_{(n)}}{b}\right).
令 Y(1)=X(1)−θY_{(1)} = X_{(1)}-\theta 和 Y(n)=X(n)−θY_{(n)} = X_{(n)}-\theta。
P(Y(1)+θ≤aθ′ and Y(n)+θ≥bθ′)P(Y_{(1)}+\theta \le a\theta' \text{ and } Y_{(n)}+\theta \ge b\theta').
P(Y(1)≤aθ′−θ and Y(n)≥bθ′−θ)P(Y_{(1)} \le a\theta'-\theta \text{ and } Y_{(n)} \ge b\theta'-\theta).

如果 a,ba, b 是常數,那麼 aθ′−θa\theta'-\theta 和 bθ′−θb\theta'-\theta 會隨著 θ\theta 的變化而變化。
這意味著 P(θ′∈C(X))P(\theta' \in C(X)) 是一個關於 θ\theta 的函數。

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 II.3.a 題5 分

Let X1,X2,…,XnX_1, X_2, \dots, X_n be a random sample from the pdf f(x∣θ)=θx−2,0<θ≤x<∞f(x|\theta) = \theta x^{-2}, 0 < \theta \le x < \infty. Find a sufficient statistic of θ\theta.

登入後即可作答並保存紀錄。

這一題的完整詳解

本題要求找出參數 θ\theta 的一個充分統計量 (sufficient statistic)。

充分統計量的定義:
一個統計量 T(X1,…,Xn)T(X_1, \dots, X_n) 被稱為參數 θ\theta 的充分統計量,如果給定 T=tT = t,樣本 X1,…,XnX_1, \dots, X_n 的條件機率分佈不依賴於 θ\theta。
即,f(x1,…,xn∣T(x1,…,xn)=t,θ)=f(x1,…,xn∣T(x1,…,xn)=t)f(x_1, \dots, x_n | T(x_1, \dots, x_n) = t, \theta) = f(x_1, \dots, x_n | T(x_1, \dots, x_n) = t)。

因子分解定理 (Factorization Theorem):
一個統計量 T(X)T(X) 是 θ\theta 的充分統計量,若且唯若聯合機率密度函數 f(x1,…,xn∣θ)f(x_1, \dots, x_n | \theta) 可以分解為:
f(x1,…,xn∣θ)=g(T(x1,…,xn),θ)⋅h(x1,…,xn)f(x_1, \dots, x_n | \theta) = g(T(x_1, \dots, x_n), \theta) \cdot h(x_1, \dots, x_n)
其中 gg 依賴於 θ\theta 和 TT,而 hh 僅僅依賴於樣本 x1,…,xnx_1, \dots, x_n。

步驟一:寫出聯合機率密度函數 (Joint pdf)
樣本 X1,…,XnX_1, \dots, X_n 是 i.i.d. 且來自 f(x∣θ)=θx−2f(x|\theta) = \theta x^{-2},其中 0<θ≤x<∞0 < \theta \le x < \infty。
聯合 pdf 為:
f(x1,…,xn∣θ)=∏i=1nf(xi∣θ)f(x_1, \dots, x_n | \theta) = \prod_{i=1}^n f(x_i|\theta)
f(x1,…,xn∣θ)=∏i=1n(θxi−2)for 0<θ≤xi<∞ for all if(x_1, \dots, x_n | \theta) = \prod_{i=1}^n (\theta x_i^{-2}) \quad \text{for } 0 < \theta \le x_i < \infty \text{ for all } i
f(x1,…,xn∣θ)=θn(∏i=1nxi−2)for θ≤min⁡(x1,…,xn)<∞f(x_1, \dots, x_n | \theta) = \theta^n \left( \prod_{i=1}^n x_i^{-2} \right) \quad \text{for } \theta \le \min(x_1, \dots, x_n) < \infty
這個條件 0<θ≤xi<∞0 < \theta \le x_i < \infty 對於所有 ii 成立,意味著 θ≤min⁡(x1,…,xn)\theta \le \min(x_1, \dots, x_n)。
令 X(1)=min⁡(x1,…,xn)X_{(1)} = \min(x_1, \dots, x_n)。
所以,
f(x1,…,xn∣θ)=θn(∏i=1nxi−2)for θ≤X(1)<∞f(x_1, \dots, x_n | \theta) = \theta^n \left( \prod_{i=1}^n x_i^{-2} \right) \quad \text{for } \theta \le X_{(1)} < \infty
並且在其他情況下為 0。

步驟二:應用因子分解定理
我們將聯合 pdf 分解為 g(T,θ)⋅h(x1,…,xn)g(T, \theta) \cdot h(x_1, \dots, x_n) 的形式。
f(x1,…,xn∣θ)=θn⋅(∏i=1nxi−2)for θ≤X(1)f(x_1, \dots, x_n | \theta) = \theta^n \cdot \left( \prod_{i=1}^n x_i^{-2} \right) \quad \text{for } \theta \le X_{(1)}
我們需要將條件 θ≤X(1)\theta \le X_{(1)} 包含在 gg 或 hh 中。

令 g(T,θ)=θn⋅I(θ≤T)g(T, \theta) = \theta^n \cdot I(\theta \le T),其中 T=X(1)T = X_{(1)},且 I(⋅)I(\cdot) 是指示函數。
令 h(x1,…,xn)=∏i=1nxi−2h(x_1, \dots, x_n) = \prod_{i=1}^n x_i^{-2}。
這個分解似乎是有效的,但我們需要確保 hh 不依賴於 θ\theta。

檢查分解:
f(x1,…,xn∣θ)=g(X(1),θ)⋅h(x1,…,xn)f(x_1, \dots, x_n | \theta) = g(X_{(1)}, \theta) \cdot h(x_1, \dots, x_n)

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 II.3.b 題5 分

Let X1,X2,…,XnX_1, X_2, \dots, X_n be a random sample from the pdf f(x∣θ)=θx−2,0<θ≤x<∞f(x|\theta) = \theta x^{-2}, 0 < \theta \le x < \infty. Find the maximum likelihood estimator (MLE) of θ\theta.

登入後即可作答並保存紀錄。

這一題的完整詳解

本題要求我們找到參數 θ\theta 的最大概似估計量 (Maximum Likelihood Estimator, MLE)。

步驟一:寫出概似函數 (Likelihood Function)
聯合機率密度函數 (joint pdf) 為:
f(x1,…,xn∣θ)=∏i=1nf(xi∣θ)=∏i=1n(θxi−2)f(x_1, \dots, x_n | \theta) = \prod_{i=1}^n f(x_i|\theta) = \prod_{i=1}^n (\theta x_i^{-2})
f(x1,…,xn∣θ)=θn(∏i=1nxi−2)f(x_1, \dots, x_n | \theta) = \theta^n \left( \prod_{i=1}^n x_i^{-2} \right)
這個函數對於 θ≤xi<∞\theta \le x_i < \infty 對所有 ii 成立。這等價於 θ≤min⁡(x1,…,xn)\theta \le \min(x_1, \dots, x_n)。
令 X(1)=min⁡(x1,…,xn)X_{(1)} = \min(x_1, \dots, x_n)。
則概似函數為:
L(θ)=θn(∏i=1nxi−2)for 0<θ≤X(1)L(\theta) = \theta^n \left( \prod_{i=1}^n x_i^{-2} \right) \quad \text{for } 0 < \theta \le X_{(1)}
且 L(θ)=0L(\theta) = 0 for θ>X(1)\theta > X_{(1)}。

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 II.3.c 題5 分

Let X1,X2,…,XnX_1, X_2, \dots, X_n be a random sample from the pdf f(x∣θ)=θx−2,0<θ≤x<∞f(x|\theta) = \theta x^{-2}, 0 < \theta \le x < \infty. Find the method of moments estimator of θ\theta.

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

方法矩估計法(Method of Moments)是令樣本矩等於母體矩。對單一未知參數 θ\theta,通常使用一階樣本矩:

X‾=1n∑i=1nXi\overline{X}=\frac{1}{n}\sum_{i=1}^n X_i

並令其等於母體期望值:

X‾=Eθ(X)\overline{X}=E_\theta(X)

本題的關鍵在於判斷 Eθ(X)E_\theta(X) 是否存在且為有限值。

解題方法

已知機率密度函數為

f(x∣θ)=θx−2,θ≤x<∞,f(x\mid\theta)=\theta x^{-2},\qquad \theta\le x<\infty,

其中 0<θ≤x0<\theta\le x。

先計算母體平均數:

Eθ(X)=∫θ∞xf(x∣θ) dx=∫θ∞x⋅θx−2 dx.E_\theta(X) =\int_\theta^\infty x f(x\mid\theta)\,dx =\int_\theta^\infty x\cdot \theta x^{-2}\,dx.

整理得

Eθ(X)=θ∫θ∞1x dx.E_\theta(X) =\theta\int_\theta^\infty \frac{1}{x}\,dx.

由於

∫θ∞1x dx=∞,\int_\theta^\infty \frac{1}{x}\,dx=\infty,

因此

Eθ(X)=∞.E_\theta(X)=\infty.
🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 II.4 題5 分

Consider the testing problem H0:θ=θ0H_0: \theta = \theta_0 versus H1:θ≠θ0H_1: \theta \neq \theta_0 based on i.i.d. X1,X2,…,XnX_1, X_2, \dots, X_n from U(0,θ)U(0, \theta). Show that the uniformly most powerful (UMP) test with rejection region X(n)>θ0α1/nX_{(n)} > \theta_0 \alpha^{1/n} or X(n)≤θ0(1−α)1/nX_{(n)} \le \theta_0 (1-\alpha)^{1/n} is a likelihood ratio (LR) test with size α\alpha.

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

本題考查:

  • U(0,θ)U(0,\theta) 母體下最大次序統計量 X(n)X_{(n)} 的分布。
  • 似然比檢定(likelihood ratio test, LR test)。
  • 複合對立假設下的 UMP 檢定。
  • 檢定大小(size)的計算。

令

M=X(n)=max⁡(X1,…,Xn).M=X_{(n)}=\max(X_1,\ldots,X_n).

若 Xi∼i.i.d.U(0,θ)X_i\overset{\text{i.i.d.}}{\sim}U(0,\theta),則

Pθ(M≤m)=Pθ(X1≤m,…,Xn≤m)=(mθ)n,0<m<θ.P_\theta(M\le m) =P_\theta(X_1\le m,\ldots,X_n\le m) =\left(\frac{m}{\theta}\right)^n, \qquad 0<m<\theta.

因此,在 H0:θ=θ0H_0:\theta=\theta_0 下,

Pθ0(M≤m)=(mθ0)n.P_{\theta_0}(M\le m)=\left(\frac{m}{\theta_0}\right)^n.

解題方法:建立似然比

樣本的聯合密度為

L(θ;x)=θ−nI(0<x(1), x(n)≤θ).L(\theta;\mathbf{x}) =\theta^{-n}I(0<x_{(1)},\ x_{(n)}\le \theta).

由於似然函數在 θ≥x(n)\theta\ge x_{(n)} 時為 θ−n\theta^{-n},且隨 θ\theta 遞減,因此最大概似估計量為

θ^=x(n).\widehat{\theta}=x_{(n)}.

考慮廣義似然比

Λ(x)=L(θ0;x)sup⁡θ>0L(θ;x).\Lambda(\mathbf{x}) = \frac{L(\theta_0;\mathbf{x})} {\sup_{\theta>0}L(\theta;\mathbf{x})}.

分兩種情況討論。

情況一:x(n)>θ0x_{(n)}>\theta_0

此時 L(θ0;x)=0L(\theta_0;\mathbf{x})=0,因為樣本中有觀察值大於 θ0\theta_0。所以

Λ(x)=0.\Lambda(\mathbf{x})=0.

因此,所有滿足

X(n)>θ0X_{(n)}>\theta_0

的樣本都應拒絕 H0H_0。

情況二:x(n)≤θ0x_{(n)}\le\theta_0

此時

L(θ0;x)=θ0−n,L(\theta_0;\mathbf{x})=\theta_0^{-n},

而

sup⁡θ>0L(θ;x)=x(n)−n.\sup_{\theta>0}L(\theta;\mathbf{x}) =x_{(n)}^{-n}.

故

Λ(x)=θ0−nx(n)−n=(x(n)θ0)n.\Lambda(\mathbf{x}) = \frac{\theta_0^{-n}}{x_{(n)}^{-n}} = \left(\frac{x_{(n)}}{\theta_0}\right)^n.

依 LR 檢定原理,Λ\Lambda 越小越應拒絕 H0H_0。若將臨界值取為 α\alpha,則在 x(n)≤θ0x_{(n)}\le\theta_0 時,

(x(n)θ0)n≤α,\left(\frac{x_{(n)}}{\theta_0}\right)^n\le\alpha,

等價於

x(n)≤θ0α1/n.x_{(n)}\le\theta_0\alpha^{1/n}.

因此,正確的 LR 拒絕域為

R={X(n)>θ0}∪{X(n)≤θ0α1/n}.\mathcal{R} = \left\{X_{(n)}>\theta_0\right\} \cup \left\{X_{(n)}\le\theta_0\alpha^{1/n}\right\}.

檢查檢定大小

在 H0:θ=θ0H_0:\theta=\theta_0 下,

Pθ0(X(n)>θ0)=0,P_{\theta_0}\left(X_{(n)}>\theta_0\right)=0,
🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

其他考古題