【问题标题】:Error in Markov-switching VAR in RR中的马尔可夫切换VAR错误
【发布时间】:2016-03-02 09:43:18
【问题描述】:

我正在尝试使用命令 msvar 估计 R 中的 Markov-switching VAR。这是我的两个时间序列的前 10 个条目。我有 798。当我尝试运行它时,我收到一条错误消息

a <- c(1.998513, 1.995302, 2.030693, 2.122130, 2.236770, 2.314639, 2.365214, 2.455784, 2.530696, 2.596537)
b <- c(0.6421369, 0.6341437, 0.6494933, 0.6760939, 0.7113511, 0.7173038, 0.7250545, 0.7812490, 0.7874657, 0.8275209)
x <- matrix (NA,10,2)
x[,1] <- a
x[,2] <- b
time.seriesx <- ts(x)
markov.switchingx <- msvar(time.seriesx, p = 2, h = 2, niterblkopt = 10)

我得到的错误信息如下:

优化错误(par = c(beta0.it), fn = llf.msar, Y = Yregmat, X = Xregmat, : 'vmmin' 中的初始值不是有限的

谁能帮帮我?谢谢

【问题讨论】:

    标签: r time-series var markov


    【解决方案1】:

    我认为你必须先运行 log-likehood 函数。我得到了同样的错误,但是当我这样做时,它可以工作。 我不确定,但我希望这可以帮助你:(我使用了我的数据,所以不要关注“M1euro”)

    library(base)
    
    data <- data.matrix(M1euro, rownames.force = NA)
    
    library(stats)
    
    ss1<-ts(data, frequency=12, start=c(2007,1), end=c(2016,4))
    
    class(ss1)
    
    length(ss1)
    
    ss <- na.approx(ss1,na.rm=F,rule=2)
    
    ss
    
    class(ss)
    
    library(MSBVAR)
    
    require(graphics)
    
    set.seed(1)
    

    【讨论】:

      猜你喜欢
      • 1970-01-01
      • 2015-11-15
      • 2016-12-14
      • 1970-01-01
      • 1970-01-01
      • 1970-01-01
      • 2017-05-15
      • 1970-01-01
      相关资源
      最近更新 更多