第25卷第1期 湖南工业大学学报 VO1.25 No.1 2011年1月 Journal of Hunan University of Technology Jan.2O11 P~混合误差下半变系数模型的相合估计 李强,陈志彬 (湖南工业大学理学院,湖南株洲412007) 摘 要:讨论了误差序列是 混合序列的变系数模型。给出变系数模型系数函数的局部多项式估计,并证 明了该模型系数函数估计是弱相合估计。 关键词:半变系数模型; 混合;局部多项式估计;弱相合 中图分类号:0212.7 文献标志码:A 文章编号:1673-9833(2011)01 ̄50-03 Consistent Estimation of Semi—Varying—Coefifcient Models under Jp—Mixing Sequences Li Qiang,Chen Zhibin (School of Science,Hunan University ofTechnology,Zhuzhou Hunan 412007,China) Abstract:Discusses the semi—varying—coefifcient models whose error sequence is ̄-imxing.Gives hte local polyomila estimation of coeficient functions of semi—varying—coefifcient models.Proves that the estimation of hte coeficient functions of this model is weak consistency. Keywords:semi—varying—coefifcient models; 一mixing;local polynomial method;weak consistency . T.J.Hastie和R.J.Tibshirani口 于1993年提出了变系 1 P混合误差的定义 数模型。J.Fan等人心 在检验函数系数是否真的变化 时,提出了半变系数模型,其定义如下: 设{ I i=1,2,…, ,…j是概率空间( , ,P)上的 】,=∑ ( ) +∑卢 Zj+s, (I) 随机变量序列, = ( li∈ScN)为 一域,在 中给 i=I i=I 定 一域 ,R,令 式中:】,为响应变量;U∈R,随机误差s独立于协变 JD ,R)=sup{lcor(xr)llx L2(F),YE )), ( , , 2,・一, ,z ,z2,’ 。,z ),¨.£∽)=0,D( = 仃 ; .(・), 1,2,…,P乃术 数系数; .乃常数 其中,co0r ,( = =_ 为相关系数。 半变系数模型是一类应用较广泛的模型,例如, 当P=1,Xi 1时,式(1)就成为部分线性模型。文 Bradley 引入如下相依系数:对k≥0令 献【3—7】等都对此进行了讨论。吕士钦等人研究了误差 (Ij})=sup{p( ̄s,簖)I有限子集S,TcN,H.dist(S,T)≥七), 独立时,半变系数模型PLS估计(profile least squares 显然O ( 1) ( 1,且 (0)=l。 estimation)的一致收敛速度盯 ,本文在此基础上讨论J9 定义1 对随机序列{ Ii=1,2,…,”,…J,如果存 混合误差条件下半变系数模型估计的相合性。 在k∈N使 ( ) 1,则称{ }是 混合序列。 收稿日期:2010-12-09 基金项目:湖南省教育厅科学研究基金资助项目(10C0656),湖南工业大学教学改革基金资助项目(09E66) 作者简介:李强(1975一),男,湖南湘潭人,湖南工业大学讲师,博士,主要从事非参数估计和高维数据分析方面的研究, E—mail:lqdr2008@126.com 第1期 李强,等 混合误差下半变系数模型的相合估计 5 1 独立序列是 混合序列的特殊情形。p混合序列要 求对某距离为k的任意有限子集之间的相关度不大, 这是一种很弱的相依性要求,因此P混合序列是一类 应用广泛且很有研究价值的随机变量序列。自文献【8】 引入P混合概念以来,它引起了许多学者的兴趣。近 年来, 混合序列极限理论的研究已取得不少成果 。 从而得: )= ,{o:w.o3~ (】,一 )。 ~ 3 主要假设 ; 首先给出证明需要的一些假设: Al 函数系数 (・)( l,2,…,p)有连续的二阶导数; A2 随机变量 有有界支撑 ,其密度函数p(・) 2系数的PLS估计 假设 {( , , ,…, , , ,…,Zt., ),k=1 2一, } 是模型(1)的随机样本,对任意给定的 ,模型(1) 可记作: =∑ ( ) + ,k=l,2,…, 。 (2) 其中, = - ̄,8sZk:。若 )( 1,2,…,p)有连续的二 阶导数,则在 。附近可记为 ,(“)≈ (“0)+ ( o)(“一glo)三a,+6『( 一go).1-l,2,…,P。 根据局部线性方法的思想,求解 r ]2 然 【 一善 j (G-Uo), 其中, : + (UI一 );K (’)是核函数,h=h 是正 数序列,称为窗宽,且当n一∞时,h一0,nh—o。。 令 =( , ,…, ) ,£= 。,£ ,…, ) , =( ,Xz,…,Xn)r, =( , ,…, ) ,Zi=(Zf。,Zi:,…, )。, =diag(K (U.一“),K (U 一”),…,K ( 一“)), y(甜):al( ),口 ( ),…, (”), 6l(“),hb2(u),…,hbp(”)) 『.aT(u.)x。 X x 且 :I【 -( ) , = ; x TU1--U A T 则式(2)可记作:y- =M+z, (4) 可得式(3)的解为: (“)=[ T D )~o:w(r一 )。 M的估计为: f o]{ . 。 D TWu (Y.Zf1)I M=I ; I(y—z )垒 l ix:ol[o ̄ w. )一 T (r-213)I S(Y一 ), (5) 将 入式(4)得:(,一s)r=(,一 ) + , (6) 对线性模型(6)应用最小二乘法,可得: :f (J— ) (,一s)2 } (J— ) (,~S)Y, 于是 (“)可修正为: (“)=( T D ) (Y-Z/ ̄), Lipschitz连续; A3 核函数 (・)是一有界对称的,且为有界紧支 撑的概率密度函数; A4 E(XX IU=“)对任意的“∈Q非奇异, E( =“),E(xxTIu:“)  ̄E(XZ Iv---U)都是 Lipschitz连续的; A5存在s>2, ̄-EIIXII <" ̄fllzll <o。; A6当n oo时,h O,nh.÷o。,h=n~,0<l<112; A7 设(sf1i=l,2,一・, )为 混合同分布序列, E(s )=0,E(e )= 。, ( 一 )( 】 C,c为常数。 4结论及证明 引理1【l2 设{ I i=1,2,…, )为 混合随机变量 序列,E( ):。,E(1 I吾)<oo,。<a<l,I l≤ 坨一, 0<6≤2/ ,∑ 2f<cn一,0>0,c为常数,则 Zan 0。 证明参见文献【l2l中定理2的证明。 引理2在假设A1~A6成立的条件下,有 一 山0。 证明 由文献[13】中定理4.1易得。 引理3在假设A1~A7成立的条件下,有 C以山0, 其中, , = 一∑( ( 一 )( + f ̄= ’zo.( ,))(\ , J 。 证明令 = , : +c ,其中 。喜f_l l,=l (哪)(、 竿n “ , I n 喜 (l )(、 , 卜 。的证明见文献【7】,下面只需证 52 湖南工业大学学报 2011正 一。。-@v i=n-JXiKh(Ui ( =∑ s,,在引理1中,令 :1, 0< =1/2-I<1/2<2/ =2,由假设A6和A7可知 、f, n2)=∑v.ie 满足引理1的条件,故有: fzl n =∑ si=l 0。从而, 2)=∑ £i=1 — 0。 同理可证,C n2。)=∑ fs,— 0。 i=l 引理4在假设A1~A6成立的条件下,有 一—— — 0,A . p( )厂( ) —— — 0, 其中, 一 了XiKh(Ui"ld j=l( ( ) ( (州 ) ]J , n 。 F(u)=E(XX = ), J (“) 。 证明类似文献【13】中引理7.2的证明。 定理1在假设A1~A7成立的条件下, (“)一 (“)——巴 0, 从而&f(甜)一 f(“)—— 0。i=1,2…,P。 证明 ( )一 (材): { T D ) (y一 一 ’,)= { ( 一 y)+( 一 + )= (A.,o-4.3)- 斟 由引理2—4可知,A 有界,且 _ 0 =D ( )—— 0, =0,1。定理得证。 . 参考文献: 【1】 Hastie T J.Tibshirani R J.Varying—Coefifcient Models[J].J. Roy.Statist.Soc.Ser.B,1993.55(4):757—796. [2】 Fan J,Zhang W Y.Simultaneous Confidence Bands Are nad Hypothesis Testing in Varying Coeffcient Models[J]. Scandinavian Journal of Statistics,2ooO,27(4):715—731. 【3】 李志强,薛留根.缺失数据下的半参数变系数模型的借补 估计[J].应用数学学报,2009,32(3):422-430. Li Zhiqiang,Xue Liugen.The Imputation Estimators of Semiparametric Varying-Coefifcient Models with Missing Data [J1.Acta Mathematicae Applicatae Sinica,2009,32(3): 422-430. [4]吕士钦,张日权,卢准炜.删失数据下半变系数模型估计 的渐近正态性[J】.应用概率统计,2009,25(6):641-648. Lv Shiqin,Zhang Riquan,Lu Zhunwei.The Asymptotic Normality of Estimation on Semivarying Coefficient Models wiht Censored Data[J].Chinese Journal ofApplied Probability nad Statistics,2009,25(6):641—648. [5】Lv Shiqin,Zhang Riquan,Lu Zhunwei.Asymptotic Normality of Constrained Profile Least Squares Estimation Oil Semivarying Coeffcient Models[J].Chinese Journal of Engineering Mathematics,2009,26(2):315—319. 【6】 富伯亭,郑刘锁.半变系数模型的M估计【J].工程数学学 报,2009,26(4):609—617. Fu Boting,Zheng Liusuo.M—Estimation on Semi—Varying Coefficient Models[J].Chinese Journal of Engineering Mathematics,2009,26(4):609—617. [7]吕士钦,张日权,蒙家福.半变系数模型PLS估计的一致 收敛速度【J1.太原科技大学学报,2006,27(6):432-434. Lv Shiqin,Zhang Riquan,Meng Jiafu.The Uniform Convergence Rates of Profile Least—Squares Estimation on Semivarying Coefficient Models[J].Journal of Taiyuan University of Science and Technology,2006,27(6):432- 434. [8 Br8】adleyRC.EqnivalentMixingCondiitonsforRandomFields 【R】.Chapel Hill:Technical Report,1990. 【9】吴群英.混合序列的概率极限理论【M】.北京:科学出版 社,2006:206-247. WuQunying.LimitTheoryofMixingSequences[M].Beijing: Science Press,2006:206—247. [1O】蔡光辉,郭宝才. 混合序列的对数律和强大数律[J】.科技 通报,2006,22(3):286—289. Cai Guanghui,Guo Baocai.Law of Logarithm and Strong Law ofLarge Numbers ofr P Mixing Sequences[J].Bulletin of Science nad Technology,2006,22(3):286—289. [11】LuZhaoyang,XuWei,NiuYujun.MarcinkiewiczStrongLaws for Weighted Sums of P—Mixing Random Variable Sequences 【J】.Chinese Journal of Engineering Mathematics,2009,26 (2):320-330. 【l2】吴群英.混合序列加权和的完全收敛性和强收敛【J].应用 数学,2002,15(1):1-4. Wu Qunying.Convergence for Weighted Surns of Mixing Random Sequences[J].Mathematica Applicata,2002,15 (1):卜4. 【13】Fan J Q,Huang T.Profile Likelihood Inferences on Semi— parametric Varying—Coefifcient Partilaly Linear Models[J]. Bernoulli,2005,1 1(6):1031—1057. (责任编辑:廖友媛)