R语言VAR模型平稳性检验求助！

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

dat<-read.csv(file="rawdata.csv")

var.mod1<-VAR(dat[,2:4],p=1,type="c")
var.mod2<-VAR(dat[,2:4],p=2,type="c")
var.mod3<-VAR(dat[,2:4],p=3,type="c")
var.mod4<-VAR(dat[,2:4],p=4,type="c")

AIC(var.mod1,var.mod2,var.mod3,var.mod4)

BIC(var.mod1,var.mod2,var.mod3,var.mod4)

var.mod<-var.mod3
summary(var.mod)

Rcode如上，现在就想请教各位大神如何检验我的var模型的平稳性，自己试了ur.df等都返回Error messages，该怎么办啊T T求助求助，不知道这个版还有没有好心人能帮一下T T

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏2 回帖

关键词：平稳性检验 VAR模型 AR模型 平稳性 VaR 模型

自己无限顶贴，一定会有大神帮忙的！

再顶一个

我自己解决了。。。

我遇到了相同的问题，请问您是怎么处理的？

我遇到了相同的问题，请问您是怎么处理的？

如果是打算绘制像Eviews或者stata里的单位圆来判断VAR系统的稳定性，可以使用vars包里的roots()函数计算出伴随矩阵的特征值，就可以进行判断了。如果打算绘制单位圆，可使用如下的自编函数实现。

### VAR系统稳定性判别图的绘制
##compute the unit roots of VAR model
varstability <- function(x) {
    # Calculate the eigenvalues & Modulus of the VAR model.
    Eigenvalue<-roots(x, modulus = FALSE)
    Modulus<-roots(x, modulus = TRUE)
    res <- data.frame(Eigenvalue, Modulus)
    class(res) <- c("varstability", "data.frame")
    return(res)
}


## plot VAR stability
plot.varstability <- function(x, ...) {
    x1 <- seq(-1,1,length=1000)
    y1 <- sqrt(1 - x1^2)
    y2 <- -sqrt(1 - x1^2)
   
    dfroots <- data.frame(Real=Re(x$Eigenvalue),
                          Imaginary=Im(x$Eigenvalue))
   
    limits <- max(abs(dfroots))
   
    if (limits < 1) {
        limits <- c(-1,1)
    } else {
        limits <- c(-ceiling(limits), ceiling(limits))
    }
   
    ggplot2::ggplot(data.frame(x=x1,y1=y1,y2=y2)) + #Using the results from earlier
        ggplot2::geom_hline(ggplot2::aes(yintercept=0)) +
        ggplot2::geom_vline(ggplot2::aes(xintercept=0)) +
        ggplot2::geom_line(ggplot2::aes(x=x,y=y1)) + # top of circle
        ggplot2::geom_line(ggplot2::aes(x=x,y=y2)) +
        ggplot2::geom_point(ggplot2::aes_(x=~Real, y=~Imaginary),
                            dfroots, size = 2,col="blue") +
        ggplot2::scale_y_continuous("Imaginary", limits=c(-1,1)) +
        ggplot2::scale_x_continuous("Real", limits=c(-1,1)) + ggplot2::ggtitle("Roots of the companion matrix")
   
}

var_stab<-varstability(out.var2)
plot(var_stab)


可以绘制出单位圆如下，其中out.vars是使用VAR()函数生成的VAR对象。