高级计量经济学 10:大样本OLS(下)

高级计量经济学 10:大样本OLS(下)

此文内容为《高级计量经济学及STATA应用》的笔记,陈强老师著,高等教育出版社出版。

我只将个人会用到的知识作了笔记,并对教材较难理解的部分做了进一步阐述。为了更易于理解,我还对教材上的一些部分(包括代码和正文)做了修改。

仅供学习参考,请勿转载,侵删!


目录

  • 5 大样本OLS
    • 5.7 大样本OLS的假定
      • 5.7.1 线性假定
      • 5.7.2 渐近独立的平稳过程
      • 5.7.3 前定解释变量
      • 5.7.4 秩条件
      • 5.7.5 鞅差分序列
      • 5.7.6 解释变量四阶矩存在
    • 5.8 OLS的大样本性质
      • 5.8.1 \pmb b 为一致估计量
      • 5.8.2 \pmb b 位渐近正态
      • 5.8.3 {\rm Avar}(\pmb b) 的一致估计
    • 5.9 线性假设的大样本检验
      • 5.9.1 单个系数的检验
      • 5.9.2 线性假设的检验

\S \text{ 第 5 章 } \S

\text{大样本OLS}


5 大样本OLS

5.7 大样本OLS的假定

5.7.1 线性假定

y_i = \pmb x_i^\prime \pmb \beta + \varepsilon_i

5.7.2 渐进独立的平稳过程(ergodic stationarity)

(K+1) 维随机过程 \{y_i,\pmb x_i\} 是渐进独立的平稳过程

5.7.3 前定解释变量(predetermined regressors)

所有解释变量均前定(predetermined),这就意味着他们与同期的扰动项正交,即 {\rm E}(x_{ik}\varepsilon_i)=0,\forall i,k 。由于方程通常都有常数项,所以总可以假设 {\rm E}(\varepsilon_i)=0 ,故:
{\rm Cov}(x_{ik}\varepsilon_i) = {\rm E}(x_{ik}\varepsilon_i) - {\rm E}(x_{ik})\underbrace{{\rm E}(\varepsilon_i)}_{=0} = 0 - 0 = 0

这个 {\rm E}(\varepsilon_i)=0 的假设跟严格外生性 {\rm E}(\varepsilon_i|{\bf X})=0 是不一样的,后者可以推出前者

这意味着 x_{ik} 与同期的扰动项 \varepsilon_i 不相关,就好像 \varepsilon_i 产生以前, x_{ik} 就已经被确定了一样,故名前定解释变量。如果定义:
\pmb g_i \equiv \pmb x_i \varepsilon_i = \left( \begin{array}{c} x_{i1} \\ \vdots\\x_{iK} \end{array} \right) \varepsilon_i
那么我们就有:
{\rm E}(\pmb g_i) \equiv {\rm E}(\pmb x_i \varepsilon_i)= \pmb 0
这一点比小样本的严格外生性的要求更低,因为后者要求 \varepsilon_i 与过去、现在、未来的所有 \pmb x_j,\forall j = 1,\cdots,n 都不相关;而现在只要求 \varepsilon_i 与现在的 \pmb x_i 不相关即可!

5.7.4 秩条件

K \times K 阶矩阵 {\rm E}(\pmb x_i \pmb x_i^\prime) 为非退化矩阵,从而 \left[{\rm E}(\pmb x_i \pmb x_i^\prime)\right]^{-1} 存在,这保证了 \pmb S_{XX}^{-1} 存在

5.7.5 鞅差分序列

要求 \pmb g_i 是一个鞅差分序列,而且其协方差矩阵
\pmb S \equiv {\rm Cov}(\pmb g_i) \equiv{\rm E}(\pmb g_i \pmb g_i^\prime)-\underbrace{{\rm E}(\pmb g_i){\rm E}(\pmb g_i^\prime)}_{=\pmb 0} = {\rm E}(\varepsilon^2_i\pmb x_i \pmb x_i^\prime)
是一个非退化矩阵。

这个假设比 5.7.3 的假设更强。前定解释变量保证了 \pmb b 是一致估计量,而鞅差分序列则进一步地保证了其渐进正态

5.7.6 解释变量的四阶矩存在

这是一个数学技巧的假定,不用太在意。即\forall i,j,k{\rm E}[x_{ik}x_{ij}] 存在且有限

我们发现,在上面的假定中,我们相比于小样本OLS,放松了:

  • 严格外生性 -> 前定解释变量
  • 正态随机扰动项 -> 独立平稳、鞅差分序列

这比小样本OLS的条件放宽了太多


5.8 OLS的大样本性质

为了便于使用渐近理论,我们需把估计量 \pmb b 写成 \{y_i, \pmb x_i\} 的函数。由于:
\pmb b = ({\bf X'X})^{-1}{\bf X'}\pmb y
于是我们把 {\bf X'X} 划为一组,{\bf X}'\pmb y 划为另外一组,并定义:
\pmb b = \left(\frac{{\bf X'X}}{n}\right)^{-1} \frac{{\bf X'}\pmb y}{n} = \pmb S_{XX}^{-1} S_{Xy}
其中,由于数据矩阵 \bf X 的结构是:
{\bf X} \equiv \left( \begin{array}{c} \pmb x_1^\prime \\ \pmb x_2^\prime \\ \vdots\\\pmb x_n^\prime \end{array} \right)
所以 \pmb S_{XX} 就是:
\pmb S_{XX} \equiv \frac{1}{n} {\bf X'X} = \frac{1}{n} (\pmb x_1 \cdots \pmb x_n) \cdot \left( \begin{array}{c} \pmb x_1^\prime \\ \vdots\\\pmb x_n^\prime \end{array} \right) = \frac{1}{n} \sum_{i=1}^n \pmb x_i \pmb x_i^\prime
同理 \pmb S_{Xy} 就是:
\pmb S_{Xy} \equiv \frac{1}{n} {\bf X}^\prime \pmb y =\frac{1}{n} (\pmb x_1 \cdots \pmb x_n) \cdot \left( \begin{array}{c} y_1 \\ \vdots\\ y_n \end{array} \right) =\frac{1}{n} \sum_{i=1}^n \pmb x_i y_i

5.8.1 \pmb b 为一致估计量

在假设 5.7.1 到 5.7.4 (不包括 \pmb g_i 是鞅差分序列的情况下)下,有:
\mathop{\rm plim}\limits_{n \to \infty} \pmb b = \pmb \beta

证明:抽样误差可以写为
\pmb b - \pmb \beta = A \pmb \varepsilon = ({\bf X'X})^{-1}{\bf X'}\pmb \varepsilon = \left(\frac{\sum\limits_{i=1}^n \pmb x_i \pmb x_i^\prime}{n}\right)^{-1} \frac{{\bf X'}\pmb \varepsilon}{n} = \pmb S_{XX}^{-1} \bar{\pmb g}
其中, \bar{\pmb g} = \frac{{\bf X'}\pmb \varepsilon}{n}\pmb g = (\pmb g_1 \cdots \pmb g_n) 的“均值”。我们的目标是证明 \pmb S_{XX}^{-1} \bar{\pmb g} \stackrel{p}\longrightarrow \pmb 0 ,所以证明的思路是分别求 \pmb S_{XX}^{-1}\bar{\pmb g} 分别依概率收敛到什么。

假定 5.7.2 意味着随机序列 \{\pmb x_i \pmb x_i^\prime\} 也是一个渐进独立的平稳序列,所以根据渐进独立定理:

如果 \{x_i\}_{i=1}^\infty 渐进独立且平稳,那么样本矩依概率收敛到总体矩

\pmb S_{XX} \equiv \frac{1}{n} \sum_{i=1}^n \pmb x_i \pmb x_i^\prime \stackrel{p}\longrightarrow {\rm E}(\pmb x_i \pmb x_i^\prime )

假定 5.7.4 意味着 \left[{\rm E}(\pmb x_i \pmb x_i^\prime )\right]^{-1} 存在,所以:
\pmb S_{XX}^{-1} \stackrel{p}\longrightarrow \left[{\rm E}(\pmb x_i \pmb x_i^\prime )\right]^{-1}
另外,由于 \pmb g_i \equiv \pmb x_i\varepsilon_i = \pmb x_i - (y_i - \pmb x_i \pmb \beta)\pmb x_iy_i 的函数。而 \{y_i, \pmb x_i\} 是一个渐进独立的平稳序列,从而 \{\pmb g_i\}_{i=1}^\infty 也是一个渐进独立的平稳序列,于是假定 5.7.3 意味着:
\bar{\pmb g} \equiv \frac{1}{n}\sum_{i=1}^n \pmb g_i \stackrel{p}\longrightarrow {\rm E}(\pmb g_i) \equiv {\rm E}(\pmb x_i \varepsilon_i)= \pmb 0
所以我们发现:
\pmb b - \pmb \beta = \pmb S_{XX}^{-1} \bar{\pmb g} \stackrel{p}\longrightarrow \left[{\rm E}(\pmb x_i \pmb x_i^\prime )\right]^{-1} \cdot \pmb 0 = \pmb 0
于是:
\mathop{\rm plim}\limits_{n \to \infty} (\pmb b -\pmb \beta) = \pmb 0 \Rightarrow \mathop{\rm plim}\limits_{n \to \infty} \pmb b =\pmb \beta
证毕。

5.8.2 \pmb b 为渐近正态

如果把 假定5.7.3{\rm E}(\pmb g_i) = 0)强化为 假定5.7.5\{\pmb g_i\}是\rm MDS),那么 \pmb b 是渐近正态分布:
\sqrt{n}(\pmb b - \pmb \beta) \stackrel2cccccccc8\longrightarrow N(0,{\rm Avar}(\pmb b))
其中 {\rm Avar}(\pmb b)\pmb b 的渐近方差:
{\rm Avar}(\pmb b) = \left[{\rm E}(\pmb x_i \pmb x_i^\prime )\right]^{-1} \pmb S \left[{\rm E}(\pmb x_i \pmb x_i^\prime )\right]^{-1}
其中 \pmb S\pmb g_i 的协方差矩阵:
\pmb S \equiv {\rm Cov}(\pmb g_i) \equiv{\rm E}(\pmb g_i \pmb g_i^\prime)-\underbrace{{\rm E}(\pmb g_i){\rm E}(\pmb g_i^\prime)}_{=\pmb 0} = {\rm E}(\varepsilon^2_i\pmb x_i \pmb x_i^\prime)

证明:由于我们知道:
\pmb{b-\beta} = \pmb S_{XX}^{-1} \bar{\pmb g}
所以
\sqrt{n}(\pmb{b-\beta}) = \pmb S_{XX}^{-1}(\sqrt{n} \bar{\pmb g})
根据 假定5.7.5\pmb g_i 是鞅差分序列,于是根据鞅差分序列的中心极限定理,我们有:

渐近独立平稳鞅差分序列 \{\pmb g_i\}_{i=1}^\infty样本均值\sqrt{n} 的速度依分布收敛为均值为 \pmb 0 协方差矩阵为 {\rm Cov}(\pmb g_i) 的正态分布

\bar{\pmb g} = \frac{1}{n}\sum_{i=1}^n \pmb g_i \Rightarrow \sqrt{n} \bar{\pmb g} \stackrel2cccccccc8\longrightarrow N(\pmb 0,\pmb S)

由于 \sqrt{n}(\pmb{b-\beta}) = \pmb S_{XX}^{-1}(\sqrt{n} \bar{\pmb g})\sqrt{n} \bar{\pmb g} 的线性组合,所以它也会以 \sqrt{n} 的速度依分布收敛为正态分布:
\sqrt{n}(\pmb{b-\beta}) = \pmb S_{XX}^{-1}(\sqrt{n} \bar{\pmb g}) \stackrel2cccccccc8\longrightarrow N(\pmb0,{\rm Avar}(\pmb b))
根据渐近正态的定义,那么 \sqrt{n}(\pmb{b-\beta}) 是渐近正态的。接下来我们计算 {\rm Avar}(\pmb b)

由于 \sqrt{n}(\pmb{b-\beta}) = \pmb S_{XX} (\sqrt{n}\bar{\pmb g}) ,因此:
{\rm Avar}(\pmb b) = {\rm Var}\left[\sqrt{n}(\pmb{b-\beta})\right] = {\rm Var}\left[(\pmb S_{XX}(\sqrt{n}\bar{\pmb g})\right]
使用夹心估计量,
{\rm Var}\left[(\pmb S_{XX}(\sqrt{n}\bar{\pmb g})\right]) = \pmb S_{XX} {\rm Var}(\sqrt{n}\bar{\pmb g}) \pmb S_{XX}^\prime
那么在 n\to \infty 时,\pmb S_{XX}{\rm Var}(\sqrt{n}\bar{\pmb g}) 分别收敛:
{\rm Avar}(\pmb) = {\rm E}(\pmb x_i \pmb x_i^\prime ) S [{\rm E}(\pmb x_i \pmb x_i^\prime )]^\prime
由于 {\rm E}(\pmb x_i \pmb x_i^\prime ) 是对称阵,所以 [{\rm E}(\pmb x_i \pmb x_i^\prime )]^\prime = {\rm E}(\pmb x_i \pmb x_i^\prime ) ,于是:
\sqrt{n}(\pmb b - \pmb \beta) \stackrel2cccccccc8\longrightarrow N(0,{\rm E}(\pmb x_i \pmb x_i^\prime ) \cdot\pmb S \cdot {\rm E}(\pmb x_i \pmb x_i^\prime ))
证毕。

我们发现,在这里不需要假设扰动项服从正态分布。这里的思路是:

  • 我们假设 \pmb g_i = \pmb x_i \varepsilon_i 是渐近独立的平稳的鞅差分序列(假定5.7.5
  • 那么 \sqrt{n}\bar{\pmb g} 依分布收敛到正态分布(鞅差分序列的中心极限定理)(这里是正态分布的来源)
  • 建立起 \sqrt{n}\pmb(b-\beta) 与的 \sqrt{n}\bar{\pmb g} 线性关系,那么 \sqrt{n}\pmb(b-\beta) 也依分布收敛到正态分布(正态分布的线性组合也是正态分布)(把要证明的目标绑到一个正态分布上)
  • \sqrt{n}\pmb(b-\beta) 也依分布收敛到正态分布本身就是 \pmb b 渐近正态的定义
  • 下一步要计算 \sqrt{n}\pmb(b-\beta) 的方差
  • 使用夹心估计量,并在 n\to \infty 分别收敛
  • 得证

5.8.3 {\rm Avar}(\pmb b) 的一致估计

如果存在 \hat{\pmb S}\pmb S 的一致估计量,那么 \pmb S_{XX}^{-1} \hat{\pmb S} \pmb S_{XX}^{-1}{\rm Avar}(\pmb b) 的一致估计,即:
\mathop{\rm plim}\limits_{n \to \infty} \widehat{{\rm Avar}}(\pmb b) = \mathop{\rm plim}\limits_{n \to \infty} \pmb S_{XX}^{-1} \hat{\pmb S} \pmb S_{XX}^{-1} = {\rm Avar}(\pmb b)

证明:如果存在 \hat{\pmb S} \stackrel{p}\longrightarrow \pmb S ,已知 \pmb S_{XX}^{-1} \equiv \frac{1}{n} \sum_{i=1}^n \pmb x_i \pmb x_i^\prime \stackrel{p}\longrightarrow \left[{\rm E}(\pmb x_i \pmb x_i^\prime )\right]^{-1} ,那么很自然地就有:
\pmb S_{XX}^{-1} \hat{\pmb S} \pmb S_{XX}^{-1} \stackrel{p}\longrightarrow \left[{\rm E}(\pmb x_i \pmb x_i^\prime )\right]^{-1} \pmb S \left[{\rm E}(\pmb x_i \pmb x_i^\prime )\right]^{-1}
证毕。

为了求解 \pmb S 的一致估计量,我们要运用 假定5.7.6 ,从而可以证明(这个证明比较难,略):
\mathop{\rm plim}\limits_{n \to \infty}\hat{\pmb S} = \mathop{\rm plim}\limits_{n \to \infty}\frac{1}{n}\sum_{i=1}^n e_i^2\pmb x_i \pmb x_i^\prime = {\rm E}(\varepsilon^2_i\pmb x_i \pmb x_i^\prime)=\pmb S
\hat{\pmb S}\pmb S 的一致估计量,其中 \{e_i\}_{i=1}^n 为最小二乘法的残差??梢越徊街っ?s^2\sigma^2 的一致估计。

性质s^2 是无条件方差 {\rm E}(\varepsilon^2_i)=\sigma^2 的一致估计量

证明:由有关消灭矩阵的知识(见《高级计量5》3.3):
s^2 \equiv \frac{\pmb{e'e}}{n-K} = \frac{\pmb{\varepsilon'M\varepsilon}}{n-K} = \frac{\pmb{\varepsilon}^\prime \left[{\bf I}_n - {\bf X(X'X)^{-1}X'}\right] \pmb{\varepsilon}}{n-K}
继续计算下去:
\begin{split} 原式 &= \frac{1}{n-K} \left[\pmb{\varepsilon}^\prime \pmb{\varepsilon} - \pmb{\varepsilon}^\prime {\bf X(X'X)^{-1}X'}\pmb{\varepsilon} \right]\\ 提一个n出来&=\frac{n}{n-K} \left[\frac{\pmb{\varepsilon}^\prime \pmb{\varepsilon}}{n} - \frac{\pmb{\varepsilon}^\prime {\bf X}}{n} \left(\frac{\bf X'X}{n}\right)^{-1} \frac{{\bf X'}\pmb{\varepsilon}}{n} \right]\\ 用定义&=\frac{n}{n-K}\left[ \frac{1}{n} \sum_{i=1}^n \varepsilon^2 - \bar{\pmb g}^\prime \pmb S_{XX}^{-1} \bar{\pmb g}\right] \quad (\star) \end{split}
由于我们知道:
\mathop{\rm plim}\limits_{n \to \infty} \frac{n}{n-K} = 1 ,\quad \mathop{\rm plim}\limits_{n \to \infty} \frac{1}{n} \sum_{i=1}^n \varepsilon_i^2 = {\rm E}(\sigma^2), \quad \mathop{\rm plim}\limits_{n \to \infty} \bar{\pmb g}^\prime \pmb S_{XX}^{-1} \bar{\pmb g} = \pmb 0 \cdot {\rm E}(\pmb x_i \pmb x_i^\prime ) \cdot \pmb 0 = 0
其中,第二个等式运用了 \{\varepsilon_i\}_{i=1}^n 是渐近独立的平稳列。 对 (\star)n \to \infty 就有:
\mathop{\rm plim}\limits_{n \to \infty} s^2 = \mathop{\rm plim}\limits_{n \to \infty} \frac{n}{n-K}\left[ \frac{1}{n} \sum_{i=1}^n \varepsilon^2 - \bar{\pmb g}^\prime \pmb S_{XX}^{-1} \bar{\pmb g}\right] = 1 \cdot \left[ \sigma^2 - 0 \right] = \sigma^2
s^2\sigma^2 的一致估计量。

证毕。


5.9 线性假设的大样本检验

5.9.1 单个系数的检验

在原假设 H_0:\beta_k = \bar{\beta}_k 成立的条件下,我们有:
\sqrt{n}(b_k - \bar\beta)k \stackrel2cccccccc8\longrightarrow N(0,{\rm Avar}(b_k))
其中,b_k\pmb b 的第 k 个元素,而 {\rm Avar}(b_k) 是协方差矩阵 {\rm Avar}(\pmb b) 的第 (k,k) 个元素。由于 {\rm Avar}(\pmb b) 我们是未知的,所以考虑用它的一致估计 \widehat{\rm Avar}(\pmb b) \stackrel{p}\longrightarrow {\rm Avar}(\pmb b) 来替代它。那么我们可以定义 t 统计量为:

注意,虽然叫 t 统计量,实际上服从的是渐近正态!

t_k \equiv \frac{\sqrt{n}(b_k-\bar{\beta}_k)}{\sqrt{\widehat{\rm Avar}(b_k)}} \stackrel2cccccccc8\longrightarrow N(0,1)

这里的思想是简单的,回忆对任意 X\sim(\mu,\sigma^2) ,都有 \frac{X-\mu}{\sigma}\sim N(0,1) ,其实两者的结构是一样的。

当然,为了便于描述,我们可以对 t_k 做进一步变形:
t_k \equiv \frac{\sqrt{n}(b_k-\bar{\beta}_k)}{\sqrt{\widehat{\rm Avar}(b_k)}} = \frac{(b_k-\bar{\beta}_k)}{\sqrt{\frac{1}{n}\widehat{\rm Avar}(b_k)}} = \frac{b_k-\bar{\beta}_k}{{\rm SE}^\star(b_k)}\stackrel2cccccccc8 \longrightarrow N(0,1)
那么我们就定义:
{\rm SE}^\star(b_k)\equiv \sqrt{\frac{1}{n}\widehat{{\rm Avar}}(b_k)} = \sqrt{\frac{1}{n}\left(\pmb S_{XX}^{-1} \pmb{\hat S}\pmb S_{XX}^{-1} \right)_{kk}}
异方差稳健的标准误(heteroskedasticiry-consistent standard erros),简称稳健的标准误(robust standard errors)。之所以这么称呼,是因为推导的过程中没有用到条件同方差的假定,所以在条件异方差下也能用。

性质:可以证明,在条件同方差的假设下,稳健的标准误退化为普通标准误

证明:如果假设 {\rm E}(\varepsilon_i^2|\pmb x_i)=\sigma^2>0 (条件同方差),那么使用迭代期望定律:
\pmb S \equiv {\rm E}(\pmb x_i^\prime \pmb x_i \varepsilon^2) = {\rm E}_{\pmb x_i} {\rm E}(\pmb x_i^\prime \pmb x_i \varepsilon^2|\pmb x_i)= {\rm E}_{\pmb x_i}[\pmb x_i^\prime \pmb x_i {\rm E}( \varepsilon^2|\pmb x_i)] = {\rm E}_{\pmb x_i}(\sigma^2 \pmb x_i^\prime \pmb x_i) = \sigma^2 {\rm E}(\pmb x_i^\prime \pmb x_i)
我们发现 s^2 \stackrel{p}\longrightarrow \sigma^2 ,\pmb S_{XX} \stackrel{p}\longrightarrow {\rm E}(\pmb x_i \pmb x_i^\prime) ,所以:
\mathop{\rm plim}\limits_{n \to \infty} \hat{\pmb S} =\pmb S=\sigma^2 {\rm E}(\pmb x_i^\prime \pmb x_i)=\mathop{\rm plim}\limits_{n \to \infty} s^2\pmb S_{XX}
于是:
\begin{split} \mathop{\rm plim}\limits_{n \to \infty} \widehat{\rm Avar}(\pmb b) &= \mathop{\rm plim}\limits_{n \to \infty} \pmb S_{XX}^{-1} \cdot \mathop{\rm plim}\limits_{n \to \infty} \pmb{\hat S}\cdot \mathop{\rm plim}\limits_{n \to \infty}\pmb S_{XX}^{-1} \\ &= \mathop{\rm plim}\limits_{n \to \infty} \pmb S_{XX}^{-1}\cdot \mathop{\rm plim}\limits_{n \to \infty}(s^2\pmb S_{XX})\cdot\mathop{\rm plim}\limits_{n \to \infty} \pmb S_{XX}^{-1}\\ &= \mathop{\rm plim}\limits_{n \to \infty} \pmb S_{XX}^{-1}\cdot (s^2\pmb S_{XX})\cdot\pmb S_{XX}^{-1}\\ &=\mathop{\rm plim}\limits_{n \to \infty} ns^2({\bf X'X})^{-1} \end{split}

写这么多 \mathop{\rm plim}\limits_{n \to \infty} 一方面是为了装逼

但更重要的是让你看明白这个代换是怎么来的??

一定一定不要错过 \mathop{\rm plim}\limits_{n \to \infty} 的细节!

于是:
\widehat{\rm Avar}(\pmb b) = ns^2({\bf X'X})^{-1}
所以:
{\rm SE}^\star(b_k)\equiv \sqrt{\frac{1}{n}\widehat{{\rm Avar}}(b_k)} = \sqrt{\frac{1}{n}ns^2({\bf X'X})^{-1}} = \sqrt{s^2({\bf X'X})^{-1}}
这就是普通的标准误了

证毕。

5.9.2 线性假设的检验

思路与小样本OLS是类似的,对零假设 H_0: \pmb{R\beta=r} ,我们实际上要衡量 \pmb{Rb-r}\pmb 0 的距离??梢灾っ鳎臣屏?W 满足:
W = \left[\sqrt{n}(\pmb{Rb-r}) \right]^\prime \left[\pmb R \widehat{\rm Avar}(\pmb b) \pmb R^\prime \right]^{-1} \left[\sqrt{n}(\pmb{Rb-r}) \right] \stackrel2cccccccc8\longrightarrow \chi^2(m)

证明:令 \pmb c_n \equiv \sqrt{n}(\pmb{Rb-r}) ,令 \pmb Q_n \equiv \pmb R \widehat{\rm Avar}(\pmb b) \pmb R^\prime ,那么现在 W = \pmb{ c_n^\prime Q_n^{-1}c_n}

如果 H_0: \pmb{R\beta=r} 成立,那么
\pmb c_n \equiv \sqrt{n}(\pmb{Rb-r}) = \sqrt{n}(\pmb{Rb-R\beta}) = \sqrt{n}\pmb R(\pmb{b-\beta}) = \pmb R[\sqrt{n}(\pmb{b-\beta})]
由于我们知道 \sqrt{n}(\pmb b - \pmb \beta) \stackrel2cccccccc8\longrightarrow N(0,{\rm Avar}(\pmb b)) ,而 \pmb c_n\sqrt{n}(\pmb b - \pmb \beta) 的线性组合,所以 \pmb c_n 也是依分布收敛于正态分布的。而且:
{\rm Var}(\pmb c_n) = {\rm Var}( \pmb R[\sqrt{n}(\pmb{b-\beta})] ) = \pmb R'{\rm Var}( \sqrt{n}(\pmb{b-\beta}) )\pmb R
\pmb c_n \stackrel2cccccccc8\longrightarrow \pmb c ,从而在 n \to \infty 时:
\pmb c \sim N(\pmb 0, \pmb R'{\rm Avar}( \pmb b) \pmb R)
定义 \pmb Q \equiv \pmb R {\rm Avar}(\pmb b) \pmb R^\prime ,由于 \widehat{{\rm Avar}}(\pmb b) \stackrel{p}\longrightarrow {\rm Avar}(\pmb b) ,所以 \pmb Q_n \stackrel{p}\longrightarrow \pmb Q

于是:
W = \pmb{ c_n^\prime Q_n^{-1}c_n} \stackrel2cccccccc8\longrightarrow \pmb{c'Q^{-1}c}=\pmb{c'}[{\rm Var}(\pmb c)]^{-1} \pmb c \sim \chi^2(m)
证毕。

?著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 214,128评论 6 493
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,316评论 3 388
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事?!?“怎么了?”我有些...
    开封第一讲书人阅读 159,737评论 0 349
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,283评论 1 287
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,384评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,458评论 1 292
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,467评论 3 412
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,251评论 0 269
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,688评论 1 306
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,980评论 2 328
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,155评论 1 342
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,818评论 4 337
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,492评论 3 322
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,142评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,382评论 1 267
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,020评论 2 365
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,044评论 2 352