3

これは、について学んだ昨日の質問の延長ですapply.weekly。zooこれはうまく機能しますが、幅の広いオブジェクトに対してこれを実行したいと考えています。apply.weeklyワイドで使用するzooと、列が合計され、週ごとの集計が実行されます。

> library(xts)
> set.seed(2001)
> zoo.daily <- zoo(data.frame(a=rnorm(20), b=rnorm(20), c=rnorm(20)), order.by=as.Date("2001-05-25") + 0:19)
> apply.weekly(zoo.daily, sum)
2001-05-27 2001-06-03 2001-06-10 2001-06-13 
  1.091999  -3.017688   3.842305   2.045370 
> apply.weekly(zoo.daily[, 1] + zoo.daily[, 2] + zoo.daily[, 3], sum) 
2001-05-27 2001-06-03 2001-06-10 2001-06-13 
  1.091999  -3.017688   3.842305   2.045370

apply一連の演算子を試しましたが、それらはzoo日付インデックスを取り除いているようです。ループで実行できますが、実際には時間がかかります (周期性の関数forよりも 4 倍以上遅くなります)。ループは次のとおりです。aggregateas.yearmonfor

week.ends <- index(zoo.daily[endpoints(zoo.daily, "weeks")[-1], ])
num.weeks <- nweeks(zoo.daily)
num.stocks <- ncol(zoo.daily)
zoo.weeks <- zoo(matrix(NA, num.weeks, num.stocks), order.by=week.ends)
for (i in seq(num.stocks)) {
    zoo.weeks[, i] <- apply.weekly(zoo.daily[, i], mean)
}

どちらが機能しますか（つまり、各ベクトルを別々に保ちます）：

2001-05-27 -0.36663040 -0.108648725  0.8392788
2001-06-03  0.33032998  0.003025018 -0.7644534
2001-06-10  0.07816992  0.620198931 -0.1494681
2001-06-13  0.02114608  0.956226189 -0.2955824

を使用してすべての列をすばやく操作する方法はありますapply.weeklyか? ありがとう！

colMeans更新: Joshua Ulrich は、列を認識する関数 (または など)が必要だと指摘していますcolSums。これを行うと、正しい答えが得られますが、転置行列として得られます。再クラスして先に進むべきですか？または、オプション/設定が間違っていますか?

> apply.weekly(zoo.daily, colSums)
        [,1]        [,2]       [,3]        [,4]
a -1.0998912  2.31230989  0.5471894  0.06343824
b -0.3259462  0.02117512  4.3413925  2.86867857
c  2.5178365 -5.35117351 -1.0462765 -0.88674717
4

1 に答える 1

5

で列対応関数を使用する必要がありますapply.weekly。たとえば、colSums代わりにsumまたはcolMeansの代わりに使用しmeanます。

xtson R-forgeの最近のリビジョンでは、以下の出力が得られます。現在 CRAN にあるバージョンは、転置されたデータを返します。

# install.packages("xts", repos="http://r-forge.r-project.org")
> apply.weekly(zoo.daily, colSums)
                     a           b          c
2001-05-27 -1.09989120 -0.32594617  2.5178365
2001-06-03  2.31230989  0.02117512 -5.3511735
2001-06-10  0.54718941  4.34139252 -1.0462765
2001-06-13  0.06343824  2.86867857 -0.8867472
> apply.weekly(zoo.daily, colMeans)
                     a            b          c
2001-05-27 -0.36663040 -0.108648725  0.8392788
2001-06-03  0.33032998  0.003025018 -0.7644534
2001-06-10  0.07816992  0.620198931 -0.1494681
2001-06-13  0.02114608  0.956226189 -0.2955824

カスタム関数を使用する必要がある場合は、 と を組み合わせて使用​​できapply.weeklyますapply。

> apply.weekly(zoo.daily, function(x) apply(x,2,mean))
                     a            b          c
2001-05-27 -0.36663040 -0.108648725  0.8392788
2001-06-03  0.33032998  0.003025018 -0.7644534
2001-06-10  0.07816992  0.620198931 -0.1494681
2001-06-13  0.02114608  0.956226189 -0.2955824
于 2010-11-30T21:16:39.043 に答える