stata马尔可夫Markov区制转移模型分析基金利率

简介: stata马尔可夫Markov区制转移模型分析基金利率

过程会随着时间的推移而发展,结果会发生变化。

考虑一下经济衰退和扩张。在衰退开始时,产出和就业率下降并保持较低水平,然后,产出和就业率增加。从统计上讲,均值,方差和其他参数在各个状态之间都在变化。我们的问题是估计方案何时更改以及与每个方案关联的参数值。询问状态何时改变等同于询问状态持续多久。

在马尔可夫模型中,除了估算每个方案的均值,方差之外,我们还估算区制变化的可能性。某些问题的估计转移概率可能如下:






from/to
state     1 2

1     0.82 0.18

2     0.75 0.25

从状态1开始。从状态1转换为状态1的概率为0.82。换句话说,一旦处于状态1,该过程便会停留在那里。但是,以0.18的概率,过程转换到状态2。状态2的持久性不那么强。在下一个时间段,过程从状态2转换为状态1的概率为0.75。

马尔可夫转换模型不限于两种状态,尽管两种状态模型是常见的。

在上面的示例中,我们将转换描述为突然的变化:概率立即改变。这种马尔可夫模型称为动态模型。马尔可夫模型还可以通过将转移概率建模为自回归过程来拟合更平滑的变化。

因此,转换可以是平稳的或突然的。

基金利率案例

让我们看一下不同状态之间的均值变化。我们分析基金利率,研究1954年至2010年底之间基金利率的变化。以下是数据:

 

我们有季度数据。高利率似乎是七十年代和八十年代的特征。我们将假定还有另一种低利率的状态,这好像是其他几十年的特征。

为了使动态模型具有两种状态



mswit
Performing gradient-based optimization:






Iteration 0:  log likelihood = -508.66031

Iteration 1:  log likelihood = -508.6382

Iteration 2:  log likelihood = -508.63592

Iteration 3:  log likelihood = -508.63592

马尔可夫转换动态回归样本:1954q3-2010q4观测值数量= 226 状态数= 2  AIC = 4,5455 无条件概率:HQIC = 4,5760  SBIC = 4,6211 对数似然= -508.63592




fedfunds     Coef. Std. Err. z P>|z| [95% Conf. Interval]

State1

_cons     3.70877 .1767083 20.99 0.000 3.362428 4.055112

State2

_cons     9.556793 .2999889 31.86 0.000 8.968826 10.14476

sigma     2.107562 .1008692 1.918851 2.314831

p11     .9820939 .0104002 .9450805 .9943119

p21     .0503587 .0268434 .0173432 .1374344

在上面的输出中

  • 两种状态的平均值_cons);
  • 整个过程的单个标准差(sigma);
  • 状态1到1和状态2到1的转移概率(p11p21)。

 

State1是中利率状态(平均值为3.71%)。State2是高利率状态(平均9.56%)。






from/to
state     1 2

1     0.98 1 - 0.98

2     0.05 1 - 0.05

两种状态都是持久性(1-> 1和2-> 2概率分别为0.98和0.95)。

估计后可以预测的包括处于各种状态的概率。我们只有两个状态,因此处于(例如)状态2的概率告诉我们两个状态的概率。我们可以获得预测的概率并将其与原始数据一起绘制成图形:

该模型在每个时间点的状态几乎没有不确定性。我们看到三个时期的高利率状态和四个时期的中利率状态。

疾病案例

让我们看一个疾病的例子,即1929年至1972年之间腮腺炎。您可能会认为疾病对应于均值变化,但是我们在数据中看到的是方差更大的变化:

 

我们绘制了变量S12.mumpspc的图表 ,这意味着在12个月内人均季节性差异性腮腺炎病例,我们将分析 S12.mumpspc

我们将假设两个状态,其中S12.mumpspc的均值和方差会 发生 变化。拟合动态模型



mswit Performing EM optimizaton:
 




Iteration 0:  log likelihood = 110.9372 (not concave)
Iteration 1:  log likelihood = 120.68028

Iteration 2:  log likelihood = 123.23244

Iteration 3:  log likelihood = 131.47084

Iteration 3:  log likelihood = 131.72182

Iteration 3:  log likelihood = 131.7225

Iteration 3:  log likelihood = 131.7225

马尔可夫转换动态回归样本:1929m2-1972m6 obs数量= 521状态数= 2 AIC = -0.4826 无条件概率:HQIC = -0.4634 SBIC = -0.4336 对数似然= 131.7225





mumspc     Coef. Std. Err. z P>|z| [95% Conf. Interval]

State1

mumpspc

LS12.     .4202751 .0167461 25.10 0.000 .3874533 .4530968

State2

mumpspc

LS12.     .9847369 .0258383 38.11 0.000 .9340947 1.035379

sigma1     .0562405 .0050954 .0470901 .067169

sigma2     .2611362 .0111191 .2402278 .2838644

p11     .762733 .0362619 .6846007 .8264175

p12     .1473767 .0257599 .1036675 .205294

报告

  • S12.mumpspc的两个状态的 平均值 (0.42和0.98);
  • 两种状态的标准偏差(0.06和0.26);
  • 状态1到状态1和状态2到状态1的转移概率(0.76和0.15)。

 

状态1是低方差状态。

转移概率如下:



from/to
state     1 2

1     0.76 1 - 0.76

2     0.15 1 - 0.15

与以前的模型一样,状态是持久的。


相关文章
|
JavaScript 前端开发 程序员
前端开发:在vue中实现按钮倒计时功能
在vue中开发中,关于计时器的使用也是比较常见的知识点,如基于移动端的短信60秒倒计时按钮、阅读协议40秒等实际需求,很常见的需求也不复杂,主要是有些细节需要注意。那么本篇博文就来分享一下关于vue中实现倒计时功能的方法,方便查阅使用,如有问题请指正。
1696 0
|
并行计算 TensorFlow 算法框架/工具
【Deepin 20系统】解决Check failed: err == cudaSuccess || err == cudaErrorInvalidValue Unexpected CUDA erro
本文介绍了在使用Nvidia RTX 2070 GPU和TensorFlow 2时,解决GPU内存不足错误的方法,包括杀死占用内存的进程、重置GPU以及重启设备等方案。
267 3
|
8月前
|
弹性计算 资源调度 搜索推荐
阿里云ECS中长期成本节省计划解析:从原理到实战,助力企业降本提效
阿里云ECS节省计划的推出为企业用户提供了一种全新的成本优化方案。通过一次性购买的方式享受长期按量付费的折扣权益,客户不仅可以大幅降低ECS资源的使用成本还可以享受更高的灵活性和便捷性。本文将从多个维度深入剖析阿里云ECS节省计划,包括其核心优势、详尽的购买使用指引、与传统付费模式的全面对比,以及一客户成功案例,以供大家了解和参考。
|
安全 PHP 开发者
Web安全-URL跳转与钓鱼
Web安全-URL跳转与钓鱼
267 8
|
Ubuntu Linux
Linux服务器的自动启动可以在哪里进行配置?
Linux服务器的自动启动可以在哪里进行配置?
417 3
|
存储 关系型数据库 MySQL
MySQL bit类型增加索引后查询结果不正确案例浅析
【8月更文挑战第17天】在MySQL中,`BIT`类型字段在添加索引后可能出现查询结果异常。表现为查询结果与预期不符,如返回错误记录或遗漏部分数据。原因包括索引使用不当、数据存储及比较问题,以及索引创建时未充分考虑`BIT`特性。解决方法涉及正确运用索引、理解`BIT`的存储和比较机制,以及合理创建索引以覆盖各种查询条件。通过`EXPLAIN`分析执行计划可帮助诊断和优化查询。
273 1
|
SQL 分布式计算 大数据
大数据平台的毕业设计01:Hadoop与离线分析
大数据平台的毕业设计01:Hadoop与离线分析
538 0
|
Prometheus 监控 Cloud Native
Prometheus实战篇:Prometheus监控mongodb
Prometheus实战篇:Prometheus监控mongodb
|
机器学习/深度学习 Go 计算机视觉
YOLOv5改进 | Neck篇 | 利用ASF-YOLO改进特征融合层(适用于分割和目标检测)
YOLOv5改进 | Neck篇 | 利用ASF-YOLO改进特征融合层(适用于分割和目标检测)
819 0
|
JavaScript 内存技术
nvm 切换版本失效?这个方法100%解决
NVM是一个工具,类似于一个“开关”,可以帮助我们在同一台电脑上安装和切换不同版本的Node.js
8038 1
nvm 切换版本失效?这个方法100%解决