|
当y<0很小时,f(y)趋于正,如果y远低于零,f(y)趋于大。换句话说,如果Yi是正的,在expection上,Yi+1将是负的,如果Yi是负的,在expection上,Yi+1将是正的。对于波动性g,其gra ph显示“U形”笑脸。3数据校准局部多项式回归[27][8]用于估计漂移(f)和波动(g)f函数。假设{Yt}是一个时间序列,我们想要建立一个ARCHmodel,使Yi=f(Yi)-1) +g(易)-1) i,其中i是i.i.d.N(0,1)个随机变量。程序如下:我们寻找一个满足g(x)=E(Yi | Yi)的函数g(x)-1=x)- E(Yi | Yi)-1=x)。-0.1-0.05 0 0.05 0.1 0.15-0.01-0.00500.0050.010.0150.020.0250.03fYi(a)f-0.1-0.05 0.05 0.1 0.150.010.0150.020.0250.030.0350.040.045gYi(b)g-0.1-0.05 0.05 0.1 0.1500.20.40.60.811.21.41.61.8x 10-3g2Yi(c)g图4:Drif t和volatility function的大小为{Yt},然后对于每个x,考虑以下两个最小化问题:[α(x),α(x),α(x)]=arg minα,α,αnXi=1Yi- α- α易-1.- xh-α易-1.- xh!K易-1.- xh,(4) [β(x),β(x),β(x)]=arg minβ,β,βnXi=1Yi- β- β易-1.- xh-β易-1.- xh!K易-1.- xh,(5) 其中,K(·)表示一个非负权重函数,h是一个正数,称为带宽。我们选择K作为标准的普通pdf。通过^g(x)=α(x)估计g(x)-β(x),f的估计由^f(x)=β(x)给出。设{Pi}为S&P500指数的日数据,其中最后一个日期为2013年5月2日,设{Yi}为收益率过程,即Yi=log Pi- 对数Pi-1.对于每个数据大小S,我们选择带宽h为h={max(Yi)-min(Yi)}/γ,其中所有Yi都属于这个集合,γ是一个选定的正常数。
|