さまざまな差異の後の回帰係数の解釈


10

時系列を差分した後に線形回帰係数を解釈する方法(単位根を除去するため)を説明する説明はほとんどありません。単純に、それを正式に述べる必要がないほど単純ですか?

(私はこの質問を知っていますが、それがどの程度一般的であるかはわかりませんでした)。

モデル興味があるとしましょう ここで、はおそらくARMA(p、q)です。興味深いのは、、、...です。具体的には、「 1単位の変化」という観点からの解釈は、、の平均変化をもたらし δ T β 1 β 2 β Pの X I YはT β I iが= 1つの... PをYt=β0+β1X1t+β2X2t++...+βpXpt+δtδtβ1β2βpXiYtβii=1...p.

ここで、ユニットルートからの非定常性が疑われるためにを差分する必要があるとしましょう(例:ADFテスト)。次に、同じように、それぞれのも異なるようにする必要があります。 X i tYtXit

次の場合、の解釈は何ですか?βi

  1. 第1の差分 から取られるとの各? Y t X i tYtYtXit
  2. 第2の差分(差分の差分)()から取られY_ {T}との各X_ {それ}Yt X i tYtXit
  3. 季節差(例えば(1B12) 毎月のデータ用)から取得されYtとの各Xit

編集1

係数の違いと解釈につ​​いて言及しているテキストを1つ見つけました。リンクされた質問とよく似ています。これは、動的回帰を使用したAlan Pankratzの予測119〜120ページからの抜粋です。

ここに画像の説明を入力してください ここに画像の説明を入力してください


時系列は月次であると思いますか?YとXが経済変数の対数変換であることは?

問題は一般的な解釈に関するものであり、おそらくARMAエラーを伴うさまざまな形式の差異が、差異のない回帰から解釈を変更する場合です。したがって、ログには記録されません:)
B_Miner

はい。しかし、解釈はのような単純なようすることができの増加で成長のにおける単位増加のための成長の。ここで、「成長」は、質問1の月ごとの成長であり、質問の「年ごと」の成長です。成長はyの絶対的な成長ですが、yがの対数変換の場合、の相対的な成長です。それはあなたが求めているそのような解釈ですか?、Y 、X 1つの Zβ1yx1z

このコメントは、このトピックに関する私の混乱を助長します。ベータ版は差分後も変更されないため、解釈がまったく変化しない例を見つけましたが、解釈が差分データに変化することを意味する(私は考える)成長という言葉を使用する必要があることを示唆しています(私は) Yの変化、Xの変化)。
B_Miner

1
ここでやや関連する答え。
Richard Hardy

回答:


10

タイプしやすいので、1つの独立変数の例を見てみましょう。

から始めると、についても同じことがます。 、Y T - 1 = β 0 + β 1 X T - 1yt=β0+β1xtyt1=β0+β1xt1

したがって、2つを減算すると、ます。そのための解釈係数ないではない変更、それは同じこれらの式の各インチβ 1 β 1Δy=β1Δx β1 β1

ただし、方程式 の解釈は、方程式の解釈と同じではありません。それが私の言いたいことです。 Δ Y = β 1 Δ Xyt=β0+β1xt Δy=β1Δx

だから、変化であるの単位変化のための、それはまたの成長の変化であるの成長の単位変化のため。、Y 、X 、Y 、Xβ1yxyx

差異の理由は「技術的」です。シリーズが非定常の場合、OLSでを推定できません。差分系列が静止している場合、私はの推定に使用できる式からの推定値としてとして式ことがあるため、同じβ 1、Δ Y = β 1つの Δ X β 1 、Y T = β 0 + β 1 X Tyt=β0+β1xtβ1Δy=β1Δxβ1yt=β0+β1xt β1

そのため、系列が非定常の場合、はで推定値を見つけるための「技術的な」トリックです。このトリックは、同じが微分方程式に現れるという事実を利用しています。、Y T = β 0 + β 1 X T β 1β1yt=β0+β1xt β1

複数の独立変数がある場合、これは明らかに違いはありません。

注:これはすべて、モデルの線形性の結果です場合、なので、は同時に単位のの変化です変更もの成長における単位変化に対するYの成長の変化は、同じである。Δ Y = α Δ X α Y 、X 、X αy=αx+βΔy=αΔxαyxxα


したがって、解釈は双方向です。しかし、重要な点は、差異がある場合(私の質問の3つのタイプのいずれか、またはその組み合わせ)、元の未差異のベータがまだ推定されているということです(そのため、関心のある元の研究の質問はまだ利用可能です)。正しい?Armaのエラーがあった場合でも、それは成り立ちますか?
B_Miner 2015

まあ、微分方程式からを推定する場合、この推定は、れていない方程式の推定でもあります(同じ)。重要なのは、推定を行う方程式では、系列が静止している必要があり、それですべてが問題ないということです(そうでない場合は、不偏性などの望ましい特性を持つ推定量が得られません)。もちろん、この方法ではを推定できないという欠点があります。そのため、推定が必要な場合は、を必要があります。β 1 β 1 β 1 β 0 β 0β1β^1β1β1β0β0

切片が関心の対象となることはめったにありませんが、より重要なのは、関心のある連続変数またはダミー変数の係数であるB1からBPです。そして明確にするために、エラーがiidでない場合、これに関して何も変更されませんが、ARMAエラーを使用しますか?私は、違いがあるかどうかに関係なく、正しい解釈を考慮する必要があると思いますか(「他のすべてが等しい」には、制御されているyの(AR付きの)遅延値が含まれるため)?
B_Miner 2015

1
ARMAエラーによって解釈が変わることはありません。唯一の技術的な問題は、差分の後に定常系列が必要であり、それ以外の場合はの推定値に偏りがあるため、ARMAエラーがある場合に差分の後に定常系列が得られるということですが、私の考えではすべて問題ありません。β1

1
季節差の場合、差分方程式でも「元の」方程式と同じが得られるため、すべてが有効なままです。実際、何を行っても、操作後に同じがあることを示すことができる限り、推論は有効なままです。β 1β1β1

0

最終的な伝達関数を取り、それを純粋な右辺の方程式として再表現します。この形式では、PDLまたはADLになります。通訳は通常通り続きます。私はそのオプションをAUTOBOXに実装し、右手側と呼びました。使用したいデータセットとモデルを投稿していただければ、結果を投稿させていただきます。

編集:等価係数の仮説をテストするための例示的な例を提示する:

ここでhttp://www.autobox.com/stack/GASX.ASCから入手できるBox-JenklinsテキストからGASXデータセット(Xが最初、次にY)を取得し、無差別シリーズの伝達関数を推定して取得しましたここに画像の説明を入力してください

次に、YとXの両方に単純な差分を導入して取得しました ここに画像の説明を入力してください。係数が同じであるという仮説は拒否されます。係数は似ていますが、まったく同じではありません。次に、MA係数(1の近く)を導入して[1-B]を乗算する代数的な演習を完了しようとしましたが、それでも差のない結果は再現されませんでした。

要約すると、答えは異なりますが、それは、差異のない場合にここに画像の説明を入力してください定数項が省略されているためかもしれません。

2つのホワイトノイズ系列(X1とY1)をシミュレートし、定数項なしでそれらのOLSモデルを推定することにし、取得しました。次に、X1とY1の両方の白いnosieシリーズを統合し、2つの新しいシリーズ(X2とY2)を取得しました。以下は、X2 AND Y2 [ ここに画像の説明を入力してください] [4のOLSモデルの結果です。結果の回帰係数はほぼ同じです(X2、Y2研究での観測が1つ少ないため、わずかな変動です。したがって、ケースが証明されたか(またはそうでないか)拒否)回帰係数は同等です。(X1対Y1)に定数を導入した場合、回帰係数は同じではありませんでした。どうやら、基本ケースに定数を組み込んではならないという要件があります(差異なし)。これらは調査結果は@f coppensに同意します。


従わない-伝達関数?どういう意味ですか?
B_Miner

一般的な伝達関数は次の形式をとります:Yt =μ+ [(ω0−ω1B1 −.....−ωsBs)/ 1−δ1B1 −...δrBr)] Xt−b + etここで、etにはいくつかの有理数構造があります
IrishStat

解釈は実際には差分によって変化するというあなたの答えからそれをとりますか?私の質問にあるものから伝達関数を構築する方法がわかりません。βi
B_Miner

差異がない場合のβiの解釈は、Yのレベルが影響を受ける一方で、差異がある場合はYの変化が影響を受けるということです。
IrishStat 2015

私の質問のリンクを見てください。ここでは、差分モデルの解釈はレベルとまったく同じであると言っているようです。これはそうではないことを示唆していますか?答えの違い(しゃれは意図されていません)のように見えるものに混乱しています。
B_Miner
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.