函数系数协整模型局部线性估计方法的改进.pdf
《函数系数协整模型局部线性估计方法的改进.pdf》由会员分享,可在线阅读,更多相关《函数系数协整模型局部线性估计方法的改进.pdf(6页珍藏版)》请在咨信网上搜索。
1、理 论 探 讨统计与决策2023年第21期总第633期0引言Cleveland等(1992)1率先研究了函数系数模型,可灵活捕获丰富的动态特征,同时避免完全非参数形式中的“维数困难”现象。随后该模型被广泛用于协整关系的建模,统称为函数系数协整模型。20世纪初,该模型的研究重点是系数的逐点估计方法。方法多基于最小二乘估计构建,此类方法在残差服从正态分布时存在一些良好的性质2。同时,Cai等(2000)3使用两步法对系数进行了局部多项式估计;Sun和Li(2011)4研究了局部常数和局部线性核估计量,且后者具有更小的均方误差。随后学者们关注到数据中的异常值,研究发现基于最小绝对离差估计准则(即L1
2、损失函数)构建的局部估计量效果更好5,6。同样,模型中同方差的设定在解决现实问题时存在很大局限。部分学者假设动态波动中的主要特征由时变因素引起。Xu和Phillips(2008)7、Tu和Wang(2020)8分别在研究自回归模型和函数系数协整模型时均提出构建自适应估计流程来减少时变波动方差对估计量产生的偏差;Zhu(2018)9扩展研究了含厚尾特征和时变波动的自回归模型,提出自适应最小绝对离差估计量(AdaptiveLeast Absolute Deviations Estimator,ALADE)。本文将关注时间序列含厚尾特征及时变波动方差时函数系数协整模型估计流程的优化问题。1局部线性A
3、LADE1.1引入模型不失一般性,构建含厚尾特征和时变波动方差的函数系数协整模型(FCCM):yt=a0(zt)+a1(zt)x1t+a2(zt)x2t+t,x1t=x1t-1+t,x2t=x2t-1+t,t=tt=g(tn)t(1)其中,xit是单位根过程,其系数依赖于单一外生变量zt变化;t、t是服从标准正态分布的新息项;t可服从厚尾分布以体现经济变量中常出现的极端事件;t是用时间的连续或间断函数表示的时变波动因子。该模型设定了易扩展至高维模型和多平滑变量的情景。1.2时变波动方差的处理Tu和Wang(2020)8将自适应估计用于含时变波动方差的函数系数协整模型,形如:yt=(zt)xt+
4、ut,ut=tet(2)提出加权核最小二乘(Weighted KLS)估计,当权重取1 t时估计量为:(z0)=t=1T-2txtxtKh(zt-z0)-1t=1T-2txtytKh(zt-z0)(3)其拥有最小方差-协方差矩阵。其中,K()=k(/h)为核函数,h是窗宽。当使用2的非参数估计值代替权重时,可得自适应核权最小二乘(Adaptive KLS)估计。1.3厚尾特征的处理Tang和Wang(2005)5研究了含异常值或厚尾特征的函数系数模型:yt=j=1pj(zt)xtj+et(4)对于任意z001,在求解系数j(z0)局部点的估计值时,周边点的系数均用线性函数近似:aj+bj(z-
5、z0)。基于绝对离差之和最小,即最小化式(5)可求系数的逐点估计值:minabt=1T|yt-j=1pajxtj-j=1pbj(zt-z0)xtjKh(zt-z0)(5)函数系数协整模型局部线性估计方法的改进曹晓舟(济南大学 数学科学学院,济南 250002)摘要:函数系数协整模型可以克服非参数建模时的“维数困难”,同时体现系数的动态变化,被广泛应用于非平稳数据间复杂问题的研究。实践中经济序列常表现出时变波动方差和厚尾特征,传统核权最小二乘估计方法将不再适用。鉴于此,文章基于L1损失函数重构估计流程,选择表现稳健的局部线性核估计方法,并引入自适应方法,提出局部线性自适应最小绝对离差估计(ALA
6、DE)。模拟结果验证了所提估计方法可提升系数估计精度,优化模型整体拟合效果,同时绝对值交叉验证方法在选取最优窗宽时优势明显。实证分析发现,所提方法可识别中英两国汇率和价差间的动态协整关系,拟合系数平滑且接近理论值。关键词:函数系数协整模型;局部线性ALADE;时变波动方差;厚尾特征中图分类号:O212.7文献标识码:A文章编号:1002-6487(2023)21-0023-06作者简介:曹晓舟(1988),女,山东济南人,博士,讲师,研究方向:数量经济学技术与应用。DOI:10.13546/ki.tjyjc.2023.21.00423理 论 探 讨统计与决策2023年第21期总第633期上述最
7、优化问题可转化为求解线性规划问题,系数的估计值aj(zt)即局部线性最小绝对离差估计。估计过程中窗宽的选择至关重要,建议采用 ACV 方法进行选择:hacv=argminACV(h)=1ht=1T|yt-j=1pa-tj(zt)xtj(6)其中,a-tj(zt)是系数的留一估计值。1.4改进的估计流程基于模型(1)构建局部线性ALADE流程。将L1最优化准则引入估计全流程,选择局部线性估计方法,最优窗宽选择采用绝对值交叉验证方法,考虑可行的自适应估计流程以减少时变波动方差的影响。估计流程见图1。mint-1T|yt-a0-b0(zt-z0)-j=12ajxtj-j=12bj(zt-z0)xtj
8、Kh1(zt-z0)h1=arg minACV(h)=1nt=1T|yt-a-t0-j=12a-tj|xtjgt=i=1nKti|i|i=1nKtiACV(b)=1nt=1T|i|-gt|mint=1T|yt-a0-b0(zt-z0)-j=12ajxtj-j=12bj(zt-z0)xtj|gtKh2(zt-z0)h2=arg minACV(h)=1nt=1T|yt-a-t0-j=12a-tjxtj|gt局部线性LADE并提取残差i窗宽h1的ACV选择时变波动方差的非参估计估计权重时窗宽b的ACV选择局部线性ALADE窗宽h2的ACV选择图1FCCM的局部线性ALADE流程第一,模型的局部线性L
9、ADE。将某点的系数值近似为:a0(zt)a0+b0(zt-z0),aj(zt)xtjajxtj+bj(zt-z0)xtj,最小化式(7)可得局部线性LADE。minabt=1T|yt-a0-b0(zt-z0)-j=12ajxtj-j=12bj(zt|-z0)xtjKh1(zt-z0)(7)第二,在估计流程的第五步中涉及广义局部线性LADE:mint=1T|yt-a0-b0(zt-zo)-j=12ajxtj-j=12bj(zt-z0)xtjtKh2(zt-z0)(8)当权重t设定为时变方差时(t=gt),估计量拥有最小的渐近方差-协方差矩阵。该方法被称为不可行的局部线性ALADE。第三,可行的
10、局部线性ALADE的权重选择。文献中常将残差的局部核估计值作为权重的替代值,即:gt=i=1nKti|ii=1nKti(9)其中,t是局部线性LADE估计后所得的残差,核函数为:Kti=k(t-itb),ti0,t=i将gt代入式(8)并求解可得到可行的局部线性ALADE。第四,窗宽的选择。在核估计的过程中,窗宽的选择起决定性作用。异常值的存在导致该点所求的残差平方和将占据主导地位,窗宽的选取将异常困难,此处更适合采用绝对值交叉验证法(ACV)进行窗宽选择。在流程中涉及三处窗宽选择,窗宽b的选择参考式(10):bacv=argminACV(b)=1nt=1T|t-g-1t(10)窗宽h的选择参
11、考式(11):hacv=argminACV(h)=1nt=1T|yt-a-t0-j=12a-tjxtj(11)其中,g-1t为式(9)中的去一估计值,a-t0和a-tj是基于局部线性LADE方法得到的去一估计值。2蒙特卡洛模拟本文将模拟含时变波动方差和厚尾特征的函数系数协整模型,生成样本数据以分析局部线性ALADE方法的优劣。考虑下列数据生成过程:yt=(zt+2e(-8z2t)+3sin(2zt)x1t+(4-13(zt-12)2)x2t+tt,t=g(t/n)(12)其中,zt服从(-0.5,0.5)上的均匀分布,x1t、x2t是初始值为0的单位根过程。t分别服从:(1)标准正态分布;(2
12、)位置参数为0、尺度参数为0.2的柯西分布。g(t/n)为时变波动方差,分三种情形讨论:情形1(g1):g21(r)=20+(21-20)I(r)。情形2(g2):g22(r)=20+(21-20)I(r1-)。情形3(g3):g23(r)=20+(21-20)rm。情形1和情形2代表了时变波动方差的突变形式,情形3代表了方差的平滑特征。单变点令取(0.2,0.8)表示突变点发生在前端和后端;双变点和1-的取值表示跳跃发生的时间点,令取(0.2,0.4)表示两次跳跃发生在序列的两端和中间处;平滑时变特征中令m取(1,2)表示变化快慢。同时令0=2,标准差前后变化之比设为=1/0,令的取值为(0
13、.2,0.5,1,5)。样本量分别取 200 和 400,模拟次数设定为 200 次。采用核权最小二乘估计(KLS)、自适应核权最小二乘估计(AKLS)、局部线性LADE(LADE)三种方法做对比,综合分析局部线性ALADE(ALADE)方法的表现。核函数选取二阶EP核,采用局部线性估计以保证方法间的可比性。2.1系数估计情况分析先分析单次样本数据下四种估计方法对模型系数拟合精度的表现。当tN(0,1)且存在时变波动方差时,基于L2准则的系数估计曲线在各类时变波动情形下更加贴近系数真实值,而基于L1准则的曲线拟合效果较差(图略)。含单变点和双变点的时变波动方差时,()取值为(0.2,0.2)时
14、四种估计方法下的系数均与实际值更加贴近,取值为(0.2,5)时各方法下的系数较实际系数曲线偏离严重。含平滑时变波动方差时,(m)取值为(1,0.2)时各估计方法下的系数均贴合实际曲线,取值为(1,5)时各估计方法下24理 论 探 讨统计与决策2023年第21期总第633期的系数偏离实际系数曲线最大。当不存在时变波动(=1)时,AKLS与KLS估计的系数曲线基本重合,说明自适应估计和传统估计效果一致,且明显优于ALADE和LADE。当tC(0,0.2)时,基于L1准则的系数估计曲线更加贴近系数真实值(图略),且ALADE估计方法在部分情况下优于LADE。含单变点和双变点的时变波动方差时,(,)取
15、值为(0.2,0.2)时四种估计方法的系数曲线接近且紧贴实际曲线,取值为(0.2,5)时基于L2方法的估计值偏离实际值较大。含平滑时变波动方差时,(m,)取值为(1,0.2)时四种估计方法的系数与实际值相近,取值为(1,5)时基于L2方法的系数估计值严重偏离实际系数曲线。综上,当残差中含厚尾特征且存在时变波动方差时:就系数估计精度而言,ALADE估计和LADE方法优于AK-LS方法;并且时变波动方差的存在会影响估计效果,且波动值越大影响越大。2.2系数整体拟合精度分析分析可知,当残差服从柯西分布时,基于L1最优准则的估计系数表现相对较好,但ALADE方法与LADE方法间的差异通过图形难以判断。
16、接下来选取加权平均误差(Weighted Average Squared Error,WASE)指标来分析系数整体拟合情况,定义如下:WASE=1nt=1T(a0(zt)-a0(zt))2(range(a0)2+(a1(zt)-a1(zt)2(range(a1)2+(a2(zt)-a2(zt)2(range(a2)2)(13)在模拟过程中,产生200和400的样本数据,并各进行200次蒙特卡洛模拟实验后取WASE的均值和标准差进行对比分析(见表1)。分析发现:第一,对比四种估计方法对系数整体的拟合效果可以发现,LADE和ALADE方法的WASE的均值明显小于KLS和AKLS估计的取值。第二,随
17、着样本的增加,ALADE相比LADE对系数的拟合情况更好:单变点除(,)取值为(0.8,0.2)外,双变点除=0.2情形外,其他各类情况下ALADE估计所得WASE的均值相比LADE估计值都更小,标准差也很接近且波动较小;含平滑变动时变波动方差的各类情况下,除(m)取值(2,0.2),即方差变小且变动较慢外,其他各类情况下ALADE的WASE的均值均更小。整体而言,当残差来自厚尾分布时,ALADE估计对系数的拟合更好且稳健。接下来,从最优窗宽的角度分析AKLS和ALADE两种估计方法对系数的整体拟合效果,通过分析WASE与窗宽h间的函数曲线图得出结论。限于篇幅,仅展示残差服从柯西分布的部分情形
18、(见下页图2)。由图2可知,当残差服从柯西分布时,两种估计方法表1tC(0,0.2)时WASE的均值和标准差T2004000.20.5150.20.550.20.5150.20.55g10.20.20.20.20.80.80.80.20.20.20.20.80.80.8LADE0.24600.35120.52900.95351.09272.74266.356811.21770.70631.98510.82822.20054.946811.38540.17080.25760.35190.56190.68861.10115.159910.66630.47230.64140.53750.78193.
19、31947.0405ALADE0.24220.34680.52610.96511.11762.70836.323511.03560.72272.02290.85412.17384.944010.30720.16840.26930.35130.56420.67361.08325.003810.18380.49790.75220.52760.76862.87686.1699KLS551.15557216.7410652.67987285.35301035.28308356.959013266.240097321.2000810.86807696.6220854.82537736.09006769.
20、874082505.83007.344118.887520.925553.105654.8497168.64241226.53705122.581037.9192171.616044.7223174.7955552.40152947.0590AKLS320.53704026.1540443.34174141.4230736.21015354.032011189.180089332.0200580.32384847.5250645.12434900.25304705.944057151.15006.172713.946118.198845.403949.1737144.7384921.42863
21、461.083027.455192.494934.178397.7403494.16062716.0790g20.20.20.20.20.40.40.40.20.20.20.20.40.40.4LADE0.84902.04470.82321.66141.09272.74263.98418.43751.01012.65931.02472.66882.24904.79210.50121.07480.54071.13600.68861.10113.41116.79800.63161.16740.63511.04031.38122.0634ALADE0.75971.51210.83251.62081.
22、11762.70834.462610.85591.02222.57081.04252.57912.24594.84330.53461.16990.52841.11880.67361.08323.00595.30490.63491.20020.61861.02721.28702.0360KLS776.49537919.7210833.13137943.50701035.28308356.95907540.596069329.89001033.13008360.30101033.48908359.08901111.81108352.241029.7678125.820536.0816128.900
23、654.8497168.6424742.56984258.253049.9748196.613352.8576197.9797295.38201763.1480AKLS478.82764628.2010534.25284671.4180736.21015354.03207219.181069061.3600733.70245356.6630733.99575355.7450805.34855358.299028.9459116.409732.9829117.362349.1737144.7384469.83972232.741038.8572127.006541.5975128.7176285
24、.63231762.6400g3m11112221111222LADE0.52910.98600.68141.40071.09272.74266.040513.13640.64191.39720.78041.77555.381512.36290.36400.54060.45730.74570.68861.10114.37149.21820.43440.64330.49990.74013.92837.9945ALADE0.54240.98700.69231.40471.11762.70835.991610.32750.66531.40930.80251.77704.96199.15000.361
25、80.54720.45220.74560.67361.08324.36259.72760.44200.65570.49480.74363.70117.9945KLS745.90297322.7290803.98027467.88501035.28308356.95908396.275076097.6800840.33017651.4250877.79317749.93306037.633063068.150019.618849.203929.191274.488454.8497168.6424984.76084259.045027.056471.311534.768690.9537819.43
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 函数 系数 模型 局部 线性 估计 方法 改进
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。