主要内容

模拟

蒙特卡洛仿真条件方差模型

描述

例子

V.=模拟(MDL.numobs.模拟A.numobs.- 从完全指定的条件方差模型中的整个条件方差路径MDL.MDL.可以是A.加油贝加奇, 或者GJR.模型。

例子

V.=模拟(MDL.numobs.名称,价值使用一个或多个指定的其他选项模拟条件方差路径名称,价值对论点。例如,您可以生成多个样本路径或指定预先创新路径。

例子

[V.y] =模拟(___此外,使用先前语法中的任何输入参数模拟响应路径。

例子

全部收缩

模拟来自GARCH(1,1)模型的条件方差和响应路径。

使用已知参数指定GARCH(1,1)模型。

mdl = garch('持续的',0.01,'GARCH',0.7,'拱',0.2);

模拟500个样品路径,每个样品路径为100个观察结果。

RNG.默认;重复性的%[v,y] =模拟(MDL,100,'numpaths',500);图形子图(2,1,1)绘制(v)标题('模拟条件差异')子图(2,1,2)绘图(Y)标题('模拟响应'

图包含2个轴。带有标题的轴1仿真条件差异包含500个类型的类型。标题模拟响应的轴2包含500个类型的型号。

模拟响应看起来像静止随机过程中的绘制。

绘制模拟条件方差的2.5,50(中位数)和97.5百分位数。

较低=分子(V,2.5,2);中间=中位数(v,2);upper = proctile(v,97.5,2);图绘图(1:100,较低,'r:',1:100,中间,'K'......1:100,上部,'r:''行宽',2)传奇('95%的间隔''中位') 标题('近似95%的间隔'

图中包含一个坐标轴。标题为“大约95%间隔”的轴包含3个类型为line的对象。这些对象代表95%区间,中值。

由于对条件方差的积极约束,间隔是不对称的。

模拟来自Egarch(1,1)模型的条件方差和响应路径。

使用已知参数指定EGARCH(1,1)模型。

mdl = egarch('持续的',0.001,'GARCH',0.7,'拱',0.2,......'杠杆作用',-0.3);

模拟500个样品路径,每个样品路径为100个观察结果。

RNG.默认;重复性的%[v,y] =模拟(MDL,100,'numpaths',500);图形子图(2,1,1)绘制(v)标题('模拟条件差异')子图(2,1,2)绘图(Y)标题('模拟响应(创新)'

图包含2个轴。带有标题的轴1仿真条件差异包含500个类型的类型。具有标题模拟响应的轴2(创新)包含500个类型的类型。

模拟响应看起来像静止随机过程中的绘制。

绘制模拟条件方差的2.5,50(中位数)和97.5百分位数。

较低=分子(V,2.5,2);中间=中位数(v,2);upper = proctile(v,97.5,2);图绘图(1:100,较低,'r:',1:100,中间,'K'......1:100,上部,'r:''行宽',2)传奇('95%的间隔''中位') 标题('近似95%的间隔'

图中包含一个坐标轴。标题为“大约95%间隔”的轴包含3个类型为line的对象。这些对象代表95%区间,中值。

由于对条件方差的积极约束,间隔是不对称的。

模拟来自GJR(1,1)模型的条件方差和响应路径。

使用已知参数指定GJR(1,1)模型。

mdl = gjr('持续的',0.001,'GARCH',0.7,'拱',0.2,......'杠杆作用',0.1);

模拟500个样品路径,每个样品路径为100个观察结果。

RNG.默认;重复性的%[v,y] =模拟(MDL,100,'numpaths',500);图形子图(2,1,1)绘制(v)标题('模拟条件差异')子图(2,1,2)绘图(Y)标题('模拟响应(创新)'

图包含2个轴。带有标题的轴1仿真条件差异包含500个类型的类型。具有标题模拟响应的轴2(创新)包含500个类型的类型。

模拟响应看起来像静止随机过程中的绘制。

绘制模拟条件方差的2.5,50(中位数)和97.5百分位数。

较低=分子(V,2.5,2);中间=中位数(v,2);upper = proctile(v,97.5,2);图绘图(1:100,较低,'r:',1:100,中间,'K'......1:100,上部,'r:''行宽',2)传奇('95%的间隔''中位') 标题('近似95%的间隔'

图中包含一个坐标轴。标题为“大约95%间隔”的轴包含3个类型为line的对象。这些对象代表95%区间,中值。

由于对条件方差的积极约束,间隔是不对称的。

模拟每日纳斯达克综合指数的条件差异返回500天。使用模拟进行预测,并近似95%的预测间隔。比较GARCH(1,1),EGARCH(1,1)和GJR(1,1)的预测。

加载工具箱中包含的NASDAQ数据。将索引转换为返回。

加载data_equityidx.纳斯达克= DataTable.nasdaq;r = price2ret(纳斯达克);t =长度(r);

适合GARCH(1,1),EGARCH(1,1)和GJR(1,1)模型到整个数据集。推断有条件的差异用作预测模拟的预先条件差异。

mdl = cell(3,1);%preallocation.MDL {1} = GARCH(1,1);mdl {2} = egarch(1,1);MDL {3} = GJR(1,1);estmdl = cellfun(@(x)estmate(x,r,'展示''离开'),MDL,......'统一输出',错误的);v0 = Cellfun(@(x)推断(x,r),estmdl,'统一输出',错误的);

estmdl.是3×1细胞矢量。每个单元是一种不同类型的估计条件方差模型,例如,estmdl {1}是估计的GARCH(1,1)模型。v0.是一个3×1个细胞矢量,每个细胞包含来自相应的估计模型的推断条件差异。

模拟1000个样本路径,每个样本都有500个观察结果。使用观察到的返回并将条件差异推断为预先数据。

VSIM = CELL(3,1);%preallocation.为了J = 1:3 RNG默认;重复性的%vsim {j} =模拟(estmdl {j},500,'numpaths',1000,'e0',r,'v0',v0 {j});结尾

VSIM是一个3×1个细胞矢量,每个单元包含从相应的估计模型产生的500×1000矩阵的模拟条件差异。

绘制模拟平均预测和近似95%的预测间隔,以及从数据推断的条件方差。

较低= Cellfun(@(x)分色(x,2.5,2),VSIM,'统一输出',错误的);Upper = Cellfun(@(x)prctile(x,97.5,2),VSIM,'统一输出',错误的);Mn = Cellfun(@(x)均值(x,2),vsim,'统一输出',错误的);Datesplot =日期(结束 -  250:结束);datesfh =日期(结束)+(1:500)';h =零(3,4);数字为了j = 1:3 col =零(1,3);Col(j)= 1;h(j,1)= plot(datesplot,v0 {j}(ex-250:结束),'颜色',col);抓住h(j,2)= plot(datesfh,mn {j},'颜色',col,'行宽',3);h(j,3:4)= plot([datesfh datesfh],[下{j} upper {j}],':'......'颜色',col,'行宽',2);结尾HGCA = GCA;绘图(DatesFH(1)* [1 1],HGCA.YLIM,'k-');dateTick;轴紧的;H = H(:,1:3);传奇(H(:),'GARCH  - 推断''egarch-推断''GJR  - 推断'......'garch  -  sim。意思''egarch  -  sim。意思''gjr  -  sim。意思'......'GARCH  -  95%。int。''egarch  -  95%。int。'......'GJR  -  95%。int。''地点''东北') 标题('模拟条件方差预测') 抓住离开

图中包含一个坐标轴。标题仿真条件方差预测的轴包含13个类型的类型。这些物体代表GARCH  - 推断,GARCH  -  SIM。平均值,garch  -  95%前。int。,egarch  - 推断,egarch  -  sim。平均值,肉食 -  95%。int。,gjr推断,gjr  -  sim。平均值,GJR  -  95%前。int ..

输入参数

全部收缩

没有任何未知参数的条件方差模型,指定为a加油贝加奇, 或者GJR.模型对象。

MDL.不能包含任何属性价值。

示例路径长度,指定为正整数。也就是说,每个输出路径生成随机观测的数量。V.ynumobs.行。

名称值对参数

指定可选的逗号分离对名称,价值论点。姓名是参数名称和价值是相应的价值。姓名必须出现在引号内。您可以以任何顺序指定多个名称和值对参数name1,value1,...,namen,valuen

例子:'numpath',1000,'e0',[0.5;0.5]指定生成1000样品路径和使用[0.5;0.5]作为每个路径的预先创新。

生成的样本路径数量,指定为逗号分隔对组成'numpaths'和一个正整数。V.ynumpaths.列。

例子:'numpath',1000

数据类型:双倍的

预先创新,指定为逗号分隔对组成'e0'和数字列向量或矩阵。预先创新为条件方差模型的创新过程提供了初始值MDL.。预先创新从平均值的分布导出。

E0.必须包含至少mdl.q.元素或行。如果E0.包含额外的行,模拟用最新mdl.q.只要。

最后一个元素或行包含最新的预先创新。

  • 如果E0.是一个列向量,它代表了底层创新系列的单个路径。模拟适用E0.每个模拟路径。

  • 如果E0.是矩阵,然后每个列代表底层创新系列的预先样路径。E0.必须至少有numpaths.列。如果E0.有更多的列,比必要,模拟使用第一个numpaths.仅限列。

默认值为:

  • 对于GARCH(P.问:)和gjr(P.问:) 楷模,模拟将任何必要的预先创新设置为具有平均零和标准偏差的独立干扰序列等于条件方差过程的无条件标准偏差。

  • 对于egarch(P.问:) 楷模,模拟将任何必要的预先创新设置为独立的干扰序列,其平均零和方差等于EGARCH方差过程的对数的指数无条件平均值。

例子:'e0',[0.5;0.5]

正面预定的条件方差路径,指定为数字矢量或矩阵。v0.为模型中的条件差异提供初始值。

  • 如果v0.是一列栏矢量,然后模拟将其应用于每个输出路径。

  • 如果v0.是一个矩阵,那么它必须至少有numpaths.列。如果v0.有更多的列,比必要,模拟使用第一个numpaths.仅限列。

  • 对于GARCH(P.问:)和gjr(P.问:) 楷模:

    • v0.必须至少有mdl.p.初始化方差方程的行。

    • 默认情况下,模拟将任何必要的预先规定差异设置为条件方差过程的无条件方差。

  • 对于egarch(P.问:) 楷模,模拟

    • v0.必须至少有max(mdl.p,mdl.q)初始化方差方程的行。

    • 默认情况下,模拟将任何必要的预先规范差异设置为eGARCH方差过程对数的指数无条件平均值。

如果行的数量v0.超过必要的数字模拟仅使用最新,所需的观察数。最后一个元素或行包含最新的观察。

例子:'v0',[1;0.5]

数据类型:双倍的

笔记

  • 如果E0.v0.列向量,模拟将它们应用于输出的每一列V.y。此应用程序允许模拟路径共享对蒙特卡罗模拟预测和预测误差分布的公共起点。

  • s表示缺失值。模拟删除缺失的值。该软件合并预定数据(E0.v0.),然后使用列表明智删除以删除包含至少一个的任何行。删除数据中的S在数据中减少了样本大小。删除纳尼斯也可以创建不规则的时间序列。

  • 模拟假设您同步预定数据,使得每个预定系列的最新观察同时发生。

输出参数

全部收缩

模拟条件方差路径与相关的平均创新y,作为数字列向量或矩阵返回。

V.是A.numobs.-经过-numpaths.矩阵,其中每列对应于模拟条件方差路径。行V.是对应于周期性的时期MDL.

模拟响应路径,作为数字列向量或矩阵返回。y通常代表一个平均零,异源间时间序列具有条件差异的创新V.(继续预先创新系列E0.)。

y还可以代表平均零的时间序列,异源型创新加上偏移量。如果MDL.包括偏移量,然后模拟将偏移量添加到底层平均值零,异源基础创新y代表了偏移调整的创新的时间序列。

y是A.numobs.-经过-numpaths.矩阵,其中每列对应于模拟响应路径。行y是对应于周期性的时期MDL.

参考文献

[1] Bollerslev,T。“广义归类疾病性瘢痕糕点”。中国经济学杂志。卷。31,1986,第307-327页。

[2] Bollerslev,T。“一个有条件的异源性时间序列模型,用于投机价格和回报率。”经济学与统计审查。卷。69,1987,第542-547页。

[3]盒子,G.E.P.,G.M. Jenkins和G. C. Reinsel。时间序列分析:预测和控制。3 ed。Englewood Cliffs,NJ:Prentice Hall,1994年。

[4]内酯,W。应用计量计量时间序列。霍博肯,NJ:John Wiley&Sons,1995。

[5] Engle,R. F.“归类化条件异性瘢痕贴膜,估计联合王国通胀差异。”Moveryetrica.。卷。50,1982,pp。987-1007。

[6] Grosten,L. R.,R.Jagannathan和D. E. Runkle。“关于预期价值与标称超额股票的波动性关系。”财务杂志。卷。48,第5,1993,第5,199,PP。1779-1801。

[7]汉密尔顿,J.D。时间序列分析。普林斯顿,新泽:普林斯顿大学出版社,1994年。

[8]纳尔逊,D.B。“资产中的条件异质娱乐性回报:一种新方法。”Moveryetrica.。卷。59,1991,第347-370页。

在R2012A介绍