115 年 國立臺灣大學經濟學系碩士班《統計學(A)》

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

第 1 題5 分

Comment on the following statement concerning the Central Limit Theorem.

"For independent and identically distributed random variables with finite mean μ\mu and variance σ2\sigma^2, the sample average Xˉ\bar{X} converges in distribution to a normal distribution with mean μ\mu and variance σ2n\frac{\sigma^2}{n} as the sample size increases."

You need to provide an explanation for why the statement is TRUE or FALSE.

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

這一題的完整詳解

這個題目在考驗對中央極限定理 (Central Limit Theorem, CLT) 的理解。CLT 的核心在於說明,當樣本數 nn 夠大時,來自任何分布(只要其一階與二階動差存在)的獨立同分布隨機變數的樣本平均值,其分布會趨近於一個常態分布。

CLT 的標準敘述是:
設 X1,X2,…,XnX_1, X_2, \dots, X_n 為獨立且同分布(i.i.d.)的隨機變數,具有有限的期望值 E(Xi)=μE(X_i) = \mu 和有限的變異數 Var(Xi)=σ2>0Var(X_i) = \sigma^2 > 0。則樣本平均值 Xˉ=1n∑i=1nXi\bar{X} = \frac{1}{n}\sum_{i=1}^n X_i 的標準化形式,當 n→∞n \to \infty 時,會趨近於標準常態分布 N(0,1)N(0, 1)。
即:
lim⁡n→∞P(Xˉ−μσ/n≤z)=Φ(z)\lim_{n \to \infty} P\left(\frac{\bar{X} - \mu}{\sigma/\sqrt{n}} \le z\right) = \Phi(z)
其中 Φ(z)\Phi(z) 是標準常態分布的累積分布函數。

從上述標準敘述,我們可以推導出 Xˉ\bar{X} 的分布:
由於 Xˉ−μσ/n→dN(0,1)\frac{\bar{X} - \mu}{\sigma/\sqrt{n}} \xrightarrow{d} N(0, 1),
這意味著 Xˉ−μ→dN(0,(σ/n)2)=N(0,σ2/n)\bar{X} - \mu \xrightarrow{d} N(0, (\sigma/\sqrt{n})^2) = N(0, \sigma^2/n)。
進一步,這表示 Xˉ→dN(μ,σ2/n)\bar{X} \xrightarrow{d} N(\mu, \sigma^2/n)。

🔒

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

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

免費註冊

第 2 題15 分

Let (X1,X2,X3)∼(X_1, X_2, X_3) \sim i.i.d. N(μ,σ2)N(\mu, \sigma^2).

(a) Define Y=X1X2Y = \frac{X_1}{X_2}. Find the distribution of YY.

(b) Define W=X1−X2X1+X2W = \frac{X_1 - X_2}{X_1 + X_2}. Find the distribution of WW.

(c) Consider the following two statistics: U=(X1+X2+X3)2U = (X_1 + X_2 + X_3)^2 and VV is the range of the random sample (X1,X2,X3)(X_1, X_2, X_3). Show that UU and VV are independent.

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

這一題的完整詳解

核心觀念

本題考查三個重點:

  1. 兩個獨立常態變數之比的分布。當分子、分母皆為零均值時,比例服從標準 Cauchy 分布;均值不為零時,通常不是標準 Cauchy 分布。
  2. 常態變數的線性轉換,以及獨立常態變數之比的密度推導。
  3. 常態樣本中,樣本平均數與離均差向量獨立;而樣本全距只依賴離均差,因此可證明其與樣本平均數平方獨立。

以下假設 σ2>0\sigma^2>0,並以 Φ\Phi 表示標準常態分布的累積分布函數。


(a) Y=X1X2Y=\dfrac{X_1}{X_2} 的分布

解題方法

使用變數轉換

Y=X1X2,T=X2.Y=\frac{X_1}{X_2},\qquad T=X_2.

因此

X1=YT,X2=T.X_1=YT,\qquad X_2=T.

Jacobian 的絕對值為

∣∂(x1,x2)∂(y,t)∣=∣t∣.\left|\frac{\partial(x_1,x_2)}{\partial(y,t)}\right|=|t|.

由於 X1,X2X_1,X_2 獨立且皆服從 N(μ,σ2)N(\mu,\sigma^2),

fX1,X2(x1,x2)=12πσ2exp⁡{−(x1−μ)2+(x2−μ)22σ2}.f_{X_1,X_2}(x_1,x_2) = \frac{1}{2\pi\sigma^2} \exp\left\{ -\frac{(x_1-\mu)^2+(x_2-\mu)^2}{2\sigma^2} \right\}.

所以

fY(y)=∫−∞∞∣t∣fX1,X2(yt,t) dt.f_Y(y) = \int_{-\infty}^{\infty} |t|f_{X_1,X_2}(yt,t)\,dt.

代入並整理指數項:

fY(y)=12πσ2∫−∞∞∣t∣exp⁡{−(1+y2)t2−2μ(1+y)t+2μ22σ2}dt.f_Y(y) = \frac{1}{2\pi\sigma^2} \int_{-\infty}^{\infty} |t| \exp\left\{ -\frac{(1+y^2)t^2-2\mu(1+y)t+2\mu^2}{2\sigma^2} \right\} dt.

完成平方後,可得 YY 的機率密度函數:

fY(y)=exp⁡{−μ2(y−1)22σ2(1+y2)}π(1+y2)f_Y(y) = \frac{ \exp\left\{ -\dfrac{\mu^2(y-1)^2}{2\sigma^2(1+y^2)} \right\} }{\pi(1+y^2)} +μ(y+1)exp⁡{−μ2(y−1)22σ2(1+y2)}σ2π(1+y2)3/2[2Φ(μ(y+1)σ1+y2)−1],−∞<y<∞.\quad+ \frac{ \mu(y+1) \exp\left\{ -\dfrac{\mu^2(y-1)^2}{2\sigma^2(1+y^2)} \right\} }{ \sigma\sqrt{2\pi}(1+y^2)^{3/2} } \left[ 2\Phi\left( \frac{\mu(y+1)}{\sigma\sqrt{1+y^2}} \right)-1 \right], \qquad -\infty<y<\infty.

因此,YY 服從「兩個同變異數、同均值常態變數之比」的 ratio-normal 分布。

特別地,若 μ=0\mu=0,上式化為

fY(y)=1π(1+y2),f_Y(y)=\frac{1}{\pi(1+y^2)},

故

Y∼Cauchy⁡(0,1).Y\sim \operatorname{Cauchy}(0,1).

(b) W=X1−X2X1+X2W=\dfrac{X_1-X_2}{X_1+X_2} 的分布

解題方法

令

A=X1−X2,B=X1+X2.A=X_1-X_2,\qquad B=X_1+X_2.

由獨立常態變數的線性組合仍為常態分布,

A∼N(0,2σ2),A\sim N(0,2\sigma^2), B∼N(2μ,2σ2).B\sim N(2\mu,2\sigma^2).

再計算共變異數:

Cov⁡(A,B)=Cov⁡(X1−X2,X1+X2)=σ2−σ2=0.\operatorname{Cov}(A,B) = \operatorname{Cov}(X_1-X_2,X_1+X_2) = \sigma^2-\sigma^2=0.

因為 (A,B)(A,B) 為聯合常態,且共變異數為零,所以 AA 與 BB 獨立。

現在

W=AB.W=\frac{A}{B}.

令 T=BT=B,則

A=WT.A=WT.

Jacobian 的絕對值為 ∣t∣|t|,故

fW(w)=∫−∞∞∣t∣fA(wt)fB(t) dt.f_W(w) = \int_{-\infty}^{\infty} |t|f_A(wt)f_B(t)\,dt.

代入 A∼N(0,2σ2)A\sim N(0,2\sigma^2) 與 B∼N(2μ,2σ2)B\sim N(2\mu,2\sigma^2) 的密度,整理後得到

fW(w)=exp⁡{−μ2w2σ2(1+w2)}π(1+w2)f_W(w) = \frac{ \exp\left\{ -\dfrac{\mu^2w^2}{\sigma^2(1+w^2)} \right\} }{\pi(1+w^2)} +μexp⁡{−μ2w2σ2(1+w2)}σπ(1+w2)3/2[2Φ(2μσ1+w2)−1],−∞<w<∞.\quad+ \frac{ \mu \exp\left\{ -\dfrac{\mu^2w^2}{\sigma^2(1+w^2)} \right\} }{ \sigma\sqrt{\pi}(1+w^2)^{3/2} } \left[ 2\Phi\left( \frac{\sqrt{2}\mu}{\sigma\sqrt{1+w^2}} \right)-1 \right], \qquad -\infty<w<\infty.
🔒

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

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

免費註冊

第 3 題15 分

Let (X1,…,Xn)(X_1, \dots, X_n) be a random sample from U[0,θ]U[0, \theta], θ>0\theta > 0.

(a) Find the method of moments estimator for θ\theta, denoted by θ^\hat{\theta}. Is θ^\hat{\theta} a consistent estimator for θ\theta?

(b) Use θ^\hat{\theta} to propose a test statistic to test H0:θ=9H_0: \theta = 9 vs. H1:θ≠9H_1: \theta \neq 9 at significance level α=0.05\alpha = 0.05. Specify the rejection region.

(c) Given the realization of the random sample (X1,X2,X3)=(6.30,2.28,7.46)(X_1, X_2, X_3) = (6.30, 2.28, 7.46), construct a 95% confidence interval for θ\theta, and determine whether you can reject the null hypothesis at significance level α=0.05\alpha = 0.05 in part (b).

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

這一題的完整詳解

此題考驗對均值矩估計法 (Method of Moments, MOM)、一致性估計量、假設檢定與信賴區間的理解與應用,特別是針對均勻分布。

核心觀念:

  • 均值矩估計法:令樣本動差等於母體動差,解出待估參數。
  • 一致性:當樣本數趨近無限大時,估計量會趨近真實參數。
  • 假設檢定:根據檢定統計量的值,判斷是否拒絕虛無假設。
  • 信賴區間:給定信賴水準,參數落入該區間的機率。
  • 均勻分布 U[0,θ]U[0, \theta] 的性質:期望值、變異數、CDF、PDF。

(a) 尋找 θ\theta 的矩估計量 θ^\hat{\theta},並判斷其一致性。

均勻分布 U[0,θ]U[0, \theta] 的機率密度函數 (PDF) 為:
f(x;θ)={1θ0≤x≤θ0otherwisef(x; \theta) = \begin{cases} \frac{1}{\theta} & 0 \le x \le \theta \\ 0 & \text{otherwise} \end{cases}

步驟 1:計算母體的一階動差 (期望值)。
E(X)=∫−∞∞xf(x;θ)dx=∫0θx⋅1θdx=1θ[x22]0θ=1θ(θ22−0)=θ2E(X) = \int_{-\infty}^{\infty} x f(x; \theta) dx = \int_0^\theta x \cdot \frac{1}{\theta} dx = \frac{1}{\theta} \left[ \frac{x^2}{2} \right]_0^\theta = \frac{1}{\theta} \left( \frac{\theta^2}{2} - 0 \right) = \frac{\theta}{2}。

步驟 2:計算樣本的一階動差。
樣本的一階動差是樣本平均值:
Xˉ=1n∑i=1nXi\bar{X} = \frac{1}{n} \sum_{i=1}^n X_i。

步驟 3:令樣本一階動差等於母體一階動差,解出 θ\theta。
Xˉ=E(X)\bar{X} = E(X)
Xˉ=θ2\bar{X} = \frac{\theta}{2}
解出 θ\theta:
θ=2Xˉ\theta = 2\bar{X}

因此,θ\theta 的矩估計量為 θ^MOM=2Xˉ\hat{\theta}_{MOM} = 2\bar{X}。

步驟 4:判斷 θ^=2Xˉ\hat{\theta} = 2\bar{X} 的一致性。
一個估計量 θ^\hat{\theta} 對於參數 θ\theta 是一致的,如果 lim⁡n→∞P(∣θ^−θ∣<ϵ)=1\lim_{n \to \infty} P(|\hat{\theta} - \theta| < \epsilon) = 1 對於任意 ϵ>0\epsilon > 0。
這等價於 lim⁡n→∞E(θ^)=θ\lim_{n \to \infty} E(\hat{\theta}) = \theta 且 lim⁡n→∞Var(θ^)=0\lim_{n \to \infty} Var(\hat{\theta}) = 0。

  • 檢查期望值:
    E(θ^)=E(2Xˉ)=2E(Xˉ)E(\hat{\theta}) = E(2\bar{X}) = 2 E(\bar{X})
    由於 Xˉ\bar{X} 是來自 U[0,θ]U[0, \theta] 的樣本平均值,其期望值為母體的期望值:E(Xˉ)=E(X)=θ2E(\bar{X}) = E(X) = \frac{\theta}{2}。
    所以,E(θ^)=2⋅θ2=θE(\hat{\theta}) = 2 \cdot \frac{\theta}{2} = \theta。
    lim⁡n→∞E(θ^)=θ\lim_{n \to \infty} E(\hat{\theta}) = \theta,所以 θ^\hat{\theta} 是漸近無偏的 (asymptotically unbiased)。

  • 檢查變異數:
    Var(θ^)=Var(2Xˉ)=22Var(Xˉ)=4Var(Xˉ)Var(\hat{\theta}) = Var(2\bar{X}) = 2^2 Var(\bar{X}) = 4 Var(\bar{X})。
    樣本平均值的變異數為母體變異數除以樣本數:Var(Xˉ)=Var(X)nVar(\bar{X}) = \frac{Var(X)}{n}。
    均勻分布 U[0,θ]U[0, \theta] 的變異數為:
    Var(X)=E(X2)−(E(X))2Var(X) = E(X^2) - (E(X))^2
    E(X2)=∫0θx2⋅1θdx=1θ[x33]0θ=1θθ33=θ23E(X^2) = \int_0^\theta x^2 \cdot \frac{1}{\theta} dx = \frac{1}{\theta} \left[ \frac{x^3}{3} \right]_0^\theta = \frac{1}{\theta} \frac{\theta^3}{3} = \frac{\theta^2}{3}。
    Var(X)=θ23−(θ2)2=θ23−θ24=4θ2−3θ212=θ212Var(X) = \frac{\theta^2}{3} - \left(\frac{\theta}{2}\right)^2 = \frac{\theta^2}{3} - \frac{\theta^2}{4} = \frac{4\theta^2 - 3\theta^2}{12} = \frac{\theta^2}{12}。
    所以,Var(Xˉ)=θ2/12n=θ212nVar(\bar{X}) = \frac{\theta^2/12}{n} = \frac{\theta^2}{12n}。
    Var(θ^)=4⋅θ212n=θ23nVar(\hat{\theta}) = 4 \cdot \frac{\theta^2}{12n} = \frac{\theta^2}{3n}。
    lim⁡n→∞Var(θ^)=lim⁡n→∞θ23n=0\lim_{n \to \infty} Var(\hat{\theta}) = \lim_{n \to \infty} \frac{\theta^2}{3n} = 0。

由於 lim⁡n→∞E(θ^)=θ\lim_{n \to \infty} E(\hat{\theta}) = \theta 且 lim⁡n→∞Var(θ^)=0\lim_{n \to \infty} Var(\hat{\theta}) = 0,根據弱大數法則 (Weak Law of Large Numbers),θ^\hat{\theta} 是一致的。

結論: θ\theta 的矩估計量為 θ^=2Xˉ\hat{\theta} = 2\bar{X},且 θ^\hat{\theta} 是一致估計量。


(b) 針對 H0:θ=9H_0: \theta = 9 vs. H1:θ≠9H_1: \theta \neq 9,以 α=0.05\alpha = 0.05 建立檢定統計量與拒絕域。

我們使用的估計量是 θ^=2Xˉ\hat{\theta} = 2\bar{X}。
在虛無假設 H0H_0 下,我們假設 θ=9\theta = 9。
因此,在 H0H_0 為真時,
E(θ^)=θ=9E(\hat{\theta}) = \theta = 9
Var(θ^)=θ23n=923n=813n=27nVar(\hat{\theta}) = \frac{\theta^2}{3n} = \frac{9^2}{3n} = \frac{81}{3n} = \frac{27}{n}。

根據中央極限定理 (CLT),當 nn 夠大時,θ^=2Xˉ\hat{\theta} = 2\bar{X} 的分布近似於常態分布:
θ^≈N(E(θ^),Var(θ^))=N(9,27n)\hat{\theta} \approx N(E(\hat{\theta}), Var(\hat{\theta})) = N(9, \frac{27}{n})。

我們可以構造一個標準化的檢定統計量 ZZ:
Z=θ^−E(θ^)Var(θ^)=2Xˉ−927/nZ = \frac{\hat{\theta} - E(\hat{\theta})}{\sqrt{Var(\hat{\theta})}} = \frac{2\bar{X} - 9}{\sqrt{27/n}}
在 H0H_0 為真時,Z≈N(0,1)Z \approx N(0, 1)。

這是雙尾檢定,顯著水準為 α=0.05\alpha = 0.05。
我們需要找到標準常態分布的臨界值 zα/2z_{\alpha/2},使得 P(∣Z∣>zα/2)=αP(|Z| > z_{\alpha/2}) = \alpha。
對於 α=0.05\alpha = 0.05,α/2=0.025\alpha/2 = 0.025。
查標準常態分布表,累積機率為 1−0.025=0.9751 - 0.025 = 0.975 的 zz 值是 1.961.96。
所以,z0.025=1.96z_{0.025} = 1.96。

拒絕域:
我們在 ∣Z∣>zα/2|Z| > z_{\alpha/2} 時拒絕 H0H_0。
即 ∣2Xˉ−927/n∣>1.96| \frac{2\bar{X} - 9}{\sqrt{27/n}} | > 1.96。

檢定統計量:
T=2Xˉ−927/nT = \frac{2\bar{X} - 9}{\sqrt{27/n}}。
當 ∣T∣>1.96|T| > 1.96 時,拒絕 H0H_0。

注意: 題目只給了 n=3n=3 的樣本,這是一個非常小的樣本數。CLT 的近似可能不夠準確。然而,在沒有其他資訊或更精確的分布資訊時,通常會基於 CLT 來進行。如果 XiX_i 是常態分布,那麼 Xˉ\bar{X} 也是常態分布,此時 ZZ 的分布是精確的標準常態分布。但這裡 XiX_i 是均勻分布,所以 CLT 的應用是近似。


(c) 給定樣本 (X1,X2,X3)=(6.30,2.28,7.46)(X_1, X_2, X_3) = (6.30, 2.28, 7.46),建構 95% 信賴區間,並判斷是否拒絕 H0H_0。

步驟 1:計算樣本平均值 Xˉ\bar{X}。
n=3n = 3
Xˉ=6.30+2.28+7.463=16.043≈5.3467\bar{X} = \frac{6.30 + 2.28 + 7.46}{3} = \frac{16.04}{3} \approx 5.3467

步驟 2:計算 θ^\hat{\theta} 的值。
θ^=2Xˉ=2×5.3467≈10.6934\hat{\theta} = 2\bar{X} = 2 \times 5.3467 \approx 10.6934

步驟 3:建構 95% 信賴區間。
信賴區間的公式為:
θ^±zα/2Var(θ^)\hat{\theta} \pm z_{\alpha/2} \sqrt{Var(\hat{\theta})}
其中 Var(θ^)=θ23nVar(\hat{\theta}) = \frac{\theta^2}{3n}。
但是,信賴區間的建構通常是基於估計量本身的變異數。當我們建構 θ\theta 的信賴區間時,我們使用 θ^\hat{\theta} 的變異數。
我們知道 Var(θ^)=θ23nVar(\hat{\theta}) = \frac{\theta^2}{3n}。在建構信賴區間時,我們需要估計這個變異數。
一個常見的做法是將 θ\theta 用其估計量 θ^\hat{\theta} 取代。
所以,我們使用 θ^23n\frac{\hat{\theta}^2}{3n} 作為變異數的估計。
θ^23n=θ^3n\sqrt{\frac{\hat{\theta}^2}{3n}} = \frac{\hat{\theta}}{\sqrt{3n}}。

95% 信賴區間為:
θ^±zα/2θ^3n\hat{\theta} \pm z_{\alpha/2} \frac{\hat{\theta}}{\sqrt{3n}}
其中 zα/2=1.96z_{\alpha/2} = 1.96 (因為 α=0.05\alpha = 0.05)。
n=3n = 3。
θ^≈10.6934\hat{\theta} \approx 10.6934。

信賴區間 = 10.6934±1.96⋅10.69343×310.6934 \pm 1.96 \cdot \frac{10.6934}{\sqrt{3 \times 3}}
= 10.6934±1.96⋅10.6934910.6934 \pm 1.96 \cdot \frac{10.6934}{\sqrt{9}}
= 10.6934±1.96⋅10.6934310.6934 \pm 1.96 \cdot \frac{10.6934}{3}
= 10.6934±1.96⋅3.564510.6934 \pm 1.96 \cdot 3.5645
= 10.6934±6.986410.6934 \pm 6.9864

信賴區間的下限:10.6934−6.9864=3.707010.6934 - 6.9864 = 3.7070
信賴區間的上限:10.6934+6.9864=17.679810.6934 + 6.9864 = 17.6798

所以,θ\theta 的 95% 信賴區間是 (3.7070,17.6798)(3.7070, 17.6798)。

另一種建構信賴區間的方法:
有些教材會直接使用 θ^=2Xˉ\hat{\theta} = 2\bar{X} 的漸近分布 N(θ,θ23n)N(\theta, \frac{\theta^2}{3n})。
在建構信賴區間時,我們需要估計 θ23n\frac{\theta^2}{3n}。
一種方法是使用 θ^23n\frac{\hat{\theta}^2}{3n}。
另一種方法是使用 plug-in estimator for variance: θ^23n\frac{\hat{\theta}^2}{3n}。
θ^23n=θ^3n\sqrt{\frac{\hat{\theta}^2}{3n}} = \frac{\hat{\theta}}{\sqrt{3n}}。
這與上面計算的結果相同。

🔒

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

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

免費註冊

第 4 題15 分

Let (X1,…,Xn)(X_1, \dots, X_n) be a random sample from a continuous distribution with distribution function F(x)F(x) and probability density function f(x)f(x). Assume that FF is invertible. Define X(n)=max⁡(X1,…,Xn)X_{(n)} = \max(X_1, \dots, X_n), and Yn=n(1−F(X(n)))Y_n = n(1-F(X_{(n)})). Find the limiting distribution of YnY_n.

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

這一題的完整詳解

核心觀念

本題考查:

  • 機率積分轉換:若 FF 為連續分布函數,則 F(Xi)∼Uniform⁡(0,1)F(X_i)\sim \operatorname{Uniform}(0,1)。

  • 最大順序統計量 X(n)X_{(n)} 的分布:

    P(X(n)≤x)=F(x)n.P(X_{(n)}\le x)=F(x)^n.

  • 極限公式:

    (1−yn)n⟶e−y.\left(1-\frac{y}{n}\right)^n\longrightarrow e^{-y}.

題目中的 Yn=n(1−F(X(n)))Y_n=n(1-F(X_{(n)})) 衡量樣本最大值右側尾端機率的放大結果,因此其極限將呈現指數分布。

解題方法

由於 FF 可逆且連續,令

Ui=F(Xi),i=1,…,n.U_i=F(X_i),\qquad i=1,\dots,n.

根據機率積分轉換,

Ui∼Uniform⁡(0,1).U_i\sim \operatorname{Uniform}(0,1).

又因為 FF 為遞增函數,

F(X(n))=F(max⁡(X1,…,Xn))=max⁡(F(X1),…,F(Xn))=max⁡(U1,…,Un).F(X_{(n)})=F\left(\max(X_1,\dots,X_n)\right) =\max(F(X_1),\dots,F(X_n)) =\max(U_1,\dots,U_n).

記

U(n)=max⁡(U1,…,Un),U_{(n)}=\max(U_1,\dots,U_n),

則

Yn=n(1−U(n)).Y_n=n(1-U_{(n)}).

對任意固定的 y<0y<0,因為 Yn≥0Y_n\ge 0,所以

P(Yn≤y)=0.P(Y_n\le y)=0.

對任意固定的 y≥0y\ge 0,當 n>yn>y 時,

P(Yn≤y)=P(n(1−U(n))≤y)=P(U(n)≥1−yn)=1−P(U(n)<1−yn).\begin{aligned} P(Y_n\le y) &=P\left(n(1-U_{(n)})\le y\right)\\ &=P\left(U_{(n)}\ge 1-\frac{y}{n}\right)\\ &=1-P\left(U_{(n)}<1-\frac{y}{n}\right). \end{aligned}
🔒

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

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

免費註冊

第 5 題25 分

Consider the model
log⁡y=β0+β1x1+β2x2+u\log y = \beta_0 + \beta_1 x_1 + \beta_2 x_2 + u
E(u∣x)=0,u∣x∼N(0,σ2)E(u|x) = 0, \quad u|x \sim N(0, \sigma^2)
where log⁡y=log⁡(y)\log y = \log(y). We run the OLS regression of log⁡y\log y on x1x_1 and x2x_2, and obtain the predicted value of log⁡y\log y for any value of the independent variables:
log⁡y^=β^0+β^1x1+β^2x2\widehat{\log y} = \hat{\beta}_0 + \hat{\beta}_1 x_1 + \hat{\beta}_2 x_2
Let u^i\hat{u}_i be the residual term for i=1,…,Ni = 1, \dots, N.

(a) Find E(exp⁡(u))E(\exp(u)).
(b) Find E(y∣x)E(y|x).
(c) Find a consistent estimator for the predicted value of yy, y^\hat{y}, in terms of log⁡y^\widehat{\log y} and σ^2\hat{\sigma}^2, where σ^2\hat{\sigma}^2 is a consistent estimator for σ2\sigma^2.
(d) If we know that uu is not normally distributed, find a consistent estimator for the predicted value of yy, y^\hat{y}, in terms of log⁡y^\widehat{\log y} and u^i\hat{u}_i.

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

這一題的完整詳解

此題考驗對線性迴歸模型中,對數轉換與指數轉換的期望值計算,以及估計量的性質,包含一致性估計量。

核心觀念:

  1. 對數轉換: 當應變數對數轉換後,模型變成線性模型,但誤差項的分布會影響原始應變數的分布。
  2. 常態分布的性質: E(eX)E(e^X) 當 X∼N(μ,σ2)X \sim N(\mu, \sigma^2) 時,其值為 eμ+σ2/2e^{\mu + \sigma^2/2}。
  3. 條件期望值: E(Y∣X)E(Y|X) 的計算。
  4. OLS 估計量: 樣本迴歸線的預測值。
  5. 一致性估計量: 當樣本數趨近無限大時,估計量趨近真實參數。
  6. 殘差項: 實際值與預測值之間的差異。

模型設定:
log⁡y=β0+β1x1+β2x2+u\log y = \beta_0 + \beta_1 x_1 + \beta_2 x_2 + u
E(u∣x)=0,u∣x∼N(0,σ2)E(u|x) = 0, \quad u|x \sim N(0, \sigma^2)
其中 x=(x1,x2)x = (x_1, x_2)。
令 log⁡y=Y′\log y = Y'。則 Y′=β0+β1x1+β2x2+uY' = \beta_0 + \beta_1 x_1 + \beta_2 x_2 + u。
令 E(Y′∣x)=β0+β1x1+β2x2E(Y'|x) = \beta_0 + \beta_1 x_1 + \beta_2 x_2。
由 E(u∣x)=0E(u|x)=0 可知。
Var(Y′∣x)=Var(β0+β1x1+β2x2+u∣x)=Var(u∣x)=σ2Var(Y'|x) = Var(\beta_0 + \beta_1 x_1 + \beta_2 x_2 + u | x) = Var(u|x) = \sigma^2。


(a) Find E(exp⁡(u))E(\exp(u)).

我們已知 u∣x∼N(0,σ2)u|x \sim N(0, \sigma^2)。
根據常態分布的性質,E(eX)=eE(X)+Var(X)/2E(e^X) = e^{E(X) + Var(X)/2}。
在這裡,X=uX=u, E(u)=0E(u) = 0 且 Var(u)=σ2Var(u) = \sigma^2。
因此,
E(exp⁡(u)∣x)=E(eu∣x)=eE(u∣x)+Var(u∣x)/2E(\exp(u)|x) = E(e^u|x) = e^{E(u|x) + Var(u|x)/2}
=e0+σ2/2=eσ2/2= e^{0 + \sigma^2/2} = e^{\sigma^2/2}。

由於 E(exp⁡(u)∣x)E(\exp(u)|x) 對所有 xx 都是常數 eσ2/2e^{\sigma^2/2},所以其無條件期望值也為 eσ2/2e^{\sigma^2/2}。
E(exp⁡(u))=E[E(exp⁡(u)∣x)]=E[eσ2/2]=eσ2/2E(\exp(u)) = E[E(\exp(u)|x)] = E[e^{\sigma^2/2}] = e^{\sigma^2/2}。

【答案】 E(exp⁡(u))=eσ2/2E(\exp(u)) = e^{\sigma^2/2}。


(b) Find E(y∣x)E(y|x).

我們有 log⁡y=β0+β1x1+β2x2+u\log y = \beta_0 + \beta_1 x_1 + \beta_2 x_2 + u。
這意味著 y=exp⁡(log⁡y)=exp⁡(β0+β1x1+β2x2+u)y = \exp(\log y) = \exp(\beta_0 + \beta_1 x_1 + \beta_2 x_2 + u)。
y=exp⁡(β0+β1x1+β2x2)⋅exp⁡(u)y = \exp(\beta_0 + \beta_1 x_1 + \beta_2 x_2) \cdot \exp(u)。
令 E′=β0+β1x1+β2x2E' = \beta_0 + \beta_1 x_1 + \beta_2 x_2。
則 y=exp⁡(E′)⋅exp⁡(u)y = \exp(E') \cdot \exp(u)。

我們需要計算條件期望值 E(y∣x)E(y|x)。
E(y∣x)=E[exp⁡(E′)⋅exp⁡(u)∣x]E(y|x) = E[\exp(E') \cdot \exp(u) | x]
由於 exp⁡(E′)\exp(E') 是關於 xx 的函數,它對於條件期望值 E(⋅∣x)E(\cdot|x) 來說是常數。
E(y∣x)=exp⁡(E′)⋅E(exp⁡(u)∣x)E(y|x) = \exp(E') \cdot E(\exp(u)|x)
E(y∣x)=exp⁡(β0+β1x1+β2x2)⋅E(exp⁡(u)∣x)E(y|x) = \exp(\beta_0 + \beta_1 x_1 + \beta_2 x_2) \cdot E(\exp(u)|x)。

從 (a) 的結果,我們知道 E(exp⁡(u)∣x)=eσ2/2E(\exp(u)|x) = e^{\sigma^2/2}。
所以,
E(y∣x)=exp⁡(β0+β1x1+β2x2)⋅eσ2/2E(y|x) = \exp(\beta_0 + \beta_1 x_1 + \beta_2 x_2) \cdot e^{\sigma^2/2}
E(y∣x)=exp⁡(β0+β1x1+β2x2+σ22)E(y|x) = \exp(\beta_0 + \beta_1 x_1 + \beta_2 x_2 + \frac{\sigma^2}{2})。

【答案】 E(y∣x)=exp⁡(β0+β1x1+β2x2+σ22)E(y|x) = \exp(\beta_0 + \beta_1 x_1 + \beta_2 x_2 + \frac{\sigma^2}{2})。


(c) Find a consistent estimator for the predicted value of yy, y^\hat{y}, in terms of log⁡y^\widehat{\log y} and σ^2\hat{\sigma}^2, where σ^2\hat{\sigma}^2 is a consistent estimator for σ2\sigma^2.

預測值 y^\hat{y} 通常是指基於樣本估計出來的參數,代入模型計算出來的預測值。
根據 (b) 的結果,真實的條件期望值是 E(y∣x)=exp⁡(β0+β1x1+β2x2+σ22)E(y|x) = \exp(\beta_0 + \beta_1 x_1 + \beta_2 x_2 + \frac{\sigma^2}{2})。
將真實參數 (β0,β1,β2,σ2)(\beta_0, \beta_1, \beta_2, \sigma^2) 用其 OLS 估計量 (β^0,β^1,β^2,σ^2)(\hat{\beta}_0, \hat{\beta}_1, \hat{\beta}_2, \hat{\sigma}^2) 取代,我們可以得到 y^\hat{y} 的預測值。
y^=exp⁡(β^0+β^1x1+β^2x2+σ^22)\hat{y} = \exp(\hat{\beta}_0 + \hat{\beta}_1 x_1 + \hat{\beta}_2 x_2 + \frac{\hat{\sigma}^2}{2})。

題目要求用 log⁡y^\widehat{\log y} 和 σ^2\hat{\sigma}^2 來表示。
我們知道 log⁡y^=β^0+β^1x1+β^2x2\widehat{\log y} = \hat{\beta}_0 + \hat{\beta}_1 x_1 + \hat{\beta}_2 x_2 是 OLS 對 log⁡y\log y 的預測值。
所以,
y^=exp⁡(log⁡y^+σ^22)\hat{y} = \exp(\widehat{\log y} + \frac{\hat{\sigma}^2}{2})。

說明一致性:
OLS 估計量 β^0,β^1,β^2\hat{\beta}_0, \hat{\beta}_1, \hat{\beta}_2 是 β0,β1,β2\beta_0, \beta_1, \beta_2 的一致估計量,因為誤差項 uu 的條件期望為零且滿足其他標準假設。
σ^2\hat{\sigma}^2 是 σ2\sigma^2 的一致估計量。
log⁡y^\widehat{\log y} 是 E(log⁡y∣x)E(\log y|x) 的一致估計量。
由於 exp⁡(⋅)\exp(\cdot) 是一個連續函數,根據連續函數保持一致性 (Slutsky's theorem or properties of consistent estimators),exp⁡(log⁡y^)\exp(\widehat{\log y}) 是 exp⁡(E(log⁡y∣x))\exp(E(\log y|x)) 的一致估計量。
同理,exp⁡(log⁡y^+σ^22)\exp(\widehat{\log y} + \frac{\hat{\sigma}^2}{2}) 是 exp⁡(E(log⁡y∣x)+σ22)\exp(E(\log y|x) + \frac{\sigma^2}{2}) 的一致估計量。
而 exp⁡(E(log⁡y∣x)+σ22)=exp⁡(β0+β1x1+β2x2+σ22)=E(y∣x)\exp(E(\log y|x) + \frac{\sigma^2}{2}) = \exp(\beta_0 + \beta_1 x_1 + \beta_2 x_2 + \frac{\sigma^2}{2}) = E(y|x)。
所以,y^=exp⁡(log⁡y^+σ^22)\hat{y} = \exp(\widehat{\log y} + \frac{\hat{\sigma}^2}{2}) 是一個一致估計量。

【答案】 y^=exp⁡(log⁡y^+σ^22)\hat{y} = \exp(\widehat{\log y} + \frac{\hat{\sigma}^2}{2})。


(d) If we know that uu is not normally distributed, find a consistent estimator for the predicted value of yy, y^\hat{y}, in terms of log⁡y^\widehat{\log y} and u^i\hat{u}_i.

模型為:
log⁡y=β0+β1x1+β2x2+u\log y = \beta_0 + \beta_1 x_1 + \beta_2 x_2 + u
E(u∣x)=0E(u|x) = 0
log⁡y^=β^0+β^1x1+β^2x2\widehat{\log y} = \hat{\beta}_0 + \hat{\beta}_1 x_1 + \hat{\beta}_2 x_2
u^i=log⁡yi−log⁡y^i\hat{u}_i = \log y_i - \widehat{\log y}_i

我們需要找到 y^\hat{y} 的一致估計量。
從 y=exp⁡(log⁡y)=exp⁡(β0+β1x1+β2x2+u)y = \exp(\log y) = \exp(\beta_0 + \beta_1 x_1 + \beta_2 x_2 + u)。
E(y∣x)=E(exp⁡(β0+β1x1+β2x2+u)∣x)E(y|x) = E(\exp(\beta_0 + \beta_1 x_1 + \beta_2 x_2 + u)|x)
E(y∣x)=exp⁡(β0+β1x1+β2x2)⋅E(exp⁡(u)∣x)E(y|x) = \exp(\beta_0 + \beta_1 x_1 + \beta_2 x_2) \cdot E(\exp(u)|x)。

在這一部分,我們不知道 uu 的分布,只知道 E(u∣x)=0E(u|x)=0。
我們不能再使用 E(exp⁡(u)∣x)=eσ2/2E(\exp(u)|x) = e^{\sigma^2/2}。
我們需要用樣本資訊來估計 E(exp⁡(u)∣x)E(\exp(u)|x)。

一個一致的估計量 y^\hat{y} 應該基於 OLS 預測值 log⁡y^\widehat{\log y} 和對 E(exp⁡(u)∣x)E(\exp(u)|x) 的估計。

🔒

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

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

免費註冊

第 6 題10 分

Suppose that you are interested in estimating the ceteris paribus relationship between yy and x1x_1. For this purpose, you collect data on yy and x1x_1, and two control variables x2x_2 and x3x_3. Let β^1\hat{\beta}_1 be the simple regression estimate from yy on x1x_1, and let β^1.23\hat{\beta}_{1.23} be the multiple regression estimate from yy on x1x_1, x2x_2 and x3x_3.

(a) If x1x_1 is almost uncorrelated with x2x_2 and x3x_3 in the sample, and x2x_2 and x3x_3 are highly correlated, would you expect β^1\hat{\beta}_1 and β^1.23\hat{\beta}_{1.23} to be similar or very different? Explain.

(b) If x1x_1 is highly correlated with x2x_2 and x3x_3, and x2x_2 and x3x_3 have small partial effects on yy, would you expect se(β^1)se(\hat{\beta}_1) or se(β^1.23)se(\hat{\beta}_{1.23}) to be smaller? Explain. (se(.) is standard error)

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

這一題的完整詳解

此題考驗對迴歸分析中,多重共線性 (multicollinearity) 與控制變數 (control variables) 對估計量標準誤與係數估計的影響。

核心觀念:

  1. 簡單迴歸 (Simple Regression): y=β0+β1x1+uy = \beta_0 + \beta_1 x_1 + u。 β^1\hat{\beta}_1 是 yy 對 x1x_1 的簡單迴歸係數。
  2. 複迴歸 (Multiple Regression): y=β0+β1x1+β2x2+β3x3+vy = \beta_0 + \beta_1 x_1 + \beta_2 x_2 + \beta_3 x_3 + v。 β^1.23\hat{\beta}_{1.23} 是 yy 對 x1,x2,x3x_1, x_2, x_3 的複迴歸係數。
  3. Ceteris Paribus: 在其他條件不變的情況下。複迴歸的係數 β^1.23\hat{\beta}_{1.23} 試圖捕捉 x1x_1 對 yy 的 ceteris paribus 效應。
  4. 多重共線性: 自變數之間高度相關。
  5. 標準誤 (Standard Error, SE): 衡量迴歸係數估計量變異性的指標。較小的標準誤表示估計量較精確。

(a) If x1x_1 is almost uncorrelated with x2x_2 and x3x_3 in the sample, and x2x_2 and x3x_3 are highly correlated, would you expect β^1\hat{\beta}_1 and β^1.23\hat{\beta}_{1.23} to be similar or very different? Explain.

解釋:

  • 簡單迴歸 β^1\hat{\beta}_1: 估計的是 x1x_1 對 yy 的總效應,其中包含了 x1x_1 直接對 yy 的影響,以及 x1x_1 間接通過 x2x_2 和 x3x_3 對 yy 的影響(如果 x1x_1 與 x2,x3x_2, x_3 相關的話)。
  • 複迴歸 β^1.23\hat{\beta}_{1.23}: 估計的是 x1x_1 對 yy 的 ceteris paribus 效應,即在控制了 x2x_2 和 x3x_3 的影響後,x1x_1 對 yy 的直接影響。

題目給定:

  1. x1x_1 與 x2,x3x_2, x_3 幾乎不相關。
  2. x2x_2 與 x3x_3 高度相關。

分析:
當 x1x_1 與 x2,x3x_2, x_3 幾乎不相關時,在複迴歸模型中,將 x2,x3x_2, x_3 加入迴歸式,對 x1x_1 的係數 β^1.23\hat{\beta}_{1.23} 的影響會很小。
這是因為 x2,x3x_2, x_3 無法很好地解釋 x1x_1 的變異。
標準迴歸理論告訴我們,複迴歸係數 β^1.23\hat{\beta}_{1.23} 可以看作是 yy 對 x1x_1 的迴歸殘差,對 x2,x3x_2, x_3 的迴歸殘差進行迴歸的係數。
如果 x1x_1 與 x2,x3x_2, x_3 無關,那麼 x1x_1 對 yy 的影響,在簡單迴歸和複迴歸中,基本上是相同的。
也就是說,控制 x2,x3x_2, x_3 的動作,並不會顯著地改變 x1x_1 對 yy 的估計效應,因為 x2,x3x_2, x_3 並沒有「解釋」或「混淆」(confound) x1x_1 對 yy 的關係。

x2x_2 和 x3x_3 高度相關(即多重共線性)主要會影響 β^2\hat{\beta}_2 和 β^3\hat{\beta}_3 的標準誤,並可能使它們變得很大,但對於 x1x_1 的係數 β^1.23\hat{\beta}_{1.23} 的影響,取決於 x1x_1 與 x2,x3x_2, x_3 的關聯度。

結論:
因為 x1x_1 與 x2,x3x_2, x_3 幾乎不相關,所以 β^1\hat{\beta}_1 和 β^1.23\hat{\beta}_{1.23} 應該是相似的。
加入 x2x_2 和 x3x_3 作為控制變數,並不會顯著改變 x1x_1 對 yy 的 ceteris paribus 效應的估計。


(b) If x1x_1 is highly correlated with x2x_2 and x3x_3, and x2x_2 and x3x_3 have small partial effects on yy, would you expect se(β^1)se(\hat{\beta}_1) or se(β^1.23)se(\hat{\beta}_{1.23}) to be smaller? Explain.

解釋:

  • se(β^1)se(\hat{\beta}_1): 簡單迴歸中 x1x_1 係數的標準誤。
🔒

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

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

免費註冊

第 7 題15 分

Consider the probit model P(y=1∣z1,z2,q)=Φ(τ1z1+γ1z2+δ1q)P(y = 1|z_1, z_2, q) = \Phi(\tau_1 z_1 + \gamma_1 z_2 + \delta_1 q), where z2z_2 is a continuous variable, qq is independent of z1z_1 and z2z_2, and z1,z2,q∼N(0,1)z_1, z_2, q \sim N(0, 1). Moreover, z1z_1 and z2z_2 are observed but qq is not. Φ(⋅)\Phi(\cdot) is c.d.f. of N(0,1)N(0, 1).

(a) Find P(y=1∣z1,z2)P(y=1|z_1, z_2).
(b) Find P(y=1∣z1)P(y=1|z_1).

[Some useful N(0,1) probabilities]
P(N(0,1)≤2.33)=0.99P(N(0,1) \le 2.33) = 0.99, P(N(0,1)≤1.96)=0.975P(N(0,1) \le 1.96) = 0.975
P(N(0,1)≤1.64)=0.95P(N(0,1) \le 1.64) = 0.95, P(N(0,1)≤1.28)=0.90P(N(0,1) \le 1.28) = 0.90

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

這一題的完整詳解

此題考驗對機率模型 (Probit Model) 的理解,以及條件機率的計算,特別是涉及隱藏變數 (latent variable) 的情況。

核心觀念:

  1. Probit 模型: 應變數 yy 是一個二元變數 (0 或 1),其發生機率由一個標準常態分布的 CDF 給出。
  2. 條件機率: P(A∣B)P(A|B) 的計算。
  3. 隱藏變數 (Latent Variable): 模型中的 qq 是一個未被觀測到的變數。
  4. 常態分布的性質: 兩個獨立常態變數的線性組合也是常態變數。
  5. 全機率法則 (Law of Total Probability): P(A)=∑iP(A∣Bi)P(Bi)P(A) = \sum_i P(A|B_i) P(B_i)。

模型設定:
P(y=1∣z1,z2,q)=Φ(τ1z1+γ1z2+δ1q)P(y = 1|z_1, z_2, q) = \Phi(\tau_1 z_1 + \gamma_1 z_2 + \delta_1 q)
其中 z1,z2,q∼N(0,1)z_1, z_2, q \sim N(0, 1) 獨立。
z1,z2z_1, z_2 可觀測,qq 不可觀測。


(a) Find P(y=1∣z1,z2)P(y=1|z_1, z_2).

我們需要計算在給定可觀測變數 z1,z2z_1, z_2 的條件下,y=1y=1 的機率。
由於 qq 是未知的,我們需要將其邊緣化 (marginalize out)。
我們使用全機率法則:
P(y=1∣z1,z2)=∫−∞∞P(y=1∣z1,z2,q)fq(q)dqP(y=1|z_1, z_2) = \int_{-\infty}^{\infty} P(y=1|z_1, z_2, q) f_q(q) dq
其中 fq(q)f_q(q) 是 qq 的機率密度函數 (PDF)。
由於 q∼N(0,1)q \sim N(0, 1) 獨立,其 PDF 為 fq(q)=12πe−q2/2f_q(q) = \frac{1}{\sqrt{2\pi}} e^{-q^2/2}。

將模型代入:
P(y=1∣z1,z2)=∫−∞∞Φ(τ1z1+γ1z2+δ1q)12πe−q2/2dqP(y=1|z_1, z_2) = \int_{-\infty}^{\infty} \Phi(\tau_1 z_1 + \gamma_1 z_2 + \delta_1 q) \frac{1}{\sqrt{2\pi}} e^{-q^2/2} dq

令 μq∣z1,z2=τ1z1+γ1z2\mu_{q|z_1, z_2} = \tau_1 z_1 + \gamma_1 z_2。
則積分變為:
P(y=1∣z1,z2)=∫−∞∞Φ(μq∣z1,z2+δ1q)12πe−q2/2dqP(y=1|z_1, z_2) = \int_{-\infty}^{\infty} \Phi(\mu_{q|z_1, z_2} + \delta_1 q) \frac{1}{\sqrt{2\pi}} e^{-q^2/2} dq

這是一個期望值的形式。令 X=μq∣z1,z2+δ1qX = \mu_{q|z_1, z_2} + \delta_1 q。
則 P(y=1∣z1,z2)=E[Φ(X)]P(y=1|z_1, z_2) = E[\Phi(X)]。
由於 q∼N(0,1)q \sim N(0, 1), δ1q∼N(0,δ12)\delta_1 q \sim N(0, \delta_1^2)。
μq∣z1,z2\mu_{q|z_1, z_2} 是關於 z1,z2z_1, z_2 的常數。
所以,X=(τ1z1+γ1z2)+δ1qX = (\tau_1 z_1 + \gamma_1 z_2) + \delta_1 q 是關於 qq 的一個線性組合。
由於 q∼N(0,1)q \sim N(0, 1), XX 是一個常態隨機變數。
E(X)=E(μq∣z1,z2+δ1q)=μq∣z1,z2+δ1E(q)=τ1z1+γ1z2+δ1⋅0=τ1z1+γ1z2E(X) = E(\mu_{q|z_1, z_2} + \delta_1 q) = \mu_{q|z_1, z_2} + \delta_1 E(q) = \tau_1 z_1 + \gamma_1 z_2 + \delta_1 \cdot 0 = \tau_1 z_1 + \gamma_1 z_2。
Var(X)=Var(μq∣z1,z2+δ1q)=Var(δ1q)=δ12Var(q)=δ12⋅1=δ12Var(X) = Var(\mu_{q|z_1, z_2} + \delta_1 q) = Var(\delta_1 q) = \delta_1^2 Var(q) = \delta_1^2 \cdot 1 = \delta_1^2。
所以,X∼N(τ1z1+γ1z2,δ12)X \sim N(\tau_1 z_1 + \gamma_1 z_2, \delta_1^2)。

我們需要計算 E[Φ(X)]E[\Phi(X)],其中 X∼N(μ,σ2)X \sim N(\mu, \sigma^2)。
這個期望值 ∫−∞∞Φ(x)12πσe−(x−μ)22σ2dx\int_{-\infty}^{\infty} \Phi(x) \frac{1}{\sqrt{2\pi}\sigma} e^{-\frac{(x-\mu)^2}{2\sigma^2}} dx 的結果是 Φ(μ1+σ2)\Phi(\frac{\mu}{\sqrt{1+\sigma^2}})。
在這裡,μ=τ1z1+γ1z2\mu = \tau_1 z_1 + \gamma_1 z_2 且 σ2=δ12\sigma^2 = \delta_1^2。
所以,
P(y=1∣z1,z2)=Φ(τ1z1+γ1z21+δ12)P(y=1|z_1, z_2) = \Phi\left(\frac{\tau_1 z_1 + \gamma_1 z_2}{\sqrt{1+\delta_1^2}}\right)。

證明這個結果:
令 Z∼N(0,1)Z \sim N(0, 1) 且 X∼N(μ,σ2)X \sim N(\mu, \sigma^2) 獨立。
我們想計算 E[Φ(aX+b)]E[\Phi(aX+b)]。
令 W=aX+bW = aX+b。 W∼N(aμ+b,a2σ2)W \sim N(a\mu+b, a^2\sigma^2)。
E[Φ(W)]=∫−∞∞Φ(w)12πa2σ2e−(w−(aμ+b))22a2σ2dwE[\Phi(W)] = \int_{-\infty}^{\infty} \Phi(w) \frac{1}{\sqrt{2\pi a^2\sigma^2}} e^{-\frac{(w-(a\mu+b))^2}{2a^2\sigma^2}} dw。
這是一個標準結果:如果 X∼N(μ,σ2)X \sim N(\mu, \sigma^2),則 E[Φ(X)]=Φ(μ1+σ2)E[\Phi(X)] = \Phi(\frac{\mu}{\sqrt{1+\sigma^2}})。
在這裡,XX 是我們上面定義的 X=μq∣z1,z2+δ1qX = \mu_{q|z_1, z_2} + \delta_1 q。
μ=E(X)=τ1z1+γ1z2\mu = E(X) = \tau_1 z_1 + \gamma_1 z_2
σ2=Var(X)=δ12\sigma^2 = Var(X) = \delta_1^2
所以,P(y=1∣z1,z2)=Φ(τ1z1+γ1z21+δ12)P(y=1|z_1, z_2) = \Phi\left(\frac{\tau_1 z_1 + \gamma_1 z_2}{\sqrt{1+\delta_1^2}}\right)。

🔒

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

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

免費註冊

其他考古題