穩定分布的参数估计.docVIP

稳定分布的参数估计 顾娟 茆诗松 2012-9-12 10:25:00　　来源：《应用概率统计》2002年第4期 　　内容提要：由于金融数据经常具有“高峰厚尾”现象，所以用稳定分布去拟合。但由于稳定分布没有密度函数显式，而且可能一阶矩或二阶矩又不存在，因此稳定分布的参数估计问题用经典方法很难处理。本文利用类似Duffie和Singleton(1993)的模拟矩法(SMM)的想法，构造了一种新的参数估计方法，并得到该估计的强相合性结果。最后举了一个实际的例子，分析了深圳成分指数的情况。　　关键词：稳定分布 模拟矩法(SMM) 强相合性 特征函数　　作者简介：顾娟广发证券股份有限公司博士后工作站，广州510075；茆诗松华东师范大学统计系，上海200062　　引言　　金融资产的收益是大量在实践上连续到来的信息和个人决策的结果。根据中心极限定理，如果经适当的平移和尺度变换，大量的独立同分布的随机变量和的极限分布一定是稳定分布族中的某一分布。因此很自然地认为，如果资产收益的积累是可加的，则应近似服从稳定分布；如果资产收益的累积是相乘的，则应近似服从对数稳定分布。　　对于稳定分布的研究从本世纪初就引起了大家的注意。正态分布属于稳定分布族。同时也是大家最为熟悉也最易处理的分布，所以在处理一般问题时大都认为资产收益的分布是近似服从正态分布或对数正态分布。在处理股票价格的经验研究中，Mandelbrot[1]和Fama[2]最先发现股票价格的对数变化明显偏离了正态分布，尤其是非常小和非常大的价格变化比在正态假设出现的频率明显偏大、对这类问题用非高斯的稳定分布去拟和优于正态分布。（近年来，很多金融数据的经验研究发现，金融数据具有一些普遍的规律，其中与平稳直接相关的就是“高峰厚尾”现象，这种现象现在已得到广泛的认识）这之后，很多统计学，经济，金融学等方面的学者对稳定分布和资产收益之间的关系进行了研究。　　我们考察了深圳成分指数从1994年1月3日到1998年11月2日的日收益率情况，发现经过作收益率变化后，日收益率的分布形态明显具有这种“高峰厚尾”现象。其中表1是回收益率的常用描述性统计量的估计值，其中的曲线是用正态分布拟和的结果。可以看出在上述时期，深圳市场的日平均收益率大于0，折算成年收益率大约为25%左右。考虑厚尾的衡量指标常常用峰度，即四阶矩。由于正态分布的峰度为3，而该数据的峰度为17.7618，说明该数据的峰度远大于正态分布的峰度，这正是“高峰”说法的来源。另外，数据比正态分布不仅峰度高，而且两端数据出现的频率也较大。所以我们试图用稳定分布去拟和深圳数据。 　　 　　　　对于稳定分布，没有显示的密度函数或分布函数形式。Levy(1924)得到了稳定分布特征函数的显示表达。经过Zolotarev[3]的再参数化变成极点形式，得到了目前比较常用的稳定分布的特征函数的表达式。　　　　由于稳定分布没有密度函数显式，因此对它的研究比普通的密度函数要困难得多。但人们还是通过很多方法得到了它的一些性质。比如：它的密度函数是连续的，并且是单峰的；它的小于α阶矩存在，大于α阶矩不存在等等。　　这些对稳定分布的研究，到七十年代有很大的进展。但由于没有显示的密度函数形式，而且可能一阶矩或二阶矩又不存在，因此稳定分布的参数估计问题受到了很大的阻碍，对平稳原理的进一步推广和使用也受到了限制。八十年代对稳定分布的研究速度逐渐慢了下来。直到九十年代，由于计算机，计算方法和计算速度有了突飞猛进的发展，世界经济的发展也日益繁荣，对稳定分布的研究又引起了人们极大的兴趣。对稳定分布的估计也有了一些新方法。　　本文利用类似Duffie和Singleton(1993)的模拟矩法(SMM)的想法，构造了一种新的参数估计方法，并研究了估计的一些性质。这个内容放在一、二、三给出了模拟结果，并用深圳成分指数的例子说明。　　一、参数估计方法　　记稳定分布为S(α,β,γ,δ)，由四个参数所决定。其中γ为位置参数，δ为尺度参数。若随机变量x～S(α,β,γ,δ)，则[(x-γ)/δ]～S(α,β,0,1)。α,β是形状参数，其中α称为特征指数，决定分布函数尾部的厚薄，α离2越远，说明分布的尾部比正态尾部越厚，比正态峰度越大，当α=2时，就是正态分布。β称为偏度系数，当β=0时，分布函数就是对称形式。在稳定分布族中，我们最熟悉的就是当α=2时的分布，也即正态分布。另外当α=1,β=0时，稳定分布就是柯西分布。经过标准化变换后的稳定分布S(α,β,0.1)称为标准稳定分布，简记为S(α,β)。我们下面考虑的是标准稳定分布的参数估计问题。　　我们的想法非常类似模拟矩方法(Simulated Moment Method，简记为SMM）。