循環共変量(GAMなど)を処理する方法は何ですか?


8

私は、いくつかの共変量が「円」上にあるモデルを構築しています。これは、それらが[0,1)の間隔で値をとり、0 = 1であるという意味です。私はこの状況に対処するためのテクニックについて疑問に思っています。1つのアイデアは、循環変数thetaを変数のペア(sin(theta)、cos(theta))として表すことです。このアプローチまたはより良いアプローチについて何か考えはありますか?

特にmgcvパッケージのGAMを使用しています。特定の加法破片が端点で同じ値を持つ必要があることをモデルに伝える方法はありますか?別のパッケージ?

ありがとう!


2
あなたが言ったように私はサインとコサインしか見たことがありません。同様の質問をするサイト上の別の質問、方向データをIVとするロジスティック回帰があります。で、この質問 SOイテレータのオーバーを指しコメントがあるcircular可能性が関心のあるパッケージを。うまくいけば、あなたはより良い反応を得るでしょう。
Andy W

2
(平方可積分)周期関数の場合、@ Andyの正弦と余弦は唯一の基底ではなく、唯一の直交基底でさえありません。彼らのメリットの1つは、しばしば物理理論が他の変数との関係にそれらを巻き込むことです。これは、DVとこれらの共変量の間の仮説または予想される関係に鋭く関心を持つべきであることを示唆しています。DavidR、これについて何を教えてくれますか?
whuber

周期的な時間的プロセスをモデル化しようとしています。私の元のデータはイベントの日付と時刻であり、これらのイベントの割合を経時的にモデル化したいと思います。手始めに、時間、曜日、年の時間に関して周期性があると思います。これらは明らかな循環変数です。ポアソンGAMから始めます。これらのそれぞれの影響を個別に調べ、将来のレート予測を行うことに興味があります。
DavidR 2011年

回答:


6

循環変数を処理する方法は2つあります。1つのハックな方法は、境界条件のいずれかの側でデータセットを手動で複製することですが、より洗練されたソリューションは、周期的な境界条件で組み込みのスプライン基底関数を使用することだと思います。 !

例えば:

bs="cc"循環3次回帰スプラインを指定します(を参照cyclic.cubic.spline)。つまり、2次微分までの、両端が一致するペナルティ付き3次回帰スプライン。

球上のスプライン

bs="sos"。これらは球上の2次元のスプラインです。引数は緯度と経度であり、球の薄板スプラインに類似しています。等方性が適切な場合に、地球の大部分でサンプリングされたデータに役立ちます。詳細Spherical.Splineについては、を参照してください。

bs="cp" P-スプラインの循環バージョンを与える


完璧ですね!これやってみます 私はハッキーな方法について考えていましたが、複数の循環変数があるため、問題のさまざまな対称性をすべて表示するには、膨大な数の重複データポイントを作成する必要があると考えていました。
DavidR

@DavidR低周波の現象に円形の3次スプラインを使用して、うまくいきました。モデル化したい最高周波数を表すのに十分である必要があります。これは、毎時と毎日の期間を除外しますが、季節性に適している可能性があります。前者の2つについては、節約は、予想される各周波数の正弦や余弦などの小さな基準から始めることを推奨します。
whuber

@whuber、あなたはDavidRのモデリング手法を誤解していると思います。Davidは、時間、曜日などに個別の変数を使用しています。そのため、それぞれに個別のスプラインを適用します。(似たようなことをしているので、あなたの考えを歓迎します)
seanv507 2017年

Whuber、私はあなたの意見に同意します。たとえば、12月のより重いショッピングをモデル化するには、月ごとの変数に高次のスプラインが必要になるということです。
seanv507 2017年

0

Gill and Hangartner(2010)を調べてみてください政治学における循環データとその取り扱い方法。彼らは、循環/時計/季節データのさまざまなモデルについて話し、Jeff Gill がインスピレーションを得るために調査できる論文のRコードを提供しています。方法論とRコードを一緒に織り込むこの資料のプレゼンテーションバージョンがあるはずです。


循環データに関する文献へのこのポインタをありがとう。それは良い出発点のようです。私はすぐにこの特定の記事を読み飛ばしましたが、循環共変量ではなく循環応答を扱っているようです。
DavidR
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.