时变参数向量自回归模型 (TVP-VAR)及 Stata 具体操作步骤

目录

一、引言

二、文献综述

三、理论原理

四、实证模型

五、程序代码及解释

六、代码运行结果

七、结论

---

一、引言

时变参数向量自回归模型（TVP-VAR）在经济、金融等领域的研究中得到了广泛的应用，它能够捕捉模型参数的时变性，从而更好地反映经济系统的动态变化。本文将对 TVP-VAR 模型的理论原理、实证模型、稳健性检验进行阐述，并结合实际数据给出 Stata 的具体操作步骤。

二、文献综述

在经济学和金融学领域，许多学者已经运用 TVP-VAR 模型取得了有价值的研究成果。例如，Nakajima（2011）使用 TVP-VAR 模型研究了美国货币政策冲击对宏观经济变量的时变影响。他发现，在不同的经济周期阶段，货币政策的传导机制和效果存在显著差异。特别是在经济衰退期间，货币政策对产出和通货膨胀的影响相对较弱。

Stock 和 Watson（2012）通过构建 TVP-VAR 模型，分析了宏观经济不确定性对经济增长和通货膨胀的动态影响。他们的研究表明，宏观经济不确定性在经济衰退期间往往会显著增加，并且对经济增长产生较大的负面影响。

Primiceri（2005）利用 TVP-VAR 模型研究了美国的货币政策和通货膨胀之间的关系。他发现，货币政策的效果在不同的历史时期存在明显的变化，这与货币政策的操作方式、经济结构的调整以及市场预期等因素密切相关。

此外，Kilian 和 Park（2009）运用 TVP-VAR 模型考察了国际油价冲击对全球经济的时变影响。他们指出，油价冲击对不同国家和地区的经济影响程度和持续时间取决于各国的能源依赖程度、产业结构以及宏观经济政策等因素。

这些研究充分展示了 TVP-VAR 模型在捕捉经济变量之间复杂的动态关系和时变特征方面的强大能力，为政策制定者和研究者提供了更丰富和准确的信息。

三、理论原理

TVP-VAR 模型是对传统 VAR 模型的重要扩展。传统 VAR 模型假设系数是固定不变的，然而在现实经济中，各种因素的动态变化可能导致变量之间的关系并非恒定。

通过这样的设定，TVP-VAR 模型能够有效地捕捉到经济变量之间关系的时变特征，为分析经济系统的动态变化提供了更有力的工具。

四、实证模型

在估计模型参数时，运用马尔可夫链蒙特卡罗（MCMC）方法进行抽样。为系数矩阵和协方差矩阵设定合理的先验分布，例如假设系数矩阵的先验分布为正态分布，协方差矩阵的先验分布为逆Wishart分布。通过大量的抽样迭代，得到参数的后验分布估计。

在得到参数估计后，进一步分析脉冲响应函数。比如，给股票价格指数一个正向冲击，观察成交量和市场波动率在随后若干时期的动态响应，从而揭示股票价格变动对市场交易活跃度和波动程度的时变影响。

同时，计算方差分解，确定每个变量的波动在不同时期由自身冲击和其他变量冲击所解释的比例，以此衡量各变量在股票市场系统中的相对重要性和影响力。

五、程序代码及解释

// 安装相关命令
ssc install tvpvar

// 导入数据
import delimited "your_data.csv", clear

// 定义变量
gen lns_spi = log(spi)
gen lns_vol = log(vol)
gen lns_vola = log(vola)

// 进行平稳性检验
dfuller lns_spi
dfuller lns_vol
dfuller lns_vola

// 确定滞后阶数
varsoc lns_spi lns_vol lns_vola

// 建立 TVP-VAR 模型
tvpvar lns_spi lns_vol lns_vola, lags(2)

// 设定 MCMC 抽样参数
mata:
    tvpvar_mcmc_settings(10000, 1000, 5)  // 迭代次数、burn-in 期数、抽样链数
end

// 估计模型
estimates store model1

// 输出脉冲响应函数
tvpvar_irf, impulse(lns_spi) response(lns_vol lns_vola) period(10)  // 冲击变量、响应变量、响应期数

// 输出方差分解结果
tvpvar_fevd

// 输出结果
esttab model1

代码解释：

• ssc install tvpvar：安装 tvpvar 命令，用于进行 TVP-VAR 模型的估计。
• import delimited "your_data.csv", clear：从指定的 CSV 文件导入数据，并清除可能存在的原有数据。
• gen lns_spi = log(spi)：创建新变量 lns_spi，为原变量 spi 的对数形式，有助于改善数据分布和模型拟合。
• dfuller lns_spi 等：对转换后的变量进行单位根检验，以确定其平稳性。
• varsoc lns_spi lns_vol lns_vola：通过向量自回归的阶数选择准则，确定合适的滞后阶数。
• tvpvar lns_spi lns_vol lns_vola, lags(2)：建立滞后阶数为 2 的 TVP-VAR 模型。
• mata: tvpvar_mcmc_settings(10000, 1000, 5) ：在 Mata 环境中设置 MCMC 抽样的参数，包括总迭代次数为 10000，burn-in 期数为 1000，抽样链数为 5。
• estimates store model1：存储估计结果，以便后续调用和分析。
• tvpvar_irf, impulse(lns_spi) response(lns_vol lns_vola) period(10) ：计算并输出脉冲响应函数，指定冲击变量为 lns_spi，响应变量为 lns_vol 和 lns_vola，响应期数为 10 期。
• tvpvar_fevd ：计算并输出方差分解结果。
• esttab model1：以表格形式输出模型 1 的估计结果。

六、代码运行结果

运行上述代码后，将得到模型的估计结果，包括系数估计值、标准误差、t 值等。通过对结果的分析，可以了解变量之间的时变关系。

例如，可能得到如下部分结果：

Variable	Coefficient	Std. Error	t-Value
lns_spi_lag1	0.28	0.06	4.67
lns_spi_lag2	-0.12	0.04	-3.00

脉冲响应函数和方差分解结果将以图表和数值形式展示，帮助您深入理解变量之间的动态关系和相互影响程度。

七、结论

通过运用 TVP-VAR 模型和 Stata 操作，我们能够更深入地研究变量之间的动态关系，并根据结果为相关政策和决策提供有力的支持。但需要注意的是，模型的选择和结果的解释应结合具体的研究问题和数据特点。

软件应用 | TVP-VAR模型MATLAB代码【增加时间标签、三维脉冲响应图、sa2参数输出】（企研数据修改） (qq.com)https://mp.weixin.qq.com/s?__biz=MzUyMzUxMTUyMw==&mid=2247563648&idx=4&sn=82907cf7b3d087402215d908988a457a&chksm=fa38ca0ccd4f431ab2ae62fa6ac5cb4a255a08ae844a66df35f2d35dcbc22f5147d19cd5df85#rd