2017-06-01 20 views
1

ヘルプ解決:Rラウンドでエラーが発生しました(周波数):VARモデルから数学関数に非数値引数

fcst.tmp <を - 予測ラウンドで(VAR.tmp) エラー(周波数):数学関数への非数値引数

コードから:

library(vars) 

tmp.matrix <- matrix(data = c(0,      -0.0018674142,      0.0046620131, 
        0.0037140247,      -0.0009272138,      0.0064725145, 
        -0.000922084,      -0.0046232167,      0.0210920245, 
        -0.0100320039,      0.0136551962,      -0.00725956, 
        0.0251811866,      -0.0089206658,      -0.000896459, 
        0.0035810244,      0.0159577854,      -0.0097216851, 
        0.0176060886,      0.0034843241,      -0.0017406445, 
        0.0198367867,      -0.0034217313,      -0.0225312824, 
        -0.0043917506,      0.0157208478,      -0.0148409533, 
        0.0234691345,      -0.0401374128,      -0.0132217389, 
        0.0269124178,      -0.006590583,      0.0255472708, 
        0.003401652,      0.0052493174,      -0.0011235182, 
        -0.0033622663,      -0.0193530068,      -0.0302611445, 
        -0.0076587663,      -0.0023782155,      -0.026358019, 
        0.0161200327,      0.0171822796,      0.0193783756, 
        0.0091528136,      -0.0101419116,      -0.0013561094, 
        0.000172569,      -0.0071663343,      -0.0119922952, 
        0.0044839696,      0.0029657989,      -0.0114150402, 
        -0.0018366463,      0.0116721141),ncol = 2) 

colnames(tmp.matrix) <- c("Col1","Col2") 

VAR.tmp <- VAR(tmp.matrix, p=1, type = "const") 

fcst.tmp <- forecast(VAR.tmp) 

PS

tmp.matrixニーズはしかし、時系列に変換される:私は同じ問題を抱えていた

ts.matrix <- as.ts(tmp.matrix) 

VAR.tmp <- VAR(ts.matrix, p=1, type = "const") 

fcst.tmp <- forecast(VAR.tmp) 

plot(fcst.tmp) 
+0

最後に、エラーは再現できません。私は 'forecast'関数を使って確かめてみました。そして' VAR.tmp'の細かいプロットを取得しました。 – Masoud

答えて

2

tmp.matrixは通過時系列に変換する必要があります。これは、forecast()の代わりにpredict()を使用すると解決できます。

ts.matrix <- as.ts(tmp.matrix) 

VAR.tmp <- VAR(ts.matrix, p=1, type = "const") 

fcst.tmp <- predict(VAR.tmp) 

plot(fcst.tmp) 
+0

これは問題を解決しませんでした – mike

0

ts.matrix <- as.ts(tmp.matrix) 

VAR.tmp <- VAR(ts.matrix, p=1, type = "const") 

fcst.tmp <- forecast(VAR.tmp) 

plot(fcst.tmp)