0

私は358行のデータセットを持っています。簡単に見ると次のようになります

>        unempts gdpts
> [1,]    7.03   4.2
> [2,]    7.17   3.7
> [3,]    6.97   3.2
> [4,]    6.83   2.9
> [5,]    6.60   2.7
> [6,]    6.27   3.4

列 1 の回帰を列 2 で実行し、60 のローリング ウィンドウですべての結果をマトリックスに保存して、同じ結果をプロットして実行できるようにする必要がありますHAC。

私は次のことを試しました：

regRO <- rollapply(datam, 60, lm, by.column = TRUE,fill=NA)

しかし、次のようにエラーが発生します

formula.default(object, env = baseenv()) のエラー: 無効な数式

何を変更する必要があるかアドバイスしてください。

4

1 に答える 1

1

Xの入力要件を満たすために、最初にデータ フレームに変換されるように、関数呼び出しを明示的に指定する必要がありlmます。

library(zoo)

set.seed(123) # ensure runif in next line uses same random numbers each time run
datam <-matrix(runif(200),ncol=2)

out <- rollapplyr(datam,
          width = 60,
          FUN = function(X) coef(lm(as.data.frame(X))),
          by.column=FALSE,
          fill = NA) # omit fill=NA if NA rows not wanted

与える:

> tail(out)
       (Intercept)         V2
[95,]    0.5071618 -0.1041905
[96,]    0.5015140 -0.1032419
[97,]    0.5196683 -0.1404855
[98,]    0.5171599 -0.1401534
[99,]    0.5122871 -0.1288700
[100,]   0.5176396 -0.1297019
于 2015-11-14T15:38:49.457 に答える