主要内容

模拟

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

描述

例子

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

例子

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

例子

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

例子

全部折叠

用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(模拟的条件方差的)子plot(2,1,2) plot(Y) title(“模拟反应”

图中包含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(模拟的条件方差的)子plot(2,1,2) plot(Y) title(的模拟反应(创新)

图中包含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(模拟的条件方差的)子plot(2,1,2) plot(Y) title(的模拟反应(创新)

图中包含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)拟合的预测结果。

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

负载Data_EquityIdx纳斯达克= 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 {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) =情节(datesFH, mn {j},“颜色”上校,“线宽”3);h(j,3:4) = plot([datesFH datesFH],[lower{j} upper{j}],“:”...“颜色”上校,“线宽”2);结束甘氨胆酸hGCA =;情节(datesFH (1) * [1], hGCA。YLim,“k——”);datetick;轴;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..

输入参数

全部折叠

不含任何未知参数的条件方差模型,指定为garchegarch,或gjr模型对象。

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

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

名称-值对的观点

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

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

要生成的样本路径的数量,指定为逗号分隔的对“NumPaths”一个正整数。VYNumPaths列。

例子:“NumPaths”,1000年

数据类型:

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

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

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

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

  • 如果E0是一个矩阵,然后每一列代表潜在创新系列的前样例路径。E0必须至少NumPaths列。如果E0列比需要的多,模拟使用第一个NumPaths只列。

默认值是:

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

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

例子:“E0”,[0.5;0.5)

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

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

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

  • GARCH (P)和GJR (P)模型:

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

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

  • EGARCH (P)模型,模拟

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

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

如果行数超过了必要的数量模拟仅使用最新的、所需的观测次数。最后一个元素或行包含最新的观察结果。

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

数据类型:

笔记

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

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

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

输出参数

全部折叠

模拟与之相关的均值-零创新的条件方差路径Y,返回为数值列向量或矩阵。

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

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

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

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

参考文献

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

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

[3] Box, G. E. P. G. M. Jenkins和G. C. Reinsel。时间序列分析:预测与控制.3版。恩格尔伍德悬崖,NJ: Prentice Hall, 1994。

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

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

[6] Glosten, L. R., R. Jagannathan, D. E. Runkle。“关于期望值与股票名义超额收益波动的关系”。金融杂志.第48卷,第5期,1993年,1779-1801页。

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

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

介绍了R2012a