R语言从经济时间序列中用HP滤波器,小波滤波和经验模态分解等提取周期性成分分析|附代码数据

发布时间 2023-06-27 23:49:49作者: 拓端tecdat

全文下载链接:http://tecdat.cn/?p=9350

最近我们被客户要求撰写关于经济时间序列的研究报告,包括一些图形和统计输出。

经济时间序列的分析通常需要提取其周期性成分。这篇文章介绍了一些方法,可用于将时间序列分解为它们的不同部分 ( 点击文末“阅读原文”获取完整代码数据******** ) 。

它基于《宏观经济学手册》中Stock和Watson(1999)关于商业周期的章节,但也介绍了一些较新的方法,例如汉密尔顿(2018)替代HP滤波器,小波滤波和_经验模态分解_。

数据

我使用从1970Q1到2016Q4的美国对数实际GDP的季度数据来说明不同的方法。时间序列是通过  Quandl  及其相应的R包获得的。

 
#加载用于数据下载和转换的软件包
library(dplyr)

#下载数据
"FRED/GDPC1", order = "asc",
               start_date = "1970-01-01", end_date = "2016-10-01")  %>%
  rename(date = Date,
         gdp = Value) %>%
  mutate(lgdp = log(gdp)) # 获取对数
 
library(ggplot2)

ggplot(data, aes(x = date, y = lgdp)) +
  geom_line() +
  theme_classic()

图片

数据有明显的增长趋势,到现在似乎逐渐变小。此外,似乎或多或少有规律地围绕这一趋势波动。与趋势之间存在相对较长的持久偏差,可以将其视为周期性波动。


点击标题查阅往期内容

图片

R语言提取时间序列的周期性成分应用EMD,小波滤波器,Baxter过滤器等

图片

左右滑动查看更多

图片

01

图片

02

图片

03

图片

04

图片

与线性趋势的偏差

从系列中提取趋势的第一种方法是在常数和趋势项上回归目标变量并获得拟合值。在下图中绘制。

 
# 添加趋势
data <- data %>%
  mutate(trend = 1:n())

# 用常数和趋势估算模型
time_detrend <- fitted(lm(lgdp ~ trend, data = data))
names(time_detrend) <- NULL

# 将系列添加到主数据框
data <- data %>%
  mutate(lin_trend = time_detrend)

# 为图创建数据框
temp <- data %>%
  select(date, lgdp, lin_trend) %>%
  gather(key = "Variable", value = "value", -date)

# 画图
ggplot(temp, aes(x = date, y = value, colour = Variable)) +

图片

此方法相对有争议,因为它假定存在恒定的线性时间趋势。正如我们在上面看到的,鉴于趋势的增长率随着时间的推移持续下降,这不太可能。但是,仍然可以采用时间趋势的其他函数形式(例如二次项)来说明趋势的特殊性。该方法的另一个缺点是,它仅排除趋势,而不排除噪声,即序列中很小的波动。

Hodrick-Prescott过滤器

Hodrick和Prescott(1981)开发了一个过滤器,将时间序列分为趋势和周期性分量。与线性趋势相反,所谓的  HP过滤器可  估算趋势,该趋势会随时间变化。研究人员手动确定允许这种趋势改变的程度,即平滑参数λλ。

文献表明季度数据的值为1600。但是,也可以选择更高的值。下图绘制了由HP过滤器获得的实际GDP周期性成分的值,并将其与线性趋势下的序列的值进行比较。

图片

尽管HP过滤器在经济学中得到了广泛的应用,但它们的某些功能也受到了广泛的批评。

基于回归的HP过滤器

汉密尔顿(2018)还提出了另一种HP过滤器的方法。它可以归结为一个简单的回归模型,其中 时间序列的_第_ h _个_前导根据时间序列的最新p值进行回归。

图片

Baxter King过滤器

Baxter和King(1994,1999)提出了一种过滤器,其产生的结果与HP过滤器非常相似。另外,它从时间序列中去除了噪声,因此可以对周期分量进行平滑估计。该方法的一个相对严重的缺点是,平滑因子导致序列开始和结束时观测值的损失。当样本量较小且当前经济状况令人关注时,这可能是一个问题。

图片

小波滤波器

Yogo(2008)提出使用小波滤波器从时间序列数据中提取业务周期。该方法的优点是该函数不仅允许提取序列的趋势,周期和噪声,而且还可以更明确地了解周期发生的时间段。

R中的方法实现也很简洁,但是在使用之前需要进行一些其他的数据转换。

 
 # 计算对数GDP的一阶差分
data <- data %>%
  mutate(dlgdp = lgdp - lag(lgdp, 1))

#获取数据
y <- na.omit(data$dlgdp)

#运行过滤器
wave_gdp <- mra(y, J = 5)

# 创建用于绘制的数据框
temp <- wave_gdp %>%
  gather(key = "imf", 

# 绘制mra输出
ggplot(temp, aes(x = date, y = value)) +
  geom_line() +

图片

 
data <- data %>%
  select(date, bk, wave) %>%
  gather(key = "Variabl

ggplot(temp, aes(x = date, y = value, colour = Variable)) +
  geom_hlin

图片

经验模态分解(EMD)

Kozic和Sever(2014)提出了经验模态分解作为商业周期提取的另一种方法,正如Huang等人(2014年)提出的那样。(1998)。 emd 函数可以在_EMD_  包中找到,  并且需要一个不同的时间序列,一个边界条件和一个指定的规则,在该点上迭代算法可以停止。滤波方法的结果与HP,BK和小波滤波相对不同。

图片

图片

Chan(2017)

初始值

 
# X_gamma
x_gamma <- cbind(2:(tt +


# H_2
h2 <- diag(1, tt)
diag(h2[-1,  
t)]) <- 1
h2h2 <- crossprod(h2)

# H_phi
h_phi <- diag(1, tt)
phi <- matrix(


# sigma tau的逆
s_tau_i <- 1 / .001

#  sigma c的逆
s_c_i <- 1 / .5

# gamma
gamma <- t(rep(y[1], 2)) # 应该接近该序列的第一个值

Gibbs 采样

图片

图片

点击文末 “阅读原文”

获取全文完整代码数据资料。

本文选自《R语言从经济时间序列中用HP滤波器,小波滤波和经验模态分解等提取周期性成分分析》。

点击标题查阅往期内容

R语言计量经济学:工具变量法(两阶段最小二乘法2SLS)线性模型分析人均食品消费时间序列数据和回归诊断
R语言用标准最小二乘OLS,广义相加模型GAM ,样条函数进行逻辑回归LOGISTIC分类
R语言中实现广义相加模型GAM和普通最小二乘(OLS)回归
R语言中的偏最小二乘回归PLS-DA
R语言中的偏最小二乘PLS回归算法
偏最小二乘回归(PLSR)和主成分回归(PCR)
R语言用线性模型进行预测:加权泊松回归,普通最小二乘,加权负二项式模型,多重插补缺失值R语言计量经济学:虚拟变量(哑变量)在线性回归模型中的应用
PYTHON用时变马尔可夫区制转换(MRS)自回归模型分析经济时间序列
R语言使用ARIMAX预测失业率经济时间序列数据
R语言经济学:动态模型平均(DMA)、动态模型选择(DMS)预测原油时间序列价格
R语言用向量自回归(VAR)进行经济数据脉冲响应研究分析
R语言从经济时间序列中用HP滤波器,小波滤波和经验模态分解等提取周期性成分分析
R语言计量经济学与有时间序列模式的机器学习预测
R语言使用灰色关联分析(Grey Relation Analysis,GRA)中国经济社会发展指标
R语言用向量自回归(VAR)进行经济数据脉冲响应研究分析
R语言从经济时间序列中用小波滤波和经验模式分解等提取周期性成分分析
R语言用向量自回归(VAR)进行经济数据脉冲响应研究分析
R语言实现偏最小二乘回归法 partial least squares (PLS)回归
R语言中的block Gibbs吉布斯采样贝叶斯多元线性回归
R语言Lasso回归模型变量选择和糖尿病发展预测模型
R语言实现贝叶斯分位数回归、lasso和自适应lasso贝叶斯分位数回归分析
Python贝叶斯回归分析住房负担能力数据集
Python用PyMC3实现贝叶斯线性回归模型