108 年 國立臺灣大學農業經濟學系碩士班《統計學(E)》

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

第 1 題12 分

  1. (12 points) Please solve the following questions:
    (a) (6 points) Bowl A contains three red and two white chips, and bowl B contains four red and three white chips. A chip is drawn at random from bowl A and transferred to bowl B. Compute the probability of then drawing a red chip from bowl B.
    (b) (6 points) An urn contains two red balls and four white balls. Sample successively five times at random and with replacement, so that the trials are independent. Compute the probability of each of the two sequences WWRWR and RWWWR.

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

這一題的完整詳解

核心觀念

本題考查:

  1. 全機率公式:若事件可依不同情況分解,總機率等於各情況機率乘上條件機率後加總。
  2. 條件機率:轉入紅色或白色籌碼後,碗 B 的籌碼總數與紅色籌碼數會改變。
  3. 獨立重複試驗:有放回抽樣使每次抽取彼此獨立,且每次紅、白球的機率固定。
  4. 特定序列的機率:直接將該序列中各結果的機率相乘,不需再乘排列數。

解題方法

(a) 從碗 B 抽到紅色籌碼的機率

碗 A 原有 3 個紅色、2 個白色籌碼,共 5 個。

抽出的籌碼有兩種情況。

情況一:從碗 A 轉入紅色籌碼

轉入紅色籌碼的機率為

P(轉入紅)=35.P(\text{轉入紅})=\frac{3}{5}.

此時碗 B 原有 4 個紅色、3 個白色籌碼。轉入一個紅色籌碼後,碗 B 共有:

  • 紅色籌碼:4+1=54+1=5 個
  • 白色籌碼:33 個
  • 總數:88 個

因此從碗 B 抽到紅色籌碼的條件機率為

P(B抽紅∣轉入紅)=58.P(\text{B抽紅}\mid \text{轉入紅})=\frac{5}{8}.

情況二:從碗 A 轉入白色籌碼

轉入白色籌碼的機率為

P(轉入白)=25.P(\text{轉入白})=\frac{2}{5}.

此時碗 B 共有:

  • 紅色籌碼:44 個
  • 白色籌碼:3+1=43+1=4 個
  • 總數:88 個

因此

P(B抽紅∣轉入白)=48.P(\text{B抽紅}\mid \text{轉入白})=\frac{4}{8}.

依全機率公式,

🔒

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

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

免費註冊

第 2 題13 分

  1. (13 points) Let X1,X2,X3X_1, X_2, X_3 denote a random sample of size n=3n = 3 from a distribution with the geometric pmf
    f(x)=(14)(34)x−1f(x) = \left(\frac{1}{4}\right) \left(\frac{3}{4}\right)^{x-1}, x=1,2,3,……x = 1, 2, 3, \dots\dots
    (a) (3 points) Compute P(X1=1,X2=3,X3=3)P(X_1 = 1, X_2 = 3, X_3 = 3)
    (b) (5 points) Determine P(X1+X2+X3=6)P(X_1 + X_2 + X_3 = 6)
    (c) (5 points) If YY equals the maximum of X1,X2,X3X_1, X_2, X_3, find P(Y≤2)=P(X1≤2)P(X2≤2)P(X3≤2)P(Y \le 2) = P(X_1 \le 2) P(X_2 \le 2) P(X_3 \le 2).

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

這一題的完整詳解

本題考驗幾何分佈的性質,包括機率計算、隨機變數和的機率以及最大值的機率。

給定的機率質量函數 (pmf) 是幾何分佈:
f(x)=(14)(34)x−1f(x) = \left(\frac{1}{4}\right) \left(\frac{3}{4}\right)^{x-1}, x=1,2,3,…x = 1, 2, 3, \dots
這表示成功的機率 p=14p = \frac{1}{4},失敗的機率 q=1−p=34q = 1-p = \frac{3}{4}。
幾何分佈定義為:第一次成功出現在第 xx 次試驗的機率。

(a)
計算 P(X1=1,X2=3,X3=3)P(X_1 = 1, X_2 = 3, X_3 = 3)。
由於 X1,X2,X3X_1, X_2, X_3 是來自同一個分佈的隨機樣本,它們是獨立同分佈的 (i.i.d.)。
因此,聯合機率等於各機率的乘積:
P(X1=1,X2=3,X3=3)=P(X1=1)×P(X2=3)×P(X3=3)P(X_1 = 1, X_2 = 3, X_3 = 3) = P(X_1 = 1) \times P(X_2 = 3) \times P(X_3 = 3)

我們需要計算 f(1)f(1), f(3)f(3):
f(1)=(14)(34)1−1=(14)(34)0=14f(1) = \left(\frac{1}{4}\right) \left(\frac{3}{4}\right)^{1-1} = \left(\frac{1}{4}\right) \left(\frac{3}{4}\right)^0 = \frac{1}{4}
f(3)=(14)(34)3−1=(14)(34)2=14×916=964f(3) = \left(\frac{1}{4}\right) \left(\frac{3}{4}\right)^{3-1} = \left(\frac{1}{4}\right) \left(\frac{3}{4}\right)^2 = \frac{1}{4} \times \frac{9}{16} = \frac{9}{64}

所以,
P(X1=1,X2=3,X3=3)=f(1)×f(3)×f(3)P(X_1 = 1, X_2 = 3, X_3 = 3) = f(1) \times f(3) \times f(3)
P(X1=1,X2=3,X3=3)=(14)×(964)×(964)P(X_1 = 1, X_2 = 3, X_3 = 3) = \left(\frac{1}{4}\right) \times \left(\frac{9}{64}\right) \times \left(\frac{9}{64}\right)
P(X1=1,X2=3,X3=3)=8116384P(X_1 = 1, X_2 = 3, X_3 = 3) = \frac{81}{16384}

【答案】8116384\frac{81}{16384}

(b)
決定 P(X1+X2+X3=6)P(X_1 + X_2 + X_3 = 6)。
這是一個關於三個幾何分佈隨機變數和的問題。
一個常見的方法是考慮其機率母函數 (Probability Generating Function, PGF) 或直接列舉所有可能的組合。
幾何分佈 X∼Geom(p)X \sim Geom(p) 的 PGF 是 GX(z)=E[zX]=∑x=1∞zxp(1−p)x−1=pz∑x=1∞(z(1−p))x−1=pz1−(1−p)zG_X(z) = E[z^X] = \sum_{x=1}^\infty z^x p (1-p)^{x-1} = pz \sum_{x=1}^\infty (z(1-p))^{x-1} = \frac{pz}{1 - (1-p)z}。
對於 X1,X2,X3X_1, X_2, X_3 獨立同分佈,它們的和 S=X1+X2+X3S = X_1 + X_2 + X_3 的 PGF 是各 PGF 的乘積:
GS(z)=(GX(z))3=(pz1−(1−p)z)3G_S(z) = (G_X(z))^3 = \left(\frac{pz}{1 - (1-p)z}\right)^3
其中 p=1/4p = 1/4, 1−p=3/41-p = 3/4。
GS(z)=(14z1−34z)3=(z4−3z)3G_S(z) = \left(\frac{\frac{1}{4}z}{1 - \frac{3}{4}z}\right)^3 = \left(\frac{z}{4 - 3z}\right)^3

我們需要找到 P(S=6)P(S=6)。這可以通過對 GS(z)G_S(z) 進行泰勒展開來找到 z6z^6 的係數。
另一種方法是考慮 XiX_i 的定義:XiX_i 是第一次成功前失敗的次數加一。
如果我們定義 Yi=Xi−1Y_i = X_i - 1,則 YiY_i 是失敗的次數,服從 Yi∼Geom0(p)Y_i \sim Geom_0(p),其 pmf 為 P(Yi=k)=(1−p)kpP(Y_i=k) = (1-p)^k p, k=0,1,2,…k=0,1,2,\dots。
那麼 X1+X2+X3=6X_1 + X_2 + X_3 = 6 等價於 (Y1+1)+(Y2+1)+(Y3+1)=6(Y_1+1) + (Y_2+1) + (Y_3+1) = 6,即 Y1+Y2+Y3=3Y_1 + Y_2 + Y_3 = 3。
Y1,Y2,Y3Y_1, Y_2, Y_3 是獨立的、服從參數為 p=1/4p=1/4 的零起始幾何分佈 (Negative Binomial distribution with r=1r=1 trials)。
Y1+Y2+Y3Y_1 + Y_2 + Y_3 服從負二項分佈 NB(r=3,p=1/4)NB(r=3, p=1/4),其 pmf 為:
P(∑Yi=k)=(k+r−1r−1)pr(1−p)kP(\sum Y_i = k) = \binom{k+r-1}{r-1} p^r (1-p)^k
在這裡 r=3r=3, k=3k=3, p=1/4p=1/4, 1−p=3/41-p=3/4。

🔒

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

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

免費註冊

第 3 題6 分

  1. (6 points) Let μ\mu and σ2\sigma^2 denote the mean and variance of the random variable XX. Determine
    (a) (3 points) E[X−μσ]E\left[\frac{X-\mu}{\sigma}\right]
    (b) (3 points) E[(X−μσ)2]E\left[\left(\frac{X-\mu}{\sigma}\right)^2\right]

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

這一題的完整詳解

本題考驗期望值的線性性質和變異數的定義。

令 Z=X−μσZ = \frac{X-\mu}{\sigma}。我們知道 ZZ 是將隨機變數 XX 標準化後的結果。
標準化隨機變數的期望值為 0,變異數為 1。

(a)
我們需要計算 E[X−μσ]E\left[\frac{X-\mu}{\sigma}\right]。
利用期望值的線性性質:
E[X−μσ]=E[1σX−μσ]E\left[\frac{X-\mu}{\sigma}\right] = E\left[\frac{1}{\sigma}X - \frac{\mu}{\sigma}\right]
=1σE[X]−μσ= \frac{1}{\sigma} E[X] - \frac{\mu}{\sigma}
由於 μ=E[X]\mu = E[X],所以:
=1σ(μ)−μσ=μσ−μσ=0= \frac{1}{\sigma} (\mu) - \frac{\mu}{\sigma} = \frac{\mu}{\sigma} - \frac{\mu}{\sigma} = 0。

這也符合標準化變數的性質,即 E[Z]=0E[Z] = 0。

【答案】00

🔒

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

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

免費註冊

第 4 題9 分

  1. (9 points) If the moment-generating function of XX is M(t)=15et+25e2t+25e3tM(t) = \frac{1}{5}e^t + \frac{2}{5}e^{2t} + \frac{2}{5}e^{3t}, find the mean, variance, and pmf of XX.

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

這一題的完整詳解

本題考驗機率母函數 (Moment-Generating Function, MGF) 的應用,透過 MGF 來找出隨機變數的機率分佈、期望值和變異數。

給定的 MGF 是 M(t)=15et+25e2t+25e3tM(t) = \frac{1}{5}e^t + \frac{2}{5}e^{2t} + \frac{2}{5}e^{3t}。
MGF 的一般形式是 M(t)=E[etX]=∑etxP(X=x)M(t) = E[e^{tX}] = \sum e^{tx} P(X=x) (對於離散型) 或 ∫etxf(x)dx\int e^{tx} f(x) dx (對於連續型)。
觀察 MGF 的形式,它是由幾個指數項的線性組合構成,每一項的係數都是正數且加總為 1:
15+25+25=55=1\frac{1}{5} + \frac{2}{5} + \frac{2}{5} = \frac{5}{5} = 1。
這表明 XX 是一個離散型隨機變數,其可能的取值是 tt 的指數項的指數部分。

1. 找出機率質量函數 (pmf)
MGF 的形式 M(t)=∑icieaitM(t) = \sum_{i} c_i e^{a_i t} 表明隨機變數 XX 的可能取值為 aia_i,且對應的機率為 P(X=ai)=ciP(X=a_i) = c_i。
在此,我們有:
M(t)=15e1⋅t+25e2⋅t+25e3⋅tM(t) = \frac{1}{5}e^{1 \cdot t} + \frac{2}{5}e^{2 \cdot t} + \frac{2}{5}e^{3 \cdot t}
因此,隨機變數 XX 的可能取值是 1, 2, 3。
其對應的機率是:
P(X=1)=15P(X=1) = \frac{1}{5}
P(X=2)=25P(X=2) = \frac{2}{5}
P(X=3)=25P(X=3) = \frac{2}{5}

所以,pmf 為:
f(x)={15if x=125if x=225if x=30otherwisef(x) = \begin{cases} \frac{1}{5} & \text{if } x=1 \\ \frac{2}{5} & \text{if } x=2 \\ \frac{2}{5} & \text{if } x=3 \\ 0 & \text{otherwise} \end{cases}

2. 計算期望值 (Mean)
期望值 E[X]E[X] 可以透過 MGF 的性質求得:E[X]=M′(0)E[X] = M'(0)。
首先,計算 M(t)M(t) 對 tt 的一階導數 M′(t)M'(t):
M′(t)=ddt(15et+25e2t+25e3t)M'(t) = \frac{d}{dt} \left(\frac{1}{5}e^t + \frac{2}{5}e^{2t} + \frac{2}{5}e^{3t}\right)
M′(t)=15et+25(2e2t)+25(3e3t)M'(t) = \frac{1}{5}e^t + \frac{2}{5}(2e^{2t}) + \frac{2}{5}(3e^{3t})
M′(t)=15et+45e2t+65e3tM'(t) = \frac{1}{5}e^t + \frac{4}{5}e^{2t} + \frac{6}{5}e^{3t}
現在,將 t=0t=0 代入 M′(t)M'(t):
E[X]=M′(0)=15e0+45e0+65e0=15+45+65=115E[X] = M'(0) = \frac{1}{5}e^0 + \frac{4}{5}e^0 + \frac{6}{5}e^0 = \frac{1}{5} + \frac{4}{5} + \frac{6}{5} = \frac{11}{5}。

或者,也可以直接根據 pmf 計算期望值:
E[X]=∑xP(X=x)=1×15+2×25+3×25=15+45+65=115E[X] = \sum x P(X=x) = 1 \times \frac{1}{5} + 2 \times \frac{2}{5} + 3 \times \frac{2}{5} = \frac{1}{5} + \frac{4}{5} + \frac{6}{5} = \frac{11}{5}。

🔒

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

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

免費註冊

第 5 題10 分

  1. (10 points) Let f(x;θ)=θxθ−1f(x;\theta) = \theta x^{\theta-1}, 0<x<10 < x < 1, 0<θ<∞0 < \theta < \infty. Let X1,X2,…,XnX_1, X_2, \dots, X_n denote a random sample of size nn from this distribution.
    (a) (4 points) Sketch the pdf of XX for (i) θ=1/2\theta = 1/2, (ii) θ=1\theta = 1.
    (b) (6 points) Show that θ^=−n∑i=1nln⁡Xi\hat{\theta} = -\frac{n}{\sum_{i=1}^n \ln X_i} is the maximum likelihood estimator of θ\theta.

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

這一題的完整詳解

本題考驗機率密度函數 (pdf) 的圖形描繪以及最大概似估計量 (Maximum Likelihood Estimator, MLE) 的推導。

給定的機率密度函數 (pdf) 為 f(x;θ)=θxθ−1f(x;\theta) = \theta x^{\theta-1},定義域為 0<x<10 < x < 1,參數 θ\theta 的範圍為 0<θ<∞0 < \theta < \infty。

(a)
繪製 pdf 的圖形。
pdf 的形式為 θxθ−1\theta x^{\theta-1}。我們需要對 θ\theta 的不同值進行分析。

(i) θ=1/2\theta = 1/2
當 θ=1/2\theta = 1/2,pdf 為 f(x;1/2)=12x12−1=12x−12=12xf(x; 1/2) = \frac{1}{2} x^{\frac{1}{2}-1} = \frac{1}{2} x^{-\frac{1}{2}} = \frac{1}{2\sqrt{x}},定義域為 0<x<10 < x < 1。
分析此函數的特性:

  • 當 x→0+x \to 0^+, f(x)→∞f(x) \to \infty。
  • 當 x=1x = 1, f(1)=121=12f(1) = \frac{1}{2\sqrt{1}} = \frac{1}{2}。
  • 此函數是遞減的。
  • 我們需要驗證其積分是否為 1 (作為 pdf):
    ∫0112x−1/2dx=12[2x1/2]01=12(21−0)=1\int_0^1 \frac{1}{2} x^{-1/2} dx = \frac{1}{2} \left[ 2x^{1/2} \right]_0^1 = \frac{1}{2} (2\sqrt{1} - 0) = 1。
    積分確實為 1,所以這是一個有效的 pdf。
    圖形會從 (0,∞)(0, \infty) 開始,遞減至 (1,1/2)(1, 1/2)。

(ii) θ=1\theta = 1
當 θ=1\theta = 1,pdf 為 f(x;1)=1⋅x1−1=x0=1f(x; 1) = 1 \cdot x^{1-1} = x^0 = 1,定義域為 0<x<10 < x < 1。
這是一個均勻分佈 U(0,1)U(0, 1) 的 pdf。

  • 函數值恆為 1。
  • 圖形是一條水平線 y=1y=1 在 0<x<10 < x < 1 的區間。
  • 積分驗證:∫011dx=[x]01=1−0=1\int_0^1 1 dx = [x]_0^1 = 1 - 0 = 1。

圖形描述:

  • 對於 θ=1/2\theta = 1/2: 圖形從 x=0x=0 處無限高開始,然後遞減,在 x=1x=1 處的值為 1/21/2。
  • 對於 θ=1\theta = 1: 圖形是一條從 x=0x=0 到 x=1x=1 的水平線,高度為 1。

(b)
推導最大概似估計量 (MLE) θ^\hat{\theta}。
給定一個隨機樣本 X1,X2,…,XnX_1, X_2, \dots, X_n,其似然函數 (Likelihood function) L(θ)L(\theta) 是所有樣本點 pdf 的乘積:
L(θ)=∏i=1nf(Xi;θ)=∏i=1nθXiθ−1L(\theta) = \prod_{i=1}^n f(X_i;\theta) = \prod_{i=1}^n \theta X_i^{\theta-1}
L(θ)=θn∏i=1nXiθ−1L(\theta) = \theta^n \prod_{i=1}^n X_i^{\theta-1}
L(θ)=θn(∏i=1nXi)θ−1L(\theta) = \theta^n \left(\prod_{i=1}^n X_i\right)^{\theta-1}

🔒

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

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

免費註冊

第 6 題50 分

  1. (50 points) Consider the standard simple regression model: y=β0+β1x+uy = \beta_0 + \beta_1 x + u under the Gauss-Markov Assumptions. The OLS estimators for β1\beta_1 and β0\beta_0 are
    β^1=∑i=1n(xi−xˉ)(yi−yˉ)∑i=1n(xi−xˉ)2=SxySxx\hat{\beta}_1 = \frac{\sum_{i=1}^n (x_i - \bar{x})(y_i - \bar{y})}{\sum_{i=1}^n (x_i - \bar{x})^2} = \frac{S_{xy}}{S_{xx}} and β^0=yˉ−β^1xˉ\hat{\beta}_0 = \bar{y} - \hat{\beta}_1 \bar{x}, where Sxx=∑i=1n(xi−xˉ)2S_{xx} = \sum_{i=1}^n (x_i - \bar{x})^2.
    The residuals are defined as: u^i=yi−β^0−β^1xi\hat{u}_i = y_i - \hat{\beta}_0 - \hat{\beta}_1 x_i, for i=1,2,…,ni = 1, 2, \dots, n.
    From normal equation we know ∑i=1nxiu^i=0\sum_{i=1}^n x_i \hat{u}_i = 0.
    (a) (5 points) Suppose that E(u)=α≠0E(u) = \alpha \ne 0. Show that the model can always rewritten with the same slope, but a new intercept and error, where the new error has a zero expected value.
    (b) (5 points) What are the five Gauss-Markov assumptions for the simple regression?
    (c) (10 points) Show that β^1=ρxyσ^yσ^x\hat{\beta}_1 = \rho_{xy} \frac{\hat{\sigma}_y}{\hat{\sigma}_x}, where ρxy\rho_{xy} is the sample correlation and σ^y\hat{\sigma}_y and σ^x\hat{\sigma}_x denote the sample standard deviations.
    (d) (10 points) Show that R2R^2 is equal to the square of the sample correlation coefficient between yy and y^\hat{y}.
    (e) (20 points) Let c1>0c_1 > 0 and c2>0c_2 > 0 be constants. Regress c1yic_1 y_i on c2xic_2 x_i. Find the new intercept and slope estimators in terms of β1\beta_1 and β0\beta_0. (10 points) Again, regress (yi/c1)(y_i/c_1) on (xi−c2)(x_i - c_2). What will be the new intercept and slope estimators? (10 points)

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

這一題的完整詳解

核心觀念

本題綜合考查簡單線性迴歸的基本性質:

  • 誤差項平均數的重新定義。
  • Gauss–Markov 五項假設。
  • OLS 斜率與樣本相關係數的關係。
  • R2R^2 與相關係數的關係。
  • 變數平移與尺度轉換對 OLS 估計量的影響。

記

Sxx=∑i=1n(xi−xˉ)2,Sxy=∑i=1n(xi−xˉ)(yi−yˉ),S_{xx}=\sum_{i=1}^n(x_i-\bar{x})^2,\qquad S_{xy}=\sum_{i=1}^n(x_i-\bar{x})(y_i-\bar{y}),

並令

Syy=∑i=1n(yi−yˉ)2.S_{yy}=\sum_{i=1}^n(y_i-\bar{y})^2.

(a) 誤差項平均數不為零時的改寫

已知

y=β0+β1x+u,E(u)=α≠0.y=\beta_0+\beta_1x+u,\qquad E(u)=\alpha\ne 0.

定義新的誤差項

u∗=u−α.u^*=u-\alpha.

則

E(u∗)=E(u−α)=E(u)−α=α−α=0.E(u^*)=E(u-\alpha)=E(u)-\alpha=\alpha-\alpha=0.

又因為 u=u∗+αu=u^*+\alpha,代回原模型:

y=β0+β1x+u∗+α.y=\beta_0+\beta_1x+u^*+\alpha.

整理為

y=(β0+α)+β1x+u∗.y=(\beta_0+\alpha)+\beta_1x+u^*.

因此,令新的截距為

β0∗=β0+α,\beta_0^*=\beta_0+\alpha,

即可寫成

y=β0∗+β1x+u∗,\boxed{y=\beta_0^*+\beta_1x+u^*},

其中

E(u∗)=0.\boxed{E(u^*)=0}.

斜率仍然是原本的 β1\beta_1,只有截距增加 α\alpha。


(b) 簡單迴歸的五項 Gauss–Markov 假設

1. 參數具有線性形式

模型對未知參數必須是線性的:

yi=β0+β1xi+ui.y_i=\beta_0+\beta_1x_i+u_i.

此處「線性」是指對 β0\beta_0、β1\beta_1 線性,而不是要求變數本身只能呈現直線形式。

2. 隨機抽樣

觀察值

{(xi,yi):i=1,…,n}\{(x_i,y_i):i=1,\dots,n\}

來自母體的隨機樣本,且各觀察值通常假設彼此獨立、具有相同分配。

3. 解釋變數具有樣本變異

必須有

Sxx=∑i=1n(xi−xˉ)2>0.S_{xx}=\sum_{i=1}^n(x_i-\bar{x})^2>0.

也就是 xix_i 不能全部相同,否則 β^1=Sxy/Sxx\hat{\beta}_1=S_{xy}/S_{xx} 無法定義。

4. 零條件平均假設

E(ui∣xi)=0.E(u_i\mid x_i)=0.

此假設表示在給定 xix_i 後,誤差項的平均值仍為零,因此 xix_i 不會與誤差項系統性相關。這是保證 OLS 估計量不偏的重要條件。

5. 同質變異數

Var⁡(ui∣xi)=σ2,\operatorname{Var}(u_i\mid x_i)=\sigma^2,

其中 σ2\sigma^2 為固定常數。此假設保證 OLS 在所有線性不偏估計量中具有最小變異數,即 OLS 是 BLUE。


(c) 證明 β^1=ρxyσ^yσ^x\hat{\beta}_1=\rho_{xy}\dfrac{\hat{\sigma}_y}{\hat{\sigma}_x}

樣本相關係數定義為

ρxy=SxySxxSyy.\rho_{xy}=\frac{S_{xy}}{\sqrt{S_{xx}S_{yy}}}.

因此

Sxy=ρxySxxSyy.S_{xy}=\rho_{xy}\sqrt{S_{xx}S_{yy}}.

OLS 斜率為

β^1=SxySxx.\hat{\beta}_1=\frac{S_{xy}}{S_{xx}}.

代入相關係數的表示式:

β^1=ρxySxxSyySxx=ρxySyySxx.\hat{\beta}_1 = \frac{\rho_{xy}\sqrt{S_{xx}S_{yy}}}{S_{xx}} = \rho_{xy}\sqrt{\frac{S_{yy}}{S_{xx}}}.

若樣本標準差採用相同的分母,例如皆採用 n−1n-1,則

σ^x=Sxxn−1,σ^y=Syyn−1.\hat{\sigma}_x=\sqrt{\frac{S_{xx}}{n-1}}, \qquad \hat{\sigma}_y=\sqrt{\frac{S_{yy}}{n-1}}.

所以

σ^yσ^x=SyySxx.\frac{\hat{\sigma}_y}{\hat{\sigma}_x} = \sqrt{\frac{S_{yy}}{S_{xx}}}.

因此得到

β^1=ρxyσ^yσ^x.\boxed{ \hat{\beta}_1 = \rho_{xy}\frac{\hat{\sigma}_y}{\hat{\sigma}_x} }.

這表示迴歸斜率的符號與 xx、yy 的樣本相關係數相同;斜率的絕對大小還會受到兩者尺度,即標準差比例的影響。


(d) 證明 R2R^2 等於 yy 與 y^\hat{y} 的樣本相關係數平方

由迴歸式

y^i=β^0+β^1xi\hat{y}_i=\hat{\beta}_0+\hat{\beta}_1x_i

以及 β^0=yˉ−β^1xˉ\hat{\beta}_0=\bar{y}-\hat{\beta}_1\bar{x},可得

y^i−yˉ=β^1(xi−xˉ).\hat{y}_i-\bar{y} = \hat{\beta}_1(x_i-\bar{x}).

因此,yy 與 y^\hat{y} 的樣本相關係數為

ρy,y^=∑i=1n(yi−yˉ)(y^i−yˉ)∑i=1n(yi−yˉ)2∑i=1n(y^i−yˉ)2.\rho_{y,\hat{y}} = \frac{\sum_{i=1}^n(y_i-\bar{y})(\hat{y}_i-\bar{y})} {\sqrt{ \sum_{i=1}^n(y_i-\bar{y})^2 \sum_{i=1}^n(\hat{y}_i-\bar{y})^2 }}.

先計算分子:

∑i=1n(yi−yˉ)(y^i−yˉ)=β^1∑i=1n(yi−yˉ)(xi−xˉ)=β^1Sxy.\sum_{i=1}^n(y_i-\bar{y})(\hat{y}_i-\bar{y}) = \hat{\beta}_1 \sum_{i=1}^n(y_i-\bar{y})(x_i-\bar{x}) = \hat{\beta}_1S_{xy}.

由 β^1=Sxy/Sxx\hat{\beta}_1=S_{xy}/S_{xx},可得

β^1Sxy=Sxy2Sxx.\hat{\beta}_1S_{xy} = \frac{S_{xy}^2}{S_{xx}}.

再計算 y^\hat{y} 的離均差平方和:

∑i=1n(y^i−yˉ)2=β^12Sxx=Sxy2Sxx.\sum_{i=1}^n(\hat{y}_i-\bar{y})^2 = \hat{\beta}_1^2S_{xx} = \frac{S_{xy}^2}{S_{xx}}.

因此

ρy,y^=Sxy2/SxxSyy⋅Sxy2/Sxx=∣Sxy∣SxxSyy=∣ρxy∣.\rho_{y,\hat{y}} = \frac{S_{xy}^2/S_{xx}} {\sqrt{ S_{yy}\cdot S_{xy}^2/S_{xx} }} = \frac{|S_{xy}|}{\sqrt{S_{xx}S_{yy}}} = |\rho_{xy}|.

所以

ρy,y^2=ρxy2.\rho_{y,\hat{y}}^2=\rho_{xy}^2.

另一方面,簡單迴歸中的決定係數為

🔒

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

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

免費註冊

其他考古題