115 年 國立政治大學統計學系碩士班《數理統計學》

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

第 1 題10 分

Let X be a random variable and p∈(0,∞)p \in (0, \infty). Assume the p-th moment of X exists. Prove that
E[∣X∣p]=∫0∞ptp−1Pr⁡(∣X∣>t)dtE[|X|^p] = \int_0^\infty pt^{p-1} \Pr(|X| > t) dt

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

這一題的完整詳解

本題旨在證明隨機變數的 pp 次絕對動差(p-th absolute moment)與其累積分布函數(CDF)的關係。核心概念是利用積分的變數替換和機率測度的性質。

證明過程如下:
令 Y=∣X∣pY = |X|^p。由於 p>0p > 0,則 Y≥0Y \ge 0。
我們知道對於一個非負隨機變數 YY,其期望值可以表示為:
E[Y]=∫0∞Pr⁡(Y>y)dyE[Y] = \int_0^\infty \Pr(Y > y) dy
將 Y=∣X∣pY = |X|^p 代入上式,得到:
E[∣X∣p]=∫0∞Pr⁡(∣X∣p>y)dyE[|X|^p] = \int_0^\infty \Pr(|X|^p > y) dy
令 y=tpy = t^p。

🔒

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

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

免費註冊

第 2 題10 分

Let μ(z)=E(Y∣Z=z)\mu(z) = E(Y | Z = z). Show that
Var(μ(Z))Var(Y)=Corr2(Y,μ(Z))=sup⁡g∈L2(Z)Corr2(Y,g(Z))\frac{\text{Var}(\mu(Z))}{\text{Var}(Y)} = \text{Corr}^2(Y, \mu(Z)) = \sup_{g \in L^2(Z)} \text{Corr}^2(Y, g(Z))
where the supremum is taken over all square-integrable functions of Z.

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

這一題的完整詳解

本題是關於條件期望的性質,特別是全期望法則(Law of Total Expectation)和條件變異數(Conditional Variance)與總變異數(Total Variance)的關係,以及其與相關係數的連結。

核心觀念:

  1. 全期望法則: E[Y]=E[E[Y∣Z]]E[Y] = E[E[Y|Z]]
  2. 全變異數法則: Var(Y)=E[Var(Y∣Z)]+Var(E[Y∣Z])\text{Var}(Y) = E[\text{Var}(Y|Z)] + \text{Var}(E[Y|Z])
  3. 相關係數定義: Corr(U,V)=Cov(U,V)Var(U)Var(V)\text{Corr}(U, V) = \frac{\text{Cov}(U, V)}{\sqrt{\text{Var}(U)\text{Var}(V)}}
  4. 最佳線性預測(BLP)與條件期望的關係: 如果 YY 和 ZZ 之間存在線性關係,則 E[Y∣Z]E[Y|Z] 是 YY 對 ZZ 的最佳預測。更一般地,對於給定的 ZZ, E[Y∣Z]E[Y|Z] 是 YY 的最佳預測。

證明步驟:

第一部分:證明 Var(μ(Z))=Corr2(Y,μ(Z))Var(Y)\text{Var}(\mu(Z)) = \text{Corr}^2(Y, \mu(Z)) \text{Var}(Y)

由全變異數法則,我們有:
Var(Y)=E[Var(Y∣Z)]+Var(E[Y∣Z])\text{Var}(Y) = E[\text{Var}(Y|Z)] + \text{Var}(E[Y|Z])
令 μ(Z)=E[Y∣Z]\mu(Z) = E[Y|Z]。則上式變為:
Var(Y)=E[Var(Y∣Z)]+Var(μ(Z))\text{Var}(Y) = E[\text{Var}(Y|Z)] + \text{Var}(\mu(Z))
由於 Var(Y∣Z)≥0\text{Var}(Y|Z) \ge 0,我們有 E[Var(Y∣Z)]≥0E[\text{Var}(Y|Z)] \ge 0。
因此,Var(Y)≥Var(μ(Z))\text{Var}(Y) \ge \text{Var}(\mu(Z))。

現在我們需要計算 Cov(Y,μ(Z))\text{Cov}(Y, \mu(Z))。
Cov(Y,μ(Z))=E[Yμ(Z)]−E[Y]E[μ(Z)]\text{Cov}(Y, \mu(Z)) = E[Y \mu(Z)] - E[Y] E[\mu(Z)]
利用全期望法則:
E[Y]=E[E[Y∣Z]]=E[μ(Z)]E[Y] = E[E[Y|Z]] = E[\mu(Z)]。
所以,
Cov(Y,μ(Z))=E[Yμ(Z)]−(E[μ(Z)])2\text{Cov}(Y, \mu(Z)) = E[Y \mu(Z)] - (E[\mu(Z)])^2
另外,
E[Yμ(Z)]=E[E[Yμ(Z)∣Z]]E[Y \mu(Z)] = E[E[Y \mu(Z) | Z]]
由於 μ(Z)\mu(Z) 是 ZZ 的函數,可以提到條件期望外面:
E[Yμ(Z)]=E[μ(Z)E[Y∣Z]]=E[μ(Z)μ(Z)]=E[μ(Z)2]E[Y \mu(Z)] = E[\mu(Z) E[Y | Z]] = E[\mu(Z) \mu(Z)] = E[\mu(Z)^2]
所以,
Cov(Y,μ(Z))=E[μ(Z)2]−(E[μ(Z)])2=Var(μ(Z))\text{Cov}(Y, \mu(Z)) = E[\mu(Z)^2] - (E[\mu(Z)])^2 = \text{Var}(\mu(Z))

現在計算相關係數的平方:
Corr2(Y,μ(Z))=Cov(Y,μ(Z))2Var(Y)Var(μ(Z))=(Var(μ(Z)))2Var(Y)Var(μ(Z))\text{Corr}^2(Y, \mu(Z)) = \frac{\text{Cov}(Y, \mu(Z))^2}{\text{Var}(Y) \text{Var}(\mu(Z))} = \frac{(\text{Var}(\mu(Z)))^2}{\text{Var}(Y) \text{Var}(\mu(Z))}
Corr2(Y,μ(Z))=Var(μ(Z))Var(Y)\text{Corr}^2(Y, \mu(Z)) = \frac{\text{Var}(\mu(Z))}{\text{Var}(Y)}
從這個等式,我們可以直接得到:
Var(μ(Z))Var(Y)=Corr2(Y,μ(Z))\frac{\text{Var}(\mu(Z))}{\text{Var}(Y)} = \text{Corr}^2(Y, \mu(Z))
這個等式成立的前提是 Var(Y)>0\text{Var}(Y) > 0 且 Var(μ(Z))>0\text{Var}(\mu(Z)) > 0。如果 Var(Y)=0\text{Var}(Y) = 0,則 YY 是常數,所有期望和變異數都為零,公式仍成立(0/0 通常定義為 1,但若 Y 是常數,Y 和任何函數 g(Z) 的相關係數也是 0,所以 Corr^2 也是 0,公式 0/0 = 0 似乎有問題,但通常會假設 Y 不是常數)。如果 Var(μ(Z))=0\text{Var}(\mu(Z)) = 0,則 μ(Z)\mu(Z) 是常數,表示 E[Y∣Z]E[Y|Z] 是常數,這意味著 YY 的值不依賴於 ZZ。此時 Var(Y)=E[Var(Y∣Z)]\text{Var}(Y) = E[\text{Var}(Y|Z)],且 Cov(Y,μ(Z))=0\text{Cov}(Y, \mu(Z)) = 0,Corr2=0\text{Corr}^2=0。公式 0Var(Y)=0\frac{0}{\text{Var}(Y)}=0 成立。

第二部分:證明 Corr2(Y,μ(Z))=sup⁡g∈L2(Z)Corr2(Y,g(Z))\text{Corr}^2(Y, \mu(Z)) = \sup_{g \in L^2(Z)} \text{Corr}^2(Y, g(Z))

🔒

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

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

免費註冊

第 3 題

Let {X~1,…,X~N}\{\tilde{X}_1, \dots, \tilde{X}_N\} and {(X1,Y1),…,(Xn,Yn)}\{(X_1, Y_1), \dots, (X_n, Y_n)\} be independent sets of random samples from the same population.
(a) (10%) Given a prediction function ff, show the estimator
θ^=1N∑i=1Nf(X~i)+1n∑j=1n(Yj−f(Xj))\hat{\theta} = \frac{1}{N} \sum_{i=1}^N f(\tilde{X}_i) + \frac{1}{n} \sum_{j=1}^n (Y_j - f(X_j))
is unbiased for θ=E(Y)\theta = E(Y).

(b) (10%) Find the condition on the correlation between YY and f(X)f(X) to make θ^\hat{\theta} have smaller MSE than Yˉn=1n∑j=1nYj\bar{Y}_n = \frac{1}{n} \sum_{j=1}^n Y_j in terms of the standard deviations of YY and f(X)f(X).

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

這一題的完整詳解

核心觀念

本題考查:

  1. 期望值的線性性與無偏估計量定義。
  2. 獨立隨機樣本所造成的變異數可加性。
  3. 均方誤差公式
    MSE⁡(θ^)=Var⁡(θ^)+Bias⁡(θ^)2.\operatorname{MSE}(\hat{\theta})=\operatorname{Var}(\hat{\theta})+\operatorname{Bias}(\hat{\theta})^2.
  4. 共變異數與相關係數的關係
    Cov⁡(Y,f(X))=ρY,f(X)σYσf(X).\operatorname{Cov}(Y,f(X))=\rho_{Y,f(X)}\sigma_Y\sigma_{f(X)}.

令

μf=E[f(X)],σf2=Var⁡(f(X)),σY2=Var⁡(Y).\mu_f=E[f(X)],\qquad \sigma_f^2=\operatorname{Var}(f(X)),\qquad \sigma_Y^2=\operatorname{Var}(Y).

由於兩組樣本來自同一母體,X~\tilde X 與 XX 具有相同分布,因此

E[f(X~)]=E[f(X)]=μf.E[f(\tilde X)]=E[f(X)]=\mu_f.

(a)證明 θ^\hat{\theta} 為 θ=E(Y)\theta=E(Y) 的無偏估計量

題目中的估計量為

θ^=1N∑i=1Nf(X~i)+1n∑j=1n(Yj−f(Xj)).\hat{\theta} = \frac{1}{N}\sum_{i=1}^N f(\tilde X_i) + \frac{1}{n}\sum_{j=1}^n\bigl(Y_j-f(X_j)\bigr).

利用期望值的線性性,

E(θ^)=E[1N∑i=1Nf(X~i)]+E[1n∑j=1n(Yj−f(Xj))].E(\hat{\theta}) = E\left[\frac{1}{N}\sum_{i=1}^N f(\tilde X_i)\right] + E\left[\frac{1}{n}\sum_{j=1}^n\bigl(Y_j-f(X_j)\bigr)\right].

因為 X~i\tilde X_i 為同分布樣本,

E[1N∑i=1Nf(X~i)]=1N∑i=1NE[f(X~i)]=μf.E\left[\frac{1}{N}\sum_{i=1}^N f(\tilde X_i)\right] = \frac{1}{N}\sum_{i=1}^N E[f(\tilde X_i)] = \mu_f.

同理,

E[1n∑j=1n(Yj−f(Xj))]=E[Y−f(X)]=E(Y)−E[f(X)]=θ−μf.E\left[\frac{1}{n}\sum_{j=1}^n\bigl(Y_j-f(X_j)\bigr)\right] = E[Y-f(X)] = E(Y)-E[f(X)] = \theta-\mu_f.

因此

E(θ^)=μf+(θ−μf)=θ.E(\hat{\theta}) = \mu_f+(\theta-\mu_f) = \theta.

所以

E(θ^)=θ\boxed{E(\hat{\theta})=\theta}

故 θ^\hat{\theta} 是 θ=E(Y)\theta=E(Y) 的無偏估計量。


(b)使 θ^\hat{\theta} 的 MSE 小於 Yˉn\bar Y_n 的相關係數條件

由(a)知 θ^\hat{\theta} 無偏,因此

MSE⁡(θ^)=Var⁡(θ^).\operatorname{MSE}(\hat{\theta}) = \operatorname{Var}(\hat{\theta}).

估計量的兩部分分別來自兩組獨立樣本,因此兩部分獨立,變異數可以直接相加:

Var⁡(θ^)=Var⁡(1N∑i=1Nf(X~i))+Var⁡(1n∑j=1n(Yj−f(Xj))).\operatorname{Var}(\hat{\theta}) = \operatorname{Var}\left(\frac{1}{N}\sum_{i=1}^N f(\tilde X_i)\right) + \operatorname{Var}\left(\frac{1}{n}\sum_{j=1}^n(Y_j-f(X_j))\right).

第一部分為

Var⁡(1N∑i=1Nf(X~i))=σf2N.\operatorname{Var}\left(\frac{1}{N}\sum_{i=1}^N f(\tilde X_i)\right) = \frac{\sigma_f^2}{N}.

第二部分為

🔒

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

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

免費註冊

第 4 題

Let X1,…,XnX_1, \dots, X_n be i.i.d. N(μ,σ2)N(\mu, \sigma^2) random variables with unknown μ∈R\mu \in \mathbb{R} and known σ2>0\sigma^2 > 0.
(a) (15%) Find the UMVUE of θ=etμ\theta = e^{t\mu} with a fixed t≠0t \neq 0.
(b) (15%) Determine whether the variance of the UMVUE in (a) attains the CRLB.

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

這一題的完整詳解

核心觀念

本題考查:

  1. 充分統計量與完備統計量。
  2. Lehmann–Scheffé 定理:若某估計量是完備充分統計量的函數,且對目標參數不偏,則必為 UMVUE。
  3. 常態分配指數函數之期望。
  4. Cramér–Rao lower bound(CRLB)及其與實際變異數的比較。

令

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

因為

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

解題方法

(a)求 θ=etμ\theta=e^{t\mu} 的 UMVUE

第一步:找完備充分統計量

樣本聯合密度為

f(x;μ)=(2πσ2)−n/2exp⁡{−12σ2∑i=1n(xi−μ)2}.f(\boldsymbol{x};\mu) =(2\pi\sigma^2)^{-n/2} \exp\left\{-\frac{1}{2\sigma^2}\sum_{i=1}^n(x_i-\mu)^2\right\}.

整理與 μ\mu 有關的部分:

f(x;μ)=h(x)exp⁡{μσ2∑i=1nxi−nμ22σ2}.f(\boldsymbol{x};\mu) =h(\boldsymbol{x}) \exp\left\{ \frac{\mu}{\sigma^2}\sum_{i=1}^n x_i -\frac{n\mu^2}{2\sigma^2} \right\}.

由 Neyman–Fisher factorization theorem,

T=∑i=1nXiT=\sum_{i=1}^n X_i

是 μ\mu 的充分統計量,因此 Xˉ\bar X 也是充分統計量。

這是自然參數空間為整個 R\mathbb{R} 的一參數滿秩指數族,因此 TT 為完備統計量;相應地,Xˉ\bar X 也是完備充分統計量。


第二步:建構不偏估計量

由

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

若 Y∼N(m,v)Y\sim N(m,v),則

E(eaY)=exp⁡(am+a2v2).E(e^{aY})=\exp\left(am+\frac{a^2v}{2}\right).

取 a=ta=t,得到

Eμ(etXˉ)=exp⁡(tμ+t2σ22n).E_\mu(e^{t\bar X}) = \exp\left( t\mu+\frac{t^2\sigma^2}{2n} \right).

因此修正常數後,

θ^=exp⁡(−t2σ22n)etXˉ=exp⁡(tXˉ−t2σ22n).\widehat\theta = \exp\left(-\frac{t^2\sigma^2}{2n}\right)e^{t\bar X} = \exp\left( t\bar X-\frac{t^2\sigma^2}{2n} \right).

其期望為

Eμ(θ^)=exp⁡(−t2σ22n)Eμ(etXˉ)=etμ=θ.E_\mu(\widehat\theta) = \exp\left(-\frac{t^2\sigma^2}{2n}\right) E_\mu(e^{t\bar X}) = e^{t\mu} = \theta.

故 θ^\widehat\theta 是 θ\theta 的不偏估計量,且它是完備充分統計量 Xˉ\bar X 的函數。依 Lehmann–Scheffé 定理,θ^\widehat\theta 即為 UMVUE。


(b)變異數是否達到 CRLB?

第一步:計算 UMVUE 的變異數

先計算二次矩:

Eμ(θ^2)=Eμ[exp⁡(2tXˉ−t2σ2n)].E_\mu(\widehat\theta^2) = E_\mu\left[ \exp\left( 2t\bar X-\frac{t^2\sigma^2}{n} \right) \right].

由常態分配指數函數期望公式,

Eμ(e2tXˉ)=exp⁡(2tμ+(2t)22σ2n)=exp⁡(2tμ+2t2σ2n).E_\mu(e^{2t\bar X}) = \exp\left( 2t\mu+\frac{(2t)^2}{2}\frac{\sigma^2}{n} \right) = \exp\left( 2t\mu+\frac{2t^2\sigma^2}{n} \right).

所以

Eμ(θ^2)=exp⁡(2tμ+t2σ2n).E_\mu(\widehat\theta^2) = \exp\left( 2t\mu+\frac{t^2\sigma^2}{n} \right).

又因為

Eμ(θ^)=etμ,E_\mu(\widehat\theta)=e^{t\mu},

故

Var⁡μ(θ^)=Eμ(θ^2)−[Eμ(θ^)]2.\operatorname{Var}_\mu(\widehat\theta) = E_\mu(\widehat\theta^2)-\left[E_\mu(\widehat\theta)\right]^2.

因此

🔒

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

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

免費註冊

第 5 題

Let X1X_1 and X2X_2 be independent Exponential random variables, with density
f(x;θi)=1θie−x/θi,x>0,i=1,2f(x; \theta_i) = \frac{1}{\theta_i} e^{-x/\theta_i}, \quad x > 0, \quad i = 1, 2
Define γ=θ1/θ2\gamma = \theta_1 / \theta_2. Find a LR test of size α\alpha for testing
(a) (15%) H0:γ=1H_0: \gamma = 1 versus H1:γ≠1H_1: \gamma \neq 1
(b) (15%) H0:γ≤1H_0: \gamma \le 1 versus H1:γ>1H_1: \gamma > 1

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

這一題的完整詳解

核心觀念

本題考查:

  1. 以概似比統計量建立 LR test。
  2. 指數分布尺度參數的最大概似估計。
  3. 利用比例統計量
    R=X1X1+X2R=\frac{X_1}{X_1+X_2}
    消除共同尺度參數。
  4. 在虛無假設為複合假設時,確認檢定的最大型一錯誤率發生在哪個邊界點。

聯合概似函數為

L(θ1,θ2)=1θ1θ2exp⁡(−x1θ1−x2θ2),x1,x2>0.L(\theta_1,\theta_2) = \frac{1}{\theta_1\theta_2} \exp\left(-\frac{x_1}{\theta_1}-\frac{x_2}{\theta_2}\right), \qquad x_1,x_2>0.

令

γ=θ1θ2,θ1=γθ2.\gamma=\frac{\theta_1}{\theta_2}, \qquad \theta_1=\gamma\theta_2.

比例統計量的分布

令

R=X1X1+X2,S=X1+X2.R=\frac{X_1}{X_1+X_2}, \qquad S=X_1+X_2.

則

X1=RS,X2=(1−R)S,X_1=RS,\qquad X_2=(1-R)S,

且轉換的 Jacobian 為 SS。

將 θ1=γθ2\theta_1=\gamma\theta_2 代入聯合密度,可得

fR,S(r,s)=sγθ22exp⁡[−sθ2(rγ+1−r)].f_{R,S}(r,s) = \frac{s}{\gamma\theta_2^2} \exp\left[ -\frac{s}{\theta_2} \left(\frac{r}{\gamma}+1-r\right) \right].

對 ss 積分後,得到

fR(r;γ)=γ[r+γ(1−r)]2,0<r<1.f_R(r;\gamma) = \frac{\gamma} {\left[r+\gamma(1-r)\right]^2}, \qquad 0<r<1.

特別地,在 γ=1\gamma=1 時,

fR(r;1)=1,0<r<1,f_R(r;1)=1,\qquad 0<r<1,

所以

R∣(γ=1)∼Uniform⁡(0,1).R\mid(\gamma=1)\sim \operatorname{Uniform}(0,1).

其累積分布函數為

FR(r;γ)=rr+γ(1−r).F_R(r;\gamma) = \frac{r}{r+\gamma(1-r)}.

因此

Pγ(R>c)=1−FR(c;γ)=γ(1−c)c+γ(1−c).P_\gamma(R>c) = 1-F_R(c;\gamma) = \frac{\gamma(1-c)} {c+\gamma(1-c)}.

(a) H0:γ=1H_0:\gamma=1 對 H1:γ≠1H_1:\gamma\neq 1

解題方法

此處使用概似比檢定:

Λ(x1,x2)=sup⁡θ1=θ2=θL(θ,θ)sup⁡θ1>0,θ2>0L(θ1,θ2).\Lambda(x_1,x_2) = \frac{\displaystyle\sup_{\theta_1=\theta_2=\theta}L(\theta,\theta)} {\displaystyle\sup_{\theta_1>0,\theta_2>0}L(\theta_1,\theta_2)}.

分母:不受限制下的最大概似

對 θ1,θ2\theta_1,\theta_2 分別最大化,得到

θ^1=X1,θ^2=X2.\widehat{\theta}_1=X_1,\qquad \widehat{\theta}_2=X_2.

所以

sup⁡L=1X1X2e−2.\sup L = \frac{1}{X_1X_2}e^{-2}.

分子:在 H0H_0 下的最大概似

在 H0:θ1=θ2=θH_0:\theta_1=\theta_2=\theta 下,

L0(θ)=1θ2exp⁡(−x1+x2θ).L_0(\theta) = \frac{1}{\theta^2} \exp\left(-\frac{x_1+x_2}{\theta}\right).

其最大概似估計為

θ^0=x1+x22.\widehat{\theta}_0=\frac{x_1+x_2}{2}.

因此

sup⁡H0L=4(x1+x2)2e−2.\sup_{H_0}L = \frac{4}{(x_1+x_2)^2}e^{-2}.

故概似比統計量為

Λ=4x1x2(x1+x2)2.\Lambda = \frac{4x_1x_2}{(x_1+x_2)^2}.

以

R=X1X1+X2R=\frac{X_1}{X_1+X_2}

表示時,

Λ=4R(1−R).\Lambda=4R(1-R).

當 RR 接近 00 或 11 時,Λ\Lambda 很小,表示兩個尺度參數差異很大,因此應拒絕 H0H_0。

臨界區

在 H0H_0 下,R∼Uniform⁡(0,1)R\sim\operatorname{Uniform}(0,1)。由於檢定為雙尾檢定,令

P0(R<c)+P0(R>1−c)=α.P_0(R<c)+P_0(R>1-c)=\alpha.

因為 RR 在 (0,1)(0,1) 上均勻分布,

2c=α,c=α2.2c=\alpha, \qquad c=\frac{\alpha}{2}.

所以拒絕域為

R<α2或R>1−α2.R<\frac{\alpha}{2} \quad\text{或}\quad R>1-\frac{\alpha}{2}.

換回 X1,X2X_1,X_2:

🔒

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

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

免費註冊

其他考古題