2012-03-01 8 views
9

xts行列のすべての期間に関数を実行したいとします。適用()非常に高速ですが、返される行列は、元のオブジェクトに比べて寸法を移調していますapply()はなぜ転置されたxts行列を返しますか?

> dim(myxts) 
[1] 7429 48 
> myxts.2 = apply(myxts, 1 , function(x) { return(x) }) 
> dim(myxts.2) 
[1] 48 7429 
> str(myxts) 
An 'xts' object from 2012-01-03 09:30:00 to 2012-01-30 16:00:00 containing: 
    Data: num [1:7429, 1:48] 4092500 4098500 4091500 4090300 4095200 ... 
- attr(*, "dimnames")=List of 2 
    ..$ : NULL 
    ..$ : chr [1:48] "Open" "High" "Low" "Close" ... 
    Indexed by objects of class: [POSIXlt,POSIXt] TZ: 
    xts Attributes: 
NULL 
> str(myxts.2) 
num [1:48, 1:7429] 4092500 4098500 4091100 4098500 0 ... 
- attr(*, "dimnames")=List of 2 
    ..$ : chr [1:48] "Open" "High" "Low" "Close" ... 
    ..$ : chr [1:7429] "2012-01-03 09:30:00" "2012-01-03 09:31:00" "2012-01-03 09:32:00" "2012-01-03 09:33:00" ... 
> nrow(myxts) 
[1] 7429 
> head(myxts) 
         Open High  Low Close 
2012-01-03 09:30:00 4092500 4098500 4091100 4098500 
2012-01-03 09:31:00 4098500 4099500 4092000 4092000 
2012-01-03 09:32:00 4091500 4095000 4090000 4090200 
2012-01-03 09:33:00 4090300 4096400 4090300 4094900 
2012-01-03 09:34:00 4095200 4100000 4095200 4099900 
2012-01-03 09:35:00 4100000 4100000 4096500 4097500 

どのように私はmyxts寸法を維持することができますか?

+1

小さなデータが役に立つので、問題を再現することができます。 – Justin

答えて

15

これは、文書番号applyです。 ?applyから:

値:あなたのケースでは

If each call to ‘FUN’ returns a vector of length ‘n’, then ‘apply’ 
returns an array of dimension ‘c(n, dim(X)[MARGIN])’ if ‘n > 1’. 

'n'=48(あなたが行をループしているので)、そうapplyは次元c(48, 7429)の配列を返します。

また、myxts.2ではなく、xtsオブジェクトです。それは規則的な配列です。あなたはカップルのオプションがあります。

  1. はあなたのXTSのオブジェクトを再作成する前に、applyの結果をトランスポーズ:

    data(sample_matrix) 
    myxts <- as.xts(sample_matrix) 
    dim(myxts) # [1] 180 4 
    myxts.2 <- apply(myxts, 1 , identity) 
    dim(myxts.2) # [1] 4 180 
    myxts.2 <- xts(t(apply(myxts, 1 , identity)), index(myxts)) 
    dim(myxts.2) # [1] 180 4 
    
  2. それは、すべての行に作用するように、あなたの関数をベクトル化すること オブジェクトをXTSとxtsオブジェクトを返します。その後、あなたは約約applyについて心配する必要はありません。

最後に、再現可能な例を提供してください。それはそれほど難しいことではなく、人々の手助けをより簡単にします。私は上記の例を提供しており、あなたが次の質問にそれを使用できることを願っています。

+0

オプション2に行きたい場合、どのように各行を自分の処理機能に渡すことができますか?たとえば、myxts $ result = myxts $ Close - myxts $は特定のフィールドで機能しますが、フル行を渡す関数をパラメータの1つとして呼び出す場合はどうすればよいですか? –

+1

@RobertKubrick:あなたの機能を見ることなく、私はそれに答えることができません。それは可能かもしれないし、そうでないかもしれない。少なくとも、オプション1を実行する別の関数で処理関数をラップすることができます。本当にベクトル化されるわけではありませんが、多くの場合はより洗練されたものになります。 –

+0

解決策1は、forループと比較して非常に高速です。理由を理解したいと思いますが、多分それは別の質問です。 –

関連する問題