3

時系列クラスをプログラムしたいと思います。アイデアは、たとえば、式と他の時系列オブジェクトを使用してオブジェクトをインスタンス化することです。

(2 つの時系列) 

x <- ts(rnorm(10), frequency = 4, start = c(1959, 2))  
y <- ts(rnorm(10), frequency = 4, start = c(1959, 2))

(x と y の合計として定義される時系列) 

z <- exprTs("x+y", parents=list(x=x, y=y))

（シリーズの一部を入手） 

window(z, start=1960, end=1960.75)

問題は、式をどのように評価できるかです。私は次のことを試しました：

#(constructor for class)  
exprTs <- function(expr, parents) {  
  res = list(expr=expr, parents=parents)  
  class(res) <- "exprTs"  
  res  
}  

#(window method)  
window.exprTs <- function(z, ...) {  
  eval(substitute(z$expr, lapply(z$parents, window, ...)))  
  #do.call(z$expr, lapply(z$parents, window, ...))  
}

ウィンドウメソッドが機能しません。

代替、eval、do.call を適切に使用する方法を教えていただければ、非常に役に立ちます。

4

3 に答える 3

6

あなたはこれについて間違った方法で考えていて、物事を複雑にしすぎているかもしれません。tsオブジェクトにはすでに追加が定義されています。

R> set.seed(42)
R> x <- ts(rnorm(10), frequency = 4, start = c(1959, 2))
R> y <- ts(rnorm(10), frequency = 4, start = c(1959, 2))
R> z <- x + y
R> cbind(x,y,z)
               x       y       z
1959 Q2  1.37096  1.3049  2.6758
1959 Q3 -0.56470  2.2866  1.7219
1959 Q4  0.36313 -1.3889 -1.0257
1960 Q1  0.63286 -0.2788  0.3541
1960 Q2  0.40427 -0.1333  0.2709
1960 Q3 -0.10612  0.6360  0.5298
1960 Q4  1.51152 -0.2843  1.2273
1961 Q1 -0.09466 -2.6565 -2.7511
1961 Q2  2.01842 -2.4405 -0.4220
1961 Q3 -0.06271  1.3201  1.2574
R>

Rオブジェクトを操作するために、式パーサーは実際には必要ありません。

于 2010-01-22T01:49:07.933 に答える
5

1つの問題は、式を文字列として指定しているため、文字列に評価されることです。parse文字列を式に解析する場合は、次のコマンドを使用する必要があります。

> "x+y"
[1] "x+y"
> parse(text="x+y")
expression(x+y)
attr(,"srcfile")
<text>

しかし、ええ、apply代わりに関数を使用してみませんか？

于 2010-01-22T01:11:51.810 に答える
1

今、私は解決策を見つけました：

window.exprTs <- function(z, ...) {  
  names(z$parents)<-NULL  
  do.call(z$expr, lapply(z$parents, window, ...))  
}  

plus <- function(x, ...) if (nargs() == 1) x else x + Recall(...)  
z <- exprTs(plus, parents=list(x=x, y=y))

今 

window(z, start=1960, end=1960.75)

望ましい結果が得られます。

于 2010-01-22T20:44:38.883 に答える