2017-05-29 40 views
0

多変量時系列データセットのウェーブレット変換を計算したいと思います。 wavethreshパッケージ、特にmodwt()機能を使用する予定です。この関数のヘルプファイルは、オブジェクトが「単変量または多変量時系列」であることを指定します。数値ベクトル、行列およびデータフレームも受け入れます。多変量XTSからTSをRに変換する

現在、私のデータセットはxts zoo形式であり、時間は15分間隔であり、tsに変換したいが、私は大きな困難を抱えている。

私は次のことを試みた

modwtCoeff <- modwt(as.ts(wideRawXTS, 
+       start = head(index(wideRawXTS), 1), 
+       end = tail(index(wideRawXTS), 1), 
+       frequency = 1), 
+      filter = "la8", 
+      n.levels = "10", 
+      boundary = "periodic", 
+      fast = TRUE) 

> class(wideRawXTS) 
[1] "xts" "zoo" 

Iは、上記の行の結果として、次のエラーを表示する場合head(index(wideRawXTS,1),1)戻り"2017-01-20 16:30:00 GMT"tail(index(wideRawXTS,1),1)戻る"2017-02-03 16:00:00 GMT"

Error in ts(coredata(x), frequency = frequency(x), ...) : 
    formal argument "frequency" matched by multiple actual arguments 

私はmodwtラッパー関数を削除し、私はエラーとしてxtsからts変換にある病気は同じエラーを得る。さらにグーグル・グーグルでは、私はこの記事を見に行きましたhttps://www.r-bloggers.com/preventing-argument-use-in-r/しかし、私はそれを完全に取得していません。私の推測では、エラーがas.ts関数でいくつかの引数を使用するのを避けるために、変換を個々のステップに分解する必要があるかもしれないということです。

変換でどこが間違っているのか、誰かが私に少しの方向性を与えることができますか?ここでは再現可能な例を提供するために、wideRawXTSオブジェクトの出力に対してlinkがあります。

答えて

1

周波数を計算するための一般的な機能は次のとおりです。

周波数= number_of_events個/ TIME_INTERVAL

あなたのデータは14日間の時間間隔のための1343行を持っているとして、周波数が何であるかに依存あなたの時間単位。

時間単位:日

この場合は、周波数は次のとおりです。

平均だ
1343/14 = 95.93 => 96 

は、あなたが一日あたり96 measurmentsを作ります。

時間単位:時間

この場合は、周波数は次のとおりです。

平均だ
1343/(14*24) = 3.99 => 4 

は、あなたが毎時4 measurmentsを作ります。

時間単位:この場合、15分

は、周波数は次のとおりです。

平均だ
1343/(14*24*4) = 0.999 => 1 

は、あなたが1 measurment 15分ごとに作ります。

関連する問題