主要内容

模拟

条件方差模型的蒙特卡罗模拟

描述

例子

V=模拟(Mdl,numObs模拟了一个numObs-周期条件方差路径从完全指定的条件方差模型MdlMdl可以是garch,egarch,或gjr模型。

例子

V=模拟(Mdl,numObs,名称,值用一个或多个指定的附加选项模拟条件方差路径名称,值配对参数。例如,您可以生成多个采样路径或指定采样前创新路径。

例子

V,Y) =模拟(___另外,使用前面语法中的任何输入参数模拟响应路径。

例子

全部折叠

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

指定一个参数已知的GARCH(1,1)模型。

Mdl = garch (“不变”,0.01,“四国”, 0.7,“拱”, 0.2);

模拟500个样本路径,每个路径有100个观测值。

rng默认的%为了再现性[V, Y] =模拟(Mdl, 100,“NumPaths”, 500);图subplot(2,1,1) plot(V) title(“模拟条件方差”)子地块(2,1,2)地块(Y)名称(“模拟反应”

图中包含2个轴对象。具有标题模拟条件方差的轴对象1包含500个line类型的对象。具有标题模拟响应的轴对象2包含500个line类型的对象。

模拟的响应看起来像是从平稳的随机过程中提取的。

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

低= prctile (V, 2.5, 2);中间=值(V, 2);上= prctile (V, 97.5, 2);图绘制(1:10 0,降低,“:”1:10 0,中间,“k”,...1:10 0,上,“:”,“线宽”2)传说(“95%间隔”,“中值”)头衔(“大约95%的间隔”

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

由于条件方差的正约束,区间是不对称的。

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

指定一个参数已知的EGARCH(1,1)模型。

Mdl = egarch (“不变”, 0.001,“四国”, 0.7,“拱”, 0.2,...“杠杆”, -0.3);

模拟500个样本路径,每个路径有100个观测值。

rng默认的%为了再现性[V, Y] =模拟(Mdl, 100,“NumPaths”, 500);图subplot(2,1,1) plot(V) title(“模拟条件方差”)子地块(2,1,2)地块(Y)名称(的模拟反应(创新)

图中包含2个轴对象。标题为“模拟条件方差”的轴对象1包含500个类型为line的对象。标题为“模拟响应(创新)”的轴对象2包含500个类型为line的对象。

模拟的响应看起来像是从平稳的随机过程中提取的。

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

低= prctile (V, 2.5, 2);中间=值(V, 2);上= prctile (V, 97.5, 2);图绘制(1:10 0,降低,“:”1:10 0,中间,“k”,...1:10 0,上,“:”,“线宽”2)传说(“95%间隔”,“中值”)头衔(“大约95%的间隔”

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

由于条件方差的正约束,区间是不对称的。

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

指定具有已知参数的GJR(1,1)模型。

Mdl=gjr(“不变”, 0.001,“四国”, 0.7,“拱”, 0.2,...“杠杆”, 0.1);

模拟500个样本路径,每个路径有100个观测值。

rng默认的%为了再现性[V, Y] =模拟(Mdl, 100,“NumPaths”, 500);图subplot(2,1,1) plot(V) title(“模拟条件方差”)子地块(2,1,2)地块(Y)名称(的模拟反应(创新)

图中包含2个轴对象。标题为“模拟条件方差”的轴对象1包含500个类型为line的对象。标题为“模拟响应(创新)”的轴对象2包含500个类型为line的对象。

模拟的响应看起来像是从平稳的随机过程中提取的。

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

低= prctile (V, 2.5, 2);中间=值(V, 2);上= prctile (V, 97.5, 2);图绘制(1:10 0,降低,“:”1:10 0,中间,“k”,...1:10 0,上,“:”,“线宽”2)传说(“95%间隔”,“中值”)头衔(“大约95%的间隔”

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

由于条件方差的正约束,区间是不对称的。

模拟纳斯达克综合指数500天的每日回报的条件方差。使用模拟来作出预测和大约95%的预测间隔。比较GARCH(1,1)、EGARCH(1,1)和GJR(1,1)拟合的预测结果。

加载工具箱中包含的纳斯达克数据。将索引转换为返回值。

负载数据设备IDX纳斯达克= DataTable.NASDAQ;r = price2ret(纳斯达克);T =长度(r);

对整个数据集拟合GARCH(1,1)、EGARCH(1,1)和GJR(1,1)模型。推断条件方差作为预测模拟的前样本条件方差。

Mdl =细胞(3,1);%预先配置Mdl {1} = garch (1,1);Mdl {2} = egarch (1,1);Mdl {3} = gjr (1,1);EstMdl = cellfun (@ (x)估计(x, r,“显示”,“关”)、Mdl...“UniformOutput”、假);v0 = cellfun (@ (x)推断(x, r), EstMdl,“UniformOutput”、假);

EstMdl是3 × 1的细胞载体。每个单元格是一种不同类型的估计条件方差模型,例如:EstMdl {1}是一个估计的GARCH(1,1)模型。是一个3乘1的单元向量,每个单元包含从相应的估计模型推断出的条件方差。

模拟1000个样本路径,每个路径有500个观察值。使用观察到的返回值和推断出的条件方差作为样本前数据。

vSim =细胞(3,1);%预先配置2 . n .复合体,复合体默认的%为了再现性vSim {j} =模拟(EstMdl {j}, 500年,“NumPaths”,1000,‘E0’R“V0”v0 {j});结束

vSim是一个3 × 1的单元向量,每个单元包含一个500 × 1000的模拟条件方差矩阵,由相应的估计模型生成。

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

下=cellfun(@(x)prctile(x,2.5,2),vSim,“UniformOutput”、假);上= cellfun (@ (x) prctile (x, 97.5, 2), vSim,“UniformOutput”、假);mn = cellfun (@ (x)的意思是(x, 2), vSim,“UniformOutput”、假);datesPlot = date (end - 250:end); / /结束日期datesFH = date (end) + (1:500)';h = 0(3、4);数字J = 1:3 col = 0 (1,3);坳(j) = 1;h (j, 1) =情节(datesPlot, v0 {j}(端- 250:端),“颜色”,卡扎菲);持有h(j,2)=绘图(日期fh,mn{j},“颜色”上校,“线宽”,3); h(j,3:4)=绘图([datesFH datesFH],[lower{j}upper{j}],':',...“颜色”上校,“线宽”2);结束hGCA=gca;绘图(日期h(1)*[1],hGCA.YLim,“k——”)坐标轴紧的;h = h (:, 1:3);传奇(h (:),“GARCH-推断”,“EGARCH-推断”,“GJR -推断”,...“GARCH - Sim卡。的意思是“,“EGARCH - Sim卡。的意思是“,“GJR - Sim卡。的意思是“,...“GARCH——95% Fore。”Int。,“EGARCH——95%的福尔。Int。,...“GJR - 95%福尔。Int。,“位置”,“东北”)头衔(“模拟条件方差预测”)举行

图中包含一个轴对象。标题为“模拟条件方差预测”的轴对象包含13个线型对象。这些对象代表GARCH - infer, GARCH - Sim。平均GARCH - 95% Fore。Int。,EGARCH - Inferred, EGARCH - Sim. Mean, EGARCH - 95% Fore. Int., GJR - Inferred, GJR - Sim. Mean, GJR - 95% Fore. Int..

输入参数

全部折叠

无任何未知参数的条件方差模型,指定为garch,egarch,或gjr模型对象。

Mdl不能包含任何具有价值。

样本路径长度,指定为正整数。也就是说,每个输出路径生成的随机观测数。VYnumObs排。

名称-值参数

指定可选的逗号分隔的对名称,值论据。的名字参数名和价值为对应值。的名字必须出现在引号内。可以以任意顺序指定多个名称和值对参数名称1,值1,…,名称,值

例子:“E0”“numPaths”,1000年,[0.5;0.5)指定生成1000示例路径和使用[0.5; 0.5]作为每条路径创新的前样例。

要生成的样本路径数,指定为逗号分隔对,由“NumPaths”和一个正整数。VY努帕斯列。

例子:“NumPaths”,1000年

数据类型:双重的

前样例创新,指定为逗号分隔对组成‘E0’和数值列向量或矩阵。采样前新值为条件方差模型的新值过程提供初始值Mdl.样本创新源于均值为0的分布。

E0必须至少包含Mdl。问元素或行。如果E0包含额外的行,模拟使用了最新的Mdl。问只有。

最后一个元素或行包含最新的前样例创新。

  • 如果E0是一个列向量,它表示潜在创新系列的单一路径。模拟适用于E0到每个模拟路径。

  • 如果E0是一个矩阵,则每列表示基础创新系列的预采样路径。E0必须至少努帕斯列。如果E0列比需要的多,模拟使用第一个努帕斯仅列。

默认值是:

  • 为了GARCH(P,)和GJR (P,)模型,模拟将任何必要的采样前创新设置为独立的干扰序列,平均值为零,标准偏差等于条件方差过程的无条件标准偏差。

  • EGARCH (P,)模型,模拟将任何必要的采样前创新设置为一个独立的干扰序列,其平均值为零,方差等于EGARCH方差过程对数的指数无条件平均值。

例子:‘E0’,[0.5;0.5]

正预采样条件方差路径,指定为数字向量或矩阵。提供模型中条件方差的初始值。

  • 如果是一个列向量,那么模拟将其应用于每个输出路径。

  • 如果是一个矩阵,那么它必须至少有努帕斯列。如果列比需要的多,模拟使用第一个努帕斯仅列。

  • 为了GARCH(P,)和GJR (P,)模型:

    • 必须至少Mdl。P行初始化方差方程。

    • 默认情况下,模拟将任何必要的样本前方差设置为条件方差过程的无条件方差。

  • EGARCH (P,)模型,模拟

    • 必须至少max (Mdl.P Mdl.Q)行初始化方差方程。

    • 默认情况下,模拟将任何必要的采样前方差设置为EGARCH方差过程对数的指数无条件平均值。

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

例子:“半”,[1;0.5)

数据类型:双重的

笔记

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

  • S表示缺失值。模拟删除缺失值。该软件合并前样本数据(E0),然后使用按列表删除的方法删除包含至少一个字段的任何行.删除数据中的S减少了样本量。删除南斯也可以创建不规则的时间序列。

  • 模拟假设您同步了预样本数据,以便对每个预样本系列的最新观察同时发生。

输出参数

全部折叠

平均零新息的模拟条件方差路径Y,返回为数值列向量或矩阵。

V是一个numObs-借-努帕斯矩阵,其中每一列对应一个模拟的条件方差路径。行V周期是否对应于的周期Mdl

模拟响应路径,作为数字列向量或矩阵返回。Y通常表示具有条件方差的创新的均值为零的异方差时间序列V(前样创新系列的延续E0).

Y也可以表示均值为零的异方差创新加上偏移量的时间序列。如果Mdl那么,包括偏移量模拟对潜在的均值为零的异方差创新进行偏移,这样Y表示经过补偿调整的创新的时间序列。

Y是一个numObs-借-努帕斯矩阵,其中每一列对应一个模拟的响应路径。行Y周期是否对应于的周期Mdl

工具书类

[1] 广义自回归条件异方差性计量经济学杂志.1986年第31卷,307-327页。

关于投机性价格和收益率的条件异方差时间序列模型。《经济学与统计评论》1987年第69卷,第542-547页。

[3] 博克斯,G.E.P.,G.M.詹金斯和G.C.莱因塞尔。时间序列分析:预测与控制.3版。恩格尔伍德悬崖,NJ: Prentice Hall, 1994。

恩德斯[4],W。应用计量经济学时间序列.霍博肯:约翰·威利父子公司,1995。

[5] Engle, R. F. <英国通货膨胀方差估计的自回归条件异方差性>。计量经济学1982年第50卷,第987-1007页。

[6] 关于股票名义超额收益率的预期值与波动性之间的关系金融杂志第48卷,第5期,1993年,第1779-1801页。

j·D·汉密尔顿时间序列分析新泽西州普林斯顿:普林斯顿大学出版社,1994年。

[8] Nelson, D. B. <资产收益的条件异方差:一种新方法>。计量经济学.1991年第59卷,第347-370页。

介绍了R2012a