本文考虑半参数变系数部分线性模型,其一般形式为 Y=XTα(U)+ZTβ+ε,(1.1) 其中 α(⋅)=(α1(⋅),⋯,αp(⋅))T 为 p×1 未知系数函数向量, β=(β1,⋯,βq)T 是 q×1 维的未知参数向量, X,Z 以及 U 为协变量,Y 为响应变量, ε 为模型误差且满足 E(ε|X,Z,U)=0, Var(ε|X,Z,U)=σ2. 为了避免``维数灾祸", 一般假设 U 是单变量. 在模型 (1.1) 中,当 α(⋅)≡α 为某个未知参数向量时, 该模型即为线性模型; 当 p=1 且 X=1 时, 该模型即为部分线性模型; 当 Z=0 时,该模型即为常见的变系数模型.
关于模型 (1.1) 的统计推断, 在独立数据和纵向数据下已有许多学者对其进行了研究. Fan 和Huang[1] 提出了profile 最小二乘估计并且基于广义似然比检验方法研究了该模型的检验问题, 最后基于模型(1.1) 分析了波士顿住房数据. You和Zhou[2] 基于经验似然方法研究了模型(1.1)的估计问题, 构造了模型中参数的经验似然置信域.Ahmad等[3]考虑模型 (1.1)在误差是条件同方差时情形, 利用级数方法给出了参数的 有效估计, 并将级数方法推广到误差是条件异方差时,研究了参数的估计. Huang和Zhang[4] 则利用经验似然方法研究了模型 (1.1) 中非参数分量 α(⋅) 的统计推断问题.
对不完全数据的研究是目前统计界的热点问题. 对于缺失数据, 如果仅用可以完全观测到的数据进行统计推断, 得到的估计量往往会产生偏差,而且一般不是渐近有效的. 为了修正缺失数据带来的偏差,常用的方法有参数似然方法[5,6], 借补方法[7,8]以及逆概率加权方法[9,10], 关于这三种方法特点的比较, 可参见Qin等[11]. 对响应变量缺失的情形, Wang等[12]研究了响应变量缺失下半参数回归模型的统计推断. Wang和Rao[8] 考虑了响应变量缺失下线性模型的经验似然推断, 利用线性回归借补方法填充所有的缺失值,从而做进一步的推断. Xue[13]在响应变量缺失下运用经验似然方法构造了线性模型 中参数的经验似然置信域. 赵培信[14] 研究了响应变量缺失下变系数部分线性模型的 经验似然估计. 对协变量缺失的情形,Liang等[15]研究了协变量随机缺 失时部分线性模型的估计问题, 给出了估计量的渐近性质和估计效率. Wang[16] 利用模型纠正方法和一种新的加权方法 给出了部分线性模型中参数和非参数分量的估计, 并证明了它们的渐近性质. 本文基于逆概率加权方法,考虑协变量缺失下变系数部分线性模型的统计推断.
假设 {(Yi,Xi,Zi,Ui),i=1,⋯,n} 是来自模型 (1.1) 的独立同分布观测样本,即 Yi=XTiα(Ui)+ZTiβ+εi, i=1,2,⋯,n,(2.1) 其中 (Yi,Ui,Xi) 可以完全观测,协变量 Zi 随机缺失, 当 δi=1 时 Zi 可观测; 即当 δi=0 时 Zi 缺失. 进一步假定 Zi 为随机缺失 (见文献 [17]),即 P(δi=1|Yi,Ui,Xi,Zi)=P(δi=1|Yi,Ui,Xi)=π(Yi,Ui,Xi).(2.2) (2.2) 式暗含在给定 Yi,Xi 和 Ui 的条件下, δi 和 Zi 是独立的. 该假定是统计分析中常用的假定且在实际应用中是合理的.
\par 在实际应用中,π(Yi,Ui,Xi)往往是未知的, 它是一个需要估计的讨厌参数. 本文考虑使用logistic回归模型作为缺失机制,即 π(Yi,Ui,Xi,θ)=exp(θ0+θ1Yi+θ2Ui+θT3Xi)1+exp(θ0+θ1Yi+θ2Ui+θT3Xi), 其中 θ=(θ0,θ1,θ2,θT3)T 是未知的参数向量. 由文献 [18] 知 θ 的拟似然估计为 ˆθ. 记 Vi=(Yi,Ui,XTi)T, 则 π(Vi,θ) 的估计为 π(Vi,ˆθ). \par 若 β 已知,则模型 (2.1) 可以转变为下面的形式 Yi−ZTiβ=XTiα(Ui)+εi, i=1,2,⋯,n.(2.3) 显然在 β 已知的的条件下,模型 (2.3) 就是一般的变系数模型. 下面采用局部逆概率 加权最小二乘方法来估计系数函数 α(⋅). 在任意给定 u0 的某个小邻域内,对 αj(u) 作~Taylor 展开 αj(u)≈αj(u0)+α′j(u0)(u−u0)≡aj+bj(u−u0), j=1,⋯,p. 记 a=(a1,⋯,ap)T,b=(b1,⋯,bp)T. 对任意给定的 β,极小化下式 n∑i=1{Yi−ZTiβ−XTi[a+b(Ui−u)]}2δiπ(Vi,ˆθ)Kh(Ui−u),(2.4) 可得 a 和 b 的估计,其中 Kh(⋅)=K(⋅/h)/h, K(⋅) 是核函数,h 是带宽. 为表示方便,引入以下记号 Y=(Y1,⋯,Yn)T,X=(X1,⋯,Xn)T,Z=(Z1,⋯,Zn)T,Wu=diag(Kh(U1−u),⋯,Kh(Un−u)),ˆΔ=diag(δ1/ˆπ1,⋯,δn/ˆπn), ˆπi=π(Vi,ˆθ),i=1,⋯,n. 且 Du=(XT1 U1−uhXT1⋮ ⋮XTn Un−uhXTn),M=(XT1α(U1)⋮XTnα(Un)). 记 ˆa,ˆb 是最小化 (2.4) 式的解,则可得 [ˆa,hˆb]T=(DTuWuˆΔDu)−1DTuWuˆΔ(Y−Zβ). 从而可得 α(u) 的估计为 ˜α(u,β)=[Ip 0p](DTuWuˆΔDu)−1DTuWuˆΔ(Y−Zβ).(2.5) 故 M 的估计可以定义为 ˆM=ˆS(Y−Zβ),其中 ˆS=((XT10){DTu1Wu1ˆΔDu1}−1DTu1Wu1ˆΔ(XT20){DTu2Wu2ˆΔDu2}−1DTu2Wu2ˆΔ⋮(XTn0){DTunWunˆΔDun}−1DTunWunˆΔ) \par 下面利用逆概率加权估计方程的思想,给出 β 的估计. 构造如下估计方程 1nn∑i=1δiπ(Vi,ˆθ)~Zi{˜Yi−~ZiTβ}=0,(2.6) 其中 ˜Z=(˜Z1,⋯,˜Zn)T=(I−ˆS)Z ,˜Y=(˜Y1,⋯,˜Yn)T=(I−ˆS)Y. 则由 (2.6) 式得 ˆβ={ZT(I−ˆS)TˆΔ(I−ˆS)Z}−1ZT(I−ˆS)TˆΔ(I−ˆS)Y.(2.7) 将 ˆβ 代入 (2.5) 式可得 α(u) 的最终估计 ˆα(u)=[Ip 0p]{DTuWuˆΔDu}−1DTuWuˆΔ(Y−Zˆβ).(2.8)
为了得到本文的结果,先给出一些正则条件.
C1: U 有有界支撑 U, 它的密度函数 f(⋅) 是二次连续可微的,且 f(⋅)>0.
C2: 矩阵 Γ(U)=E(XXT|U) 和 E(ZZT∣U) 为正定矩阵. Γ(U),Γ−1(U) 和 Φ(U)=E(XZT|U) 均为Lipschitz 连续的.
C3: 存在某个 s>2.5,使得 E‖X‖2s<∞ 且 E‖Z‖2s<∞, 且存在某个 ϵ<2−s−1, 使得 n2ϵ−1h→∞.
C4: {αj(u),j=1,⋯,p} 关于 u 有连续的二阶导数.
C5: 核函数 K(⋅) 是对称概率密度函数, 满足二阶核条件,且有 ∫u4K(u)du<∞.
C6: 带宽 h=Cn−1/5,其中 C>0 为某给定常数.
定理 1在条件C1--C6 成立下, 若 β 是真实参数,则 √n(ˆβ−β)L⟶N(0,Σ1), 其中 "L⟶" 表示依分布收敛,且 Σ1=Ω−1uΛΩ−1u, Ωu=E(ZZT)−E(ΦT(U)Γ−1(U)Φ(U)), Λ=σ2{E[π−1(V)ZZT]−E[π−1(V)ΦT(U)Γ−1(U)Φ(U)]}.
注 1定理 1 表明利用逆概率加权思想得到的~profile 最小二乘估计 ˆβ 具有渐近正态性. 特别地, 当 Z 可以被精确观测时,π(V)=I, 从而 ˆβ 与Fan 和 Huang[1] 提出的估计量具有相同的渐近方差.
定理2在条件C1--C6 成立下, 当 n→∞ 时,有 √nh(ˆα(u)−α(u)−12h2μ2α″(u))L⟶N(0,Σ2), 其中 Σ2=σ2ϑ0f−1(u)Γ−1(U)E(π−1(V)XXT|U)Γ−1(U), ϑ0=∫K2(u)du, μ2=∫u2K(u)du.
虽然由定理 1 可以构造参数 β 的置信域, 但首先需要估计渐近方差,这是比较复杂的. 因此下面考虑经验似然方法 来构造 β 的置信域. 记 ˆR(u)=[Ip 0p]{DTuωuˆΔDu}−1DTuωuˆΔ=(ˆR1(u),⋯,ˆRn(u)), 为构造 β 的经验似然比函数,首先引入一个估计的辅助随机变量 ˆηi(β)=δiπ(Vi,ˆθ)^Zi{^Yi−^ZiTβ},(2.9) 其中 ^Yi=Yi−XTiˆg(Ui), ^Zi=Zi−ˆφ(Ui)TXi, ˆg(u)=n∑i=1ˆRk(u)Yk, ˆφ(u)=n∑i=1ˆRk(u)ZTk. 由 (2.9) 式可知 E{ˆηi(β)}=op(1), 因此关于 β 的估计经验对数似然比函数可定义为 ˆL(β)=−2max{n∑i=1log(npi)|pi≥0,n∑i=1pi=1,n∑i=1piˆηi(β)=0},(2.10) 由 Lagrange 乘子法,ˆL(β) 可以表示为 ˆL(β)=2n∑i=1log(1+λTˆηi(β)),(2.11) 其中 λ 为Lagrange乘子,满足 1nn∑i=1ˆηi(β)1+λTˆηi(β)=0. 极小化 ˆL(β) 可得 β 的估计 ˜β, 称为最大经验似然估计 (MELE). 由MELE ˜β 及 (2.5) 式也可以得到系数函数 α(u) 的估计.
定理 3在条件C1--C6 下, 如果 β 为参数真值,那么有 ˆL(β)L⟶χ2q. 其中"L⟶" 表示依分布收敛, χ2q 为自由度为 q 的卡方分布.
以 χ2q(1−α) 记 χ2q 的 1−α 分位数, 其中 0<α<1. 基于定理 3,可以构造 β 的置信域 Cα(β)={β:ˆL(β)≤χ2q(1−α)}.
本节通过数据模拟研究所提方法在有限样本下的表现.考虑如下半参数变系数部分线性模型 Y=sin(2πU)X+2Z1+1.5Z2+1.5Z3+ε,(3.1) 其中 β=(2,1.5,1.5)τ,U 服从 (0,1) 上的均匀分布,X∼N(1,1),Z1∼N(1,1), Z2∼N(2,1),Z3∼N(0,1),误差 ε∼N(0,0.52). 在整个模拟研究中,我们选用~Epanechnikov 核函数, 即 K(u)=0.75(1−u2)+,并采用交叉核实法来选取最优带宽 hCV, 使其满足下式达到最小 CV(h)=1nn∑i=1δiπi{Yi−XTiˆα[i](Ui)−ZTiˆβ[i]}2, 其中 ˆα[i](Ui) 和 ˆβ[i] 分别为去掉第 i 个观测值后的 α(u) 和 β 的估计. 考虑如下三种缺失概率
(1)~ π1(y,u,x)=P(δ=1∣Y=y,X=x,U=u)={1+0.02exp(y−x−u−3.5)}−1;
(2)~ π2(y,u,x)=P(δ=1∣Y=y,X=x,U=u)={1+0.08exp(y−x−u−3)}−1;
(3)~ π3(y,u,x)=P(δ=1∣Y=y,X=x,U=u)={1.25+exp(y−x−u−6)}−1.
以上三种情况对应的数据 Z 的平均缺失概率分别为0.10,0.25 和0.40. 对于参数分量,本节模拟研究比较了两种方法: 本文提出的逆概率加权方法 (ANA) 和 忽略掉缺失值的估计方法 (NNA). 在模拟过程中,样本容量分别取 n=50,100,150 和~200,重复模拟 1000 次. 在三种不同的缺失水平下,给出了两种估计量的偏差、标准差和均方误差. 模拟结果见表 1 和表 2. 另外,在样本量 n=200 和协变量 Z 的平均缺失概率为0.25 下, 我们在图 1 中给出了系数函数 α(⋅) 的估计曲线.
从表 1 和表 2 中可以得到以下结论: (1) 对给定的缺失概率, 随着样本量的增加,两种估计方法给出的估计量的偏差和均方误差都在减小; (2) 对给定样本,随着缺失比例的增加, 两种估计方法给出的估计量的标准差和均方误差均增加; (3) 对给定缺失概率和样本量的情况下, 基于~ANA 方法所得 β 估计的绝对偏差,标准差和均方误差均小于~NNA 方法, 这说明本文所提估计方法是有效的, 原因在于本文的方法考虑了缺失值的信息.
另外,对比图 1 中估计曲线和真实曲线,可以看出本文所提方法是理想的.
下面利用经验似然方法考虑模型 (3.1) 中参数 β 的置信域精度. 在三种不同的缺失水平下,利用本文提出的经验似然 (IEI) 和正态 逼近 (INA) 两种方法比较了 β 的95\% 置信区间的平均长度以及对应的覆盖概率. 在模拟过程中,样本容量分别取 n=100,150 和200. 模拟结果在表 3 中给出.
从表 3 中可以看出: (1) IEL 明显优于INA. IEL 比INA 给出了更短的置信区间和更大的覆盖概率, 在缺失概率较小和样本量较大情况下,IEL 给出的覆盖概率接近95\%; (2) 对于给定的缺失概率,随着样本量的增加, 经验似然和正态逼近方法得到的置信区间长度均缩短. 相反,对固定的样本量, 随着缺失概率的增加,经验似然和正态逼近给出的置信区间长度均增加. 因此, 缺失概率的大小会影响 置信域的精度.
下面我们应用本文提出的方法分析来自某AIDS 病研究中心的一组数据. 数据 集包括从 1984 年到 1991 年期间感染 HIV 病毒的283 个同性恋者的 HIV 状况. 该数据曾被许多学者利用变系数模型来进行研究 (例如,Huang 等[19], Xue 和Zhu[20]等). 该数据集的目标是研究患者的 CD4 浓度随时间的变化趋势, 以及患者的抽烟状态,HIV 感染前的 CD4 浓度, 和 HIV 感染时患者的年龄对 HIV 感染后的患者的 CD4 浓度的影响. Huang 等[19] 研究表明,在显著水平为0.05 的情况下, HIV 感染前 CD4 浓度对 HIV 感染后 CD4 浓度的影响为常数效应,不随时间变化, 并且患者感染 HIV 时的年龄以及患者的 抽烟状态对 HIV 感染后 CD4 的浓度均没有显著影响. 该数据是纵向数据, 在此我们不考虑其相关结构. 因此,可以建立如下模型 Y(t)=α0(t)+α1(t)X+β1Z1+β2Z2+ε(t), 其中 Z1 为 HIV 感染前标准化的 CD4 浓度, Z2 为患者的抽烟状态 (1 = 抽烟,0 = 不抽烟), Y(t) 为时刻 t 患者的 CD4 浓度,X 为 HIV 感染时患者的年龄, α0(t) 为 CD4 浓度基准函数,表示 HIV 感染 t 年后, 患者的平均 CD4 浓度.
为了应用该数据说明本文的方法,假定20\% 的 Z1 和 Z2 值缺失, 模拟中 δ 的产生是随机的,核函数取为 K(t)=0.75(1−t2)+,并通过交叉核实方法选择带宽 hCV=0.405. 利用本文提出的EL 方法所得的参数分量 (β1,β2) 的95\% 置信域见图 2.
从图 2 中可以看出 β2 的 95% 置信区间覆盖零点, 且其中 ‘‘∘" 表示 β 的MELE,˜β=(2.993,0.109). 以上说明患者的抽烟状态 (Z2) 对 HIV 感染后的患者的 CD4 浓度没有显著影响,这与 Huang 等 [19] 研究所得的结果基本一致.
记 cn=h2+{log(1/h)nh}1/2,μk=∫ukK(u)du,ϑk=∫ukK2(u)du, k=0,1,2,4. τi=(1,Yi,Ui,XTi)T, λmin 表示 n∑i=1τiτTi 的最小特征根.
引理 1 假设 (X1,Y1),⋯,(Xn,Yn) 是 n 个独立同分布的随机向量, 其中 Yi 是一维随机变量. 若进 一步假设 E|Y1|s<∞ 且 supx∫|y|sf(x,y)dy<∞, 其中 f(⋅,⋅) 表示 (X,Y) 的联合概率密度函数. 令 K(⋅) 表示满足Lipschitz 条件 且有有界支撑的有界正函数. 那么给定某个 ε<1−s−1, 使得 n2ε−1h→∞ 时,则有 supx|1nn∑i=1{Kh(Xi−x)Yi−E[Kh(Xi−x)Yi]}|=Op({log(1/h)nh}1/2).(5.1)
证由Mack 和Silverman[21] 可以得到本引理的证明.
引理2在条件 supi≥1‖τi‖<∞ 且 λmin→∞ 下, θ 的拟似然估计 ˆθ=(ˆθ0,ˆθ1,ˆθ2,ˆθT3)T, 满足 √n(ˆθ−θ)=A−1n−1/2n∑i=1τi(δi−πi)+op(1),(5.2) 其中 A=E(τ1τT1π1(1−π1)).
证本引理的证明可参见文献 [18].
引理 3假设条件C1--C6 成立, 则在 U∈U 上,当 n→∞ 时,有 DTuWuΔDu=nf(U)Γ(U)⊗(1 00 μ2){1+Op(cn)},(5.3) DTuWuΔZ=nf(U)Φ(U)⊗(1,0)T{1+Op(cn)},(5.4) DTuWuΔM=nf(U)Γ(U)⊗(1,0)Tα(U){1+Op(cn)},(5.5) DTuWuΔ(XT1(U1−uh)2⋯XTn(Un−uh)2)=nμ2f(U)Γ(U)⊗(1,0)T{1+Op(cn)}.(5.6) 证由于 (5.3)--(5.6)式的证明过程类似. 这里我们仅给出 (5.3) 式的证明. 经过一系列简单代数运算,可以得到 DTuWuΔDu=(laystyle∑ni=1laystyleδiπiXiXTiKh(Ui−u) laystyle∑ni=1laystyleδiπiXiXTi(Ui−uh)Kh(Ui−u)laystyle∑ni=1laystyleδiπiXiXTi(Ui−uh)Kh(Ui−u) laystyle∑ni=1laystyleδiπiXiXTi(Ui−uh)2Kh(Ui−u)), 由引理 1 及条件期望的平滑性,易证 (5.3) 式成立.
引理4在条件 C1--C6 成立下,有 n−1ZT(I−S)TΔ(I−S)ZP⟶Ωu,(5.7) 其中 Ωu 在定理 1 中定义.
证结合引理 1--3, 类似Fan 和 Huang[1] 中引理2 的证明可证本引理成立.
引理5在条件C1--C6 成立下,有 n−1ZT(I−S)TΔ(I−S)M=Op(c2n).(5.8)
{\heiti证结合引理 1--3, 类似Fan 和 Huang[1] 中引理4 的证明可证本引理成立.
引理6在条件 C1--C6 成立下,有 1n(DTuWuˆΔDu−DTuWuΔDu)=op(n−12+12s).(5.9)
证简单计算可得 1nDTuWu(ˆΔ−Δ)Du=(laystyle1nn∑i=1ψiXiXTiKh(Ui−u) laystyle1nn∑i=1ψiXiXTi(Ui−uh)Kh(Ui−u)laystyle1nn∑i=1ψiXiXTi(Ui−uh)Kh(Ui−u) laystyle1nn∑i=1ψiXiXTi(Ui−uh)2Kh(Ui−u)), 其中 ψi=δiπ(Vi,ˆθ)−δiπ(Vi,θ). 由引理2 知 ˆθ−θ=Op(n−1/2). 将 1π(Vi,ˆθ) 在点 θ 处 Taylor 展开得 δiπ(Vi,ˆθ)−δiπ(Vi,θ)=−δiπ′(Vi,θ)π2(Vi,θ)(ˆθ−θ)+op(n−1/2)=−δi(1−π(Vi,θ))π(Vi,θ)τTi(ˆθ−θ)+op(n−1/2)=op(n−12+12s),(5.10) 其中 π′(Vi,θ)=∂π(Vi,θ)/∂θ. 则由 (5.10) 式可得 1nn∑i=1[δiπi(ˆθ)−δiπi(θ)]XiXTiKh(Ui−u)=1nn∑i=1−δi(1−πi(θ))πi(θ)τTi(ˆθ−θ)XiXTiKh(Ui−u)+op(n−1/2)=op(n−12+12s), 从而可得 1n(DTuWuˆΔDu−DTuWuΔDu)=op(n−12+12s).
引理7在条件 C1--C6 成立下,有 ˆS−S=op(n−12+12s).(5.11)
证类似引理6 的证明可证明本引理,此处省略.
引理8在条件 C1--C6 成立下,有 sup0<u<1‖ˆφ(u)−Γ−1(u)Φ(u)‖=Op(cn),(5.12) sup0<u<1‖ˆg(u)−Γ−1(u)Φ(u)β−α(u)‖=Op(cn),(5.13) 其中 ˆφ(u) 和 ˆg(u) 由 (2.9) 式定义,并且 cn=h2+(log(1/h)nh)1/2.
证类似引理6 的证明可得 DTuWu(ˆΔ−Δ)Z=op(n−(12−1s)). 于是 ˆφ(u)=[Ip 0p]{DTuWuˆΔDu}−1DTuWuˆΔZ=[Ip 0p]{DTuWuΔDu}−1DTuWuΔZ+op(1). 因此由引理 3,对 u∈(0,1),一致有 ˆφ(u)=Γ−1(u)Φ(u){1+Op(cn)}. 因此 (5.12) 式成立. 利用类似方法可以证明 (5.13) 式, 这就完成了本引理的证明.
引理9在条件 C1--C6 成立下, 如果 β 是参数真值,则有 1√nn∑i=1ˆηi(β)L⟶N(0,Ω),(5.14) 其中 Ω=σ2{E[π(V)−1ZZT]−E[π(V)−1ΦT(U)Γ−1(U)Φ(U)]}.
证简单计算可得 1√nn∑i=1ˆηi(β)=1√nn∑i=1δiπ(Vi,ˆθ){Zi−φT(Ui)Xi}εi+1√nn∑i=1δiπ(Vi,ˆθ)(φ(Ui)−ˆφ(Ui))TXiεi+1√nn∑i=1δiπ(Vi,ˆθ){Zi−φT(Ui)Xi}XTi{α(Ui)−ˆg(Ui)+ˆφ(Ui)β}+1√nn∑i=1δiπ(Vi,ˆθ)(φ(Ui)−ˆφ(Ui))TXiXTi{α(Ui)−ˆg(Ui)+ˆφ(Ui)β}=:A1+A2+A3+A4. 将 1π(Vi,ˆθ) 在 θ 点 Taylor展开,可得 A1=1√nn∑i=1δiπ(Vi,θ){Zi−φT(Ui)Xi}εi−1√nn∑i=1δiπ′(Vi,θ)π2(Vi,θ)(ˆθ−θ){Zi−φT(Ui)Xi}εi+op(1)=1√nn∑i=1Ii1+1√nn∑i=1Ii2+op(1). 注意到 E{Ii1}=0, Var{Ii1}=Ω, 故由中心极限定理知 1√nn∑i=1Ii1L⟶N(0,Ω).
下证 1√nn∑i=1Ii2P⟶0. 记 ai=−δi(1−π(Vi,θ))π(Vi,θ){Zi−φT(Ui)Xi}εi, 则 Ii2=aiτTi(ˆθ−θ). 设 Ii2,s 为 Ii2 的第 s 个分量,ai,s 为 ai 的第 s 个分量. 由 √n(ˆθ−θ)=Op(1), 结合引理8 以及文献 [14,引理2.6.1,2.6.2]可得 |1√nn∑i=1Ii2,s|=1√n|n∑i=1τTi(ˆθ−θ)ai,s|≤1√nmax1≤i≤n|n∑i=1ai,s|max1≤i≤nτTi(ˆθ−θ)=1√nop(n−1/2)Op(√nlogn)=op(1), 因此有 A1L⟶N(0,Ω).
进一步,由 E(Xiεi)=0, E{[Zi−φT(ui)Xi]XTi}=0, 类似引理8,可证 AvP⟶0,v=2,3. 又由文献 [14,引理2.6.1,2.6.2]可得 ‖A4‖≤Op(√nc2n)=op(1). 综合上述,本引理得证.
引理 10在条件 C1--C6 成立下, 如果 β 是参数真值,则有 1nn∑i=1ˆηi(β)ˆηTi(β)P⟶Ω.(5.15)
证证明过程类似文献 [14,引理4.5.3]的证明,此处省略.
定理 1 的证明简单计算可得 √n(ˆβ−β)=√n{ZT(I−ˆS)TˆΔ(I−ˆS)Z}−1ZT(I−ˆS)TˆΔ(I−ˆS)(M+ε). 首先考虑1nZT(I−ˆS)TˆΔ(I−ˆS)Z=1nZT(I−S)TΔ(I−S)Z+1nZT(I−S)TΔ(S−ˆS)Z+1nZT(I−S)T(ˆΔ−Δ)(I−S)Z+1nZT(I−S)T(ˆΔ−Δ)(S−ˆS)Z+1nZT(S−ˆS)TΔ(I−S)Z+1nZT(S−ˆS)TΔ(S−ˆS)Z+1nZT(S−ˆS)T(ˆΔ−Δ)(I−S)Z+1nZT(S−ˆS)T(ˆΔ−Δ)(S−ˆS)Z=:T1+T2+T3+T4+T5+T6+T7+T8. 下面证明 Ti=op(n−1/2),i=2,⋯,8. 由引理7 得 ‖T2‖=1n‖ZT(I−S)TΔ(S−ˆS)Z‖≤n−1⋅n1s⋅n−(12−1s)⋅n12s=op(n−1/2). 结合 (5.10) 式可证 ‖T3‖=1n‖ZT(I−S)T(ˆΔ−Δ)(I−S)Z‖≤n−1⋅n1s⋅n−(12−12s)⋅n1s=op(n−1/2). 类似可证 Ti=op(n−1/2),i=4,⋯,8. 又由引理4 知 T1P⟶Ωu,从而可得 laystyle1nZT(I−ˆS)TˆΔ(I−ˆS)ZP⟶Ωu.
其次考虑1√nZT(I−ˆS)TˆΔ(I−ˆS)(M+ε)=1√nZT(I−S)TΔ(I−S)(M+ε)+1√nZT(I−S)TΔ(S−ˆS)(M+ε)+1√nZT(I−S)T(ˆΔ−Δ)(I−S)(M+ε)+1√nZT(I−S)T(ˆΔ−Δ)(S−ˆS)(M+ε)+1√nZT(S−ˆS)TΔ(I−S)(M+ε)+1√nZT(S−ˆS)TΔ(S−ˆS)(M+ε)+1√nZT(S−ˆS)T(ˆΔ−Δ)(I−S)(M+ε)+1√nZT(S−ˆS)T(ˆΔ−Δ)(S−ˆS)(M+ε)=:R1+R2+R3+R4+R5+R6+R7+R8.
下面证明 Ri=op(1),i=2⋯,8. 对于R2=1√nZT(I−S)TΔ(S−ˆS)M+1√nZT(I−S)TΔ(S−ˆS)ε=:R21+R22. 由条件 C3 及引理~7 知 ‖R21‖≤n−1/2⋅n1s⋅n−(12−1s)⋅n12s=op(1). 类似可以证明 ‖R22‖≤n−1/2⋅n1s⋅n−(12−1s)⋅n12s=op(1) 及 Ri=op(1),i=3,⋯,8. 由此结合引理4 和引理~5 可得√n(ˆβ−β)=n−1/2Ω−1un∑i=1[Zi−ΦT(Ui)Γ−1(Ui)Xi]δiπiεi{1+op(1)}. 由中心极限定理和 Slutsky 定理可得 √n(ˆβ−β)L⟶N(0,Σ1). 定理得证.
定理2 的证明简单计算可得 ˆα(u)=[Ip 0p]{DTuWuΔDu}−1DTuWuΔ(M+ε+Z(β−ˆβ))+[Ip 0p]{DTuWuΔDu}−1DTuWu(ˆΔ−Δ)M+[Ip 0p]{DTuWuΔDu}−1DTuWu(ˆΔ−Δ)ε+[Ip 0p]{DTuWuΔDu}−1DTuWu(ˆΔ−Δ)Z(β−ˆβ)=:S1+S2+S3+S4.
首先考虑 S2,由引理~2 和引理 3 可得 S2=[Ip 0p]{DTuWuΔDu}−1DTuWu(ˆΔ−Δ)αT(U)X=f−1(U)Γ−1(U)1nn∑i=1−δi(1−πi(θ))πi(θ)τTi(ˆθ−θ)XiXTiαT(U)Kh(Ui−u)+op(n−1/2)=op(n−12+12s). 同理可证 S3=op(n−1/2), S4=op(n−1). 下面考虑 S1, S1=[Ip 0p]{DTuWuΔDu}−1DTuWuΔM+[Ip 0p]{DTuWuΔDu}−1DTuWuΔε+[Ip 0p]{DTuWuΔDu}−1DTuWuΔZ(β−ˆβ)=:S11+S12+S13. 由 α(Ui) 在点 u 处的 Taylor 展开, 结合引理 3,简单计算可得 S11=α(u)+h22μ2α″(u)+op(h2), √nh S12=√nh[Ip 0p]{DTuWuΔDu}−1DTuWuΔε=f−1(U)Γ−1(U)(√nh1nn∑i=1δiπiXiKh(Ui−u)εi)+op(1). 由中心极限定理和 Slutsky 定理可得 √nhS12L⟶N(0,Σ2).
进一步,由定理 1 和引理 3 可得 √nhS13=√nhΓ−1(U)Φ(U)(ˆβ−β)=op(1).
综上所述,定理2 得证.
定理 3 的证明由引理9,经简单计算可得 max1≤i≤n‖ˆηi(β)‖=op(n−1/2), λ=Op(n−1/2).
对 (2.11) 式,利用Taylor 公式可得 ˆL(β)=2n∑i=1{λTˆηi(β)−12(λTˆηi(β))2}+op(1).(5.16) 进而可得 0=1nn∑i=1ˆηi(β)1+λTˆηi(β)=1nn∑i=1ˆηi(β)−1nn∑i=1ˆηi(β)ˆηTi(β)λ+1nn∑i=1ˆηi(β){λTˆηi(β)}21+λTˆηi(β). 由引理9,10 得 n∑i=1[λTˆηi(β)]2=n∑i=1λTˆηi(β)+op(1).(5.17) λ={n∑i=1ˆηi(β)ˆηTi(β)}−1n∑i=1ˆηi(β)+op(n−1/2).(5.18) 结合 (5.16)--(5.18) 式,简单计算可得 ˆL(β)={1√nn∑i=1ˆηTi(β)}TΩ−1{1√nn∑i=1ˆηTi(β)}+op(1), 因此由引理9,10,可知定理成立.