元データ分析の会社で働いていた人の四方山話スマートフォン用の表示で見る元データ分析の会社で働いていた人の四方山話この広告は、90日以上更新していないブログに表示しています。2009-07-25Rでウェーブレット変換をしたい!R結論から言うとまだできていませんよ。期待して読まない方がいいですよ。諸般の事情で、Rでウェーブレット変換をすることになりました。「信号処理」が研究グループの名前に冠されており、研究室のボスは信号処理の大御所らしいのですが、早々と信号処理からはドロップアウトして現在信号処理とは関係ないことをしているというゆとりっぷり全開な状況ですので、現在のスキルは以下のようになります。ウェーブレット変換は名前を知っている。最近Rをちょいちょい勉強しているのでRレベル3くらい(自称)このスキルである問題点は、ドキュメントを読んでもウェーブレット用語が理解できない可能性が高いということですね。この関数は○○を××します、と言われても処理の対象である○○の意味が分からず、従って処理内容である××が何を行っているのか分からない、という2重苦に陥りそうで怖い。それはともかく、ざっと調べたところあまり良さそうな入門記事が見当たらない。唯一日本語でまともにRでウェーブレット変換について触れていたのはここ。ウェーブレット解析4Rwaveが主に紹介されていたみたいので、色々調べてみたところ、こいつは連続ウェーブレット変換用のパッケージらしい。今欲しいのは離散ウェーブレット変換用のパッケージなんだぜ、ということで、waveletsとかいうパッケージを見つけた。transforms:MODWT)及び、それらの逆変換を計算し、プロットするための関数郡を含んでいる。多重解像度分析と同様に、上述の分解において使用されるウェーブレット変換のフィルターを計算しプロットするための機能を含んでいる。これだ!というわけで関数とかを調べていく。方針DescriptとUsageをざっと見ていって、それっぽいのを探す。values)を用いてクラスdwt,modwt,dwptそして,modwptのウェーブレット(ハイパス)とスケーリング(ローパス)係数の両方を揃える。valuesが何かよく分からない上に何か二つ係数があるのかー。2つの係数を揃えるらしい。それがどういう意味があるのかはよく分からない。method)を使用するかどうかの理論値inverseもとの場所へ揃えられたオブジェクトのウェーブレットとスケーリング係数をシフトするかどうかの理論値dwt:離散ウェーブレット変換オブジェクトDescription離散ウェーブレット変換の情報を含んでいるS4オブジェクト未だにS4オブジェクトがよく分かっていない><ああ、これはオブジェクトの話だな。こっちが関数の方ですね。dwt:離散ウェーブレット変換Description一変量あるいは多変量の時系列データの離散ウェーブレット変換係数を計算する。fact=TRUE)ArgumentsX一変量あるいは多変量時系列データ。数値のベクトル、マトリクス、データフレームが受け付けられる。filterwt.filterオブジェクト、分解においてどのウェーブレットフィルターを用いるかという文字列、あるいは、ウェーブレット係数(スケーリング係数ではない)数値のベクトル。受け付けられるフィルター名のためにhelp(wt.filter)を参照のこと。n.levels分解レベルを特定する整数値。デフォルトでは、これはXの長さが少なくとも、レベルJのウェーブレットフィルターの長さが同じ大きさとなる最小のJであるが、レベルJ+1のウェーブレットフィルターの長さよりは小さい。