RS类方法估计Hurst指数的有效性检验Word文档格式.docx
- 文档编号:16041159
- 上传时间:2022-11-17
- 格式:DOCX
- 页数:14
- 大小:114.52KB
RS类方法估计Hurst指数的有效性检验Word文档格式.docx
《RS类方法估计Hurst指数的有效性检验Word文档格式.docx》由会员分享,可在线阅读,更多相关《RS类方法估计Hurst指数的有效性检验Word文档格式.docx(14页珍藏版)》请在冰豆网上搜索。
ThispaperstudiestheefficiencyofestimatingHurstindexwithR/Stypeincludingclassicalrescaledrangeanalysis,modifiedrescaledrangeanalysisandrescaledvarianceanalysisbysimulatingFGNserieswithDHM.Thereareevidentdifferencesbetweentheeffectsofthelengthofseries,thetypeofdealingwithshort-rangdependence,strengthofaddedwhitenoiseinseriesonthreeestimators,moreover,theyhaveseparateapplicableintervalofestimatingHurstindex.
Keywords:
HurstIndex;
ClassicalRescaledRangeAnalysis;
ModifiedRescaledRangeAnalysis;
RescaledVarianceanalysis;
Efficiency.
①资助项目:
国家自然科学基金项目(项目名称:
资本市场的分形结构及其应用研究,项目编号:
70501034,项目负责人:
黄诒蓉。
1引言
Hurst(1951[1]在长期的水文研究工作中发现,河流流量存在较强的长记忆性。
后来,许多研究发现,该特性不仅存在于自然界,而且广泛存在于经济与管理领域的数据中。
金融时间序列长记忆性的检测与建模目前已经成为金融计量领域研究的重要内容。
长记忆性通常从具有双曲率缓慢衰减形式的自相关函数和在零频率处趋于无限值的谱密度函数等两个角度进行刻画,而且均通过Hurst指数(下文简称H指数来表征长记忆性程度。
因此,H指数的估计与检验是长记忆性研究的关键工作,其估计与检验的有效性直接影响对长记忆性的甄别。
尽管现有文献提出多种估计H指数的方法,但是R/S类方法由于其简洁性而一直受到研究者的青睐,迄今为止仍是估计H指数最常用的方法。
该类方法最早是由Hurst(1951[1]和Mandelbrot&
Wallis(1969[2]提出的经典R/S分析方法(ClassicalRescaledRange,下文简称CRS,后来Lo(1991[3]通过考虑序列的短期相关性对CRS进行了修正而提出修正R/S分析方法(ModifiedRescaledRange,下文简称MRS,Giraitisetal.(2003[4]利用部分和序列的方差替代部分和序列的极差对CRS进行了修正而提出V/S分析方法(RescaledVariance,下文简称VS。
虽然MRS和VS均是从不同方面对CRS进行了某种程度的修正,但是能否真正提高H指数的估计有效性呢?
尽管许多研究文献对R/S类分析方法的有效性提出许多质疑,但是目前大多数的研究都是这些估计方法在某些领域的直接应用,而对它们估计H指数的有效性的关注和研究却很少。
R/S类分析方法由于产生于Hurst的长期实践工作过程而没有严格的数学推导,因此,我们难以使用严格数学证明来检验其有效性,但是可以通过模拟研究方法来综合评价它们的有效性。
鉴于此,为检验R/S类方法估计H指数的有效性,下面将通过预先设定不同的H真实值并模拟出一系列的FGN序列,然后利用CRS、MRS和VS三种方法分别估计H指数,从H估计值的均值和标准差两方面反映序列长度、短期相关性处理、白噪声成分等因素对H指数估计的影响,并综合比较三种估计量的有效性,为R/S类分析方法使用者提供有益的参考。
2研究方法
2.1R/S类分析方法
通常将CRS、MRS和VS统称为R/S类分析方法,原因在于它们具有共同的基本原理,我们可将它们估计H指数的基本步骤概括如下:
第1步,将原始序列划分成若干子序列。
合适选定某标度长度n,将原始序列{}T
1tty=划分成M个互不重叠的长度为n的子序列:
y,,y,y(n21,y,,y,y(n22n1n++,…,y,,y,y(Mn2n1M(1n1M(+-+-,其中,]n/T[M=,[]⋅表示取整,下同;
第2步,计算出与标度长度n所对应的标度统计量。
首先计算出根据标度长度n所划分出的所有子序列的局部统计量m,nLS,然后利用所有局部统计量计算出标度长度n对应的标度统计量nSS,此处标度统计量是所有局部统计量的平均值,即为∑==M
1
mm,nnLSM1SS;
第3步,改变标度长度n,重复前面的步骤,这样得到一系列的标度长度n及其相应的标度统计量序列nSS;
第4步,若标度统计量序列与标度长度序列存在如下的标度关系:
HnCn~SS,其中,C为某常数,则可以采用如下的双对数回归方式估计H指数:
nlog(HClog(SSlog(n+=。
(1
CRS(Hurst,1951[1];
Mandelbrot,1969[2]使用如下的局部统计量为:
m,nm,nm,nS/RS/R(=,(2其中,]yy(minyy(max[Rm,nk
1tm,tnk1m,nk1tm,tnk1m,n---=∑∑=≤≤=≤≤为第m个子样本的极差,2m,nn1
tm,tm,nyy(n1S-=∑=为第m个子样本的标准差,∑==n1jm,tm,nyn1y为第m个子样本的均值。
MRS(Lo,1991[3]使用如下的局部统计量:
q(/R/Rm,nm,nnσ=σ(,(3
其中,m,nR同上述CRS中的定义,q(m,nσ由下式定义:
jq
jjy
njtmnmjtmnmtqjjntmnmtqyyyyqnyynqγωσωσ(2(((2(1(121,,,,121,,2
mn,∑∑∑∑=+=-==+=⎥⎦⎤⎢⎣⎡--+-=,(3其中,2n1tm,nm,t2y
yy(n1ˆ∑=-=σ,当nq<
时,1qj1q(j+-≡ω为Bartlett权重,2
m,nσ和jγ为第m个子样本{}n1
tm,ty=的样本方差和j阶样本自协方差。
(nq的选择至关重要,其直接影响到检验效果。
Lo建议使用如下的q值最优选择公式:
⎥⎥⎦
⎤⎢⎢⎣⎡⎪⎪⎭⎫⎝⎛-⎪⎭⎫⎝⎛=3/23/1ˆ1ˆ223ρρnqopt,(4其中,n为子样本长度,ρ
ˆ为一阶自相关系数的估计值,[]⋅表示取整。
当0q=时,MRS与CRS方法一致。
VS(Giraitisetal.,2003[4]使用的局部统计量为:
m
nm
nm,nSVS/V(=,(5其中,m,nV为第m个子序列部分和序列{}n1kkS=的方差,具体展开形式为:
⎥⎥⎦⎤⎢⎢⎣⎡⎪⎭⎫⎝⎛--⎪⎭⎫⎝⎛-=-=∑∑∑∑∑=====2n1kn1kk1tm,nm,t2k1tm,nm,tn1k2nkm
nyy(n1yy(n1SS(n1V,其它与上述MRS中的定义相同。
2.2模拟设计
由Mandelbrot&
VanNess(1968[5]提出的分形高斯噪声(FractionalGaussian
Noise,简称FGN是第一个完整的长记忆模型,通常定义为具有H指数的分形布朗运动t(BH的一阶差分过程{}∞
=-+=0tHHttt(B1t(By|y,其第k阶自协方差
函数为:
(H2H2H22y|1k||k|2|1k|2
k(++--σ=γ,(6其中,,1,0,1,k-=,y(Vart2=σ为任意正数,1H0<
<
为H指数。
为模拟上述定义的FGN序列,许多学者提出了多种模拟算法,比如,基于条件分布的Durbin-Levinson方法(Brockwell&
Davis,1991[6]、基于循环嵌入矩阵和傅立叶变换的模拟方法(Davies&
Harte,1987[7];
Wood&
Chan,1994[8]、基于小波的合成模拟方法(Abtry&
Sellan,1996[9]等。
我们通过反复模拟比较发现,这些模拟方法的效果基本一致。
因此,本文选择Davies&
Harte(1987
[7]提出的算法(Davies&
HarteMethod,简称DHM模拟FGN序列。
在式(1给定自协方差序列T
0kk}{=γ的条件下,DHM算法模拟长度为T的FGN样本序列
T
1tt}y{=的基本步骤如下:
第1步,计算序列1T2
0kT,k}A{-=:
T/kji1
T21TjjT2T/kjiT0jjT,kee
Aπ--+=-π-=∑∑γ+γ=,(7
其中,1i2-=,该步其实是对自协方差序列11TT1T10,,,,,,,γγγγγγ--实施离散傅立叶变换(DFT;
第2步,检查非负条件:
对所有k,0AT,k≥,该步也很重要,若不满足非负条件,则模拟序列是无效的;
第3步,模拟产生均值为0、方差为1的独立高斯随机变量序列120}{-=TttZ;
第4步,计算复值序列1T2
0kk}V{-=:
⎪⎪⎪⎩⎪⎪⎪⎨⎧-≤<
=<
≤+==---,
1T2kT,V,Tk,ZTA2,Tk1,iZZ(TA,0k,ZTA2V*kT21T2T,Tk21k2T,k0T,0k,(8
其中,*kT2V-表示kT2V-的共轭复数。
第5步,通过离散傅里叶逆变换(IDFT可得到模拟序列T1tt}y{=:
∑-=-π=1
T20
为综合对比3种估计量估计H指数的有效性,首先需要对模拟参数做出设定。
我们选择0.1、0.2、0.3、0.4、0.5、0.6、0.7、0.8、0.9等9个数值作为H指数的真实值,这些取值范围基本覆盖了序列的各种特性:
当H真实值介于0与0.5之间时,序列具有反持续性(或均值回复;
当H真实值等于0.5时,序列没有记忆性;
当H真实值介于0.5与1之间时,序列具有长记忆性(或状态持续性。
当考察序列长度对H估计值的影响时,我们选定72、8
2、9
2、10
2、11
2、122、13
2、
142等8种长度,这些长度基本可覆盖目前大多数分析数据的长度范围,当考察
相关性处理、白噪声成分对H估计值的影响时,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- RS 方法 估计 Hurst 指数 有效性 检验