肖仕维,赵培信,2

1.重庆工商大学 数学与统计学院,重庆 400067 2.经济社会应用统计重庆市重点实验室,重庆 400067

空间自回归模型[1]凭借其对于数据空间自相关性的有效解释,能够很好的弥补经典自回归模型在空间数据违背独立性假设情况下的数据解释问题。由于数据多元化时代的来临,该模型已经在统计学领域受到了广泛的关注,并积极的应用于经济政策分析的各类领域,特别是实产和房地产经济、环境与资源经济等[2]此类呈现强烈自相关性的空间数据。而部分线性空间自回归模型[3-4]通过引入未知非参数项,能够更加灵活的解释某些潜在协变量对于响应变量的影响,其模型结构为：

Y=ρWY+Xβ+m(Z)+ε。

(1)

其中,Y=(y1,y2,…yn)T表示响应变量的n×1维观测向量;W表示预先设定的n×n的空间权重矩阵;WY表示响应变量Y的空间滞后项;X=(x1,x2,…xn)T表示协变量的n×q的观测矩阵,其中xi=(xi1,xi2,…,xiq)T为q×1维列向量,i=1,2,…,n;m(Z)=(m(z1),m(z2),…m(zn))T表示n×1的未知函数向量,其中Z=(z1,z2,…zn)T表示某个潜在变量的n×1维观测向量,m(·)是一个未知光滑函数;ρ表示空间自回归参数且满足|ρ|<1;β=(β1,β2,…βq)T表示q×1的未知参数向量。此外,ε=(ε1,ε2,…εn)T是一个n×1的残差向量其中εi,i=1,2,…,n相互独立同分布。

目前对于此模型未知函数m(Z)有两种主流估计方法,一是由Fan和Gijbels[5]提出的局部多项式拟合来对m(Z)进行非参数估计;二是由Schumaker[6]提出的B样条(B-spline)逼近技术对m(Z)进行非参数估计。

在对模型(1)的检验研究过程中,本文则采用Schumaker提出的B样条逼近技术对m(Z)进行处理。对模型(1)的估计问题目前已有部分文献进行了研究,比如,Su和Jin[7]利用剖面最大拟似然估计方法研究了模型(1)的估计问题。Koch和Krisztin[8]结合B样条逼近技术和遗传算法研究了模型(1)的估计问题。Su[9]则利用广义矩估计方法(GMM)研究了模型(1)的估计问题。

但是关于模型(1)的检验问题,目前还没有相关文献进行系统研究。为此,本文将利用广义似然比检验方法研究模型(1)的检验问题。具体地,本文首先结合B样条逼近方法以及矩阵的正交投影技术消除模型的非参数分量对模型参数分量检验精度的影响,进而构造了一个关于模型参数分量的广义似然比检验统计量。另外由于模型结构的复杂性,理论上很难得到检验统计量的渐近分布。为此,本文利用Bootstrap[10]自助抽样法给出了一个计算检验统计量p值的迭代算法。

1 基于B样条投影矩阵的模型参数估计

B样条是数学的子学科数值分析里样条曲线一种特殊的表示形式,其分片多项式的构造、保持一定的连续性等性质使其拟合性能优越。具体地,记t0,t1,…tkn为区间[0,1]内的节点,并记π(t)=(B1(t),B2(t),…Bkn+l+1(t))T为一组l阶B样条基函数,则非参数函数m(t)可渐近表示为：

(2)

其中α=(α1,α2,…αkn+l+1)T为样条系数向量。结合式(2)可得:

m(Z)=(π(z1),π(z2),…,π(zn))Tα≜Πα。

进而模型(1)可以转化为：

Y=ρWY+Xβ+Πα+ε+e,

(3)

其中e=m(Z)-Πα表示B样条拟合误差。令P=Π(ΠTΠ)-1ΠT表示Π的正交投影矩阵,则(In-P)Π=Π-PΠ=Π-Π=0。在模型(3)两边同乘以矩阵S=In-P,则模型(3)可以转化为：

SY=S(ρWY+Xβ)+S(ε+e),

(4)

(5)

注意到式(5)为经典线性模型,利用最小二乘方法可得如下估计：

(6)

接下来考虑σ2的估计参数,设ε+e服从零均值同方差的正态分布,可得关于σ2的正态拟对数似然函数为：

(7)

最大化式(7)可得σ2的估计为：

(8)

2 基于似然函数的模型检验及Bootstrap实现

考虑如下参数显著性检验问题：

H0:Aθ=bvsH1:Aθ≠b,

(9)

其中,A为给定的d×(q+1)的列满秩矩阵,b为一个已知的d×1维向量。此检验的目的是检验空间自回归系数ρ和回归参数向量β是否显著。

进而可得在原假设H0下的最大拟对数似然函数为：

(10)

(11)

因此结合式(10)和式(11),并利用与Fan和Jiang[11]类似的方法,可定义广义似然比检验统计量为：

(12)

相应地,检验p值的计算为：

p=PH0(T>t),

(13)

其中,PH0表示在原假设H0下的累积概率,t表示检验统计量T的观测值。对任一给定的检验显著性水平α,如果p<α则拒绝原假设,否则就不能拒绝原假设。但是由于模型结构的复杂性,在计算p值时,检验统计量T的渐进分布很难直接的得到。因此本文采用Bootstrap自助法对检验p值的计算给出一个迭代算法,具体步骤如下：

(14)

第三步：由式(14)得到一个新的因变量向量Y*,则通过自助样本{Y*,X,Z}可以得到一个新的检验统计量T的自助版本的检验统计量的值T*。

(15)

其中#{}表示集合中元素的个数。