相互作用効果を得るために係数を追加する-SEで何をすべきか?


13

相互作用を含む多変量回帰があります。たとえば、最貧五分位の治療効果の推定値を取得するには、治療回帰変数の係数を相互作用変数(治療と五分位1と相互作用する)の係数に追加する必要があります。回帰から2つの係数を追加するとき、どのように標準誤差を取得しますか?2つの係数から標準誤差を追加することは可能ですか?t-statsはどうですか?これらも追加することは可能ですか?私は推測していませんが、これに関するガイダンスを見つけることができません。

よろしくお願いします!


これは本当に役に立ちます!Rで似たようなことをしようとしていますが、グループ間でサンプルサイズが少し異なります。同じ式を使用して2つのエラーを組み合わせて、新しい標準を取得できますか。エラー?助けてくれてありがとうクリスタル
クリスタル

1
ちょっと@クリスタル-サイトへようこそ!これは良い質問ですが、新しい質問として提示する必要があります(右上の[質問]ボタン)。今、あなたはこの古い質問への「答え」としてそれを提出しました。この質問のURLをコピーして新しい質問に貼り付けるだけで、あなたが話していることをすべて理解できます。
マットパーカー

回答:


10

これは表現だと思います。SEbnew

SE12+SE22+2Cov(b1,b2)

この新しい標準エラーを使用して、H oをテストするための新しいテスト統計を見つけることができますβ = 0Ho:β=0


こんにちはサラ、答えられたと思うならこの質問を閉じてください。
suncoolsu

こんにちは-回答ありがとうございます。Stataを使用していることを忘れていました。(Stataからの出力を使用して)2つの係数を一緒に追加するとき、標準エラーも追加できますか もしそうなら、係数の合計を標準誤差の合計で割ることによって標準誤差を取得できるはずです。同意しますか?再度、感謝します。
サラ

サラ、Stataでは、「lincom」機能を使用します。変数var1とvar2があり、var1の係数の3倍とvar2の係数の2倍を追加するとします。「lincom 3 * var1 + 2 * var2」と入力します。これにより、この推定の標準誤差と信頼区間が得られます。
チャーリー

5

「多変量」ではなく「多変数」回帰を意味すると仮定します。「多変量」とは、複数の従属変数を持つことを指します。

連続的な予測変数を取得し、それを間隔に分割することは、受け入れられる統計的手法とは見なされません。これにより、残留交絡が発生し、いくつかの相互作用が主効果の一部の適合性の欠如(ここでは不適合)を反映するため、相互作用が誤解を招くほど重要になります。外側の五分位には多くの原因不明のバリエーションがあります。さらに、「五分位効果」を正確に解釈することは実際には不可能です。

関心のある比較のために、それらを予測値の違いとして想定するのが最も簡単です。R rmsパッケージを使用した例を次に示します。

require(rms)
f <- ols(y ~ x1 + rcs(x2,3)*treat)  # or lrm, cph, psm, Rq, Gls, Glm, ...
# This model allows nonlinearity in x2 and interaction between x2 and treat.
# x2 is modeled as two separate restricted cubic spline functions with 3
# knots or join points in common (one function for the reference treatment
# and one function for the difference in curves between the 2 treatments)
contrast(f, list(treat='B', x2=c(.2, .4)),
            list(treat='A', x2=c(.2, .4)))
# Provides a comparison of treatments at 2 values of x2
anova(f) # provides 2 d.f. interaction test and test of whether treatment
# is effective at ANY value of x2 (combined treat main effect + treat x x2
# interaction - this has 3 d.f. here)

2

RRβRV^RV^Rχr2r 行列内の行数です(行が線形に独立していると仮定)。


ありがとう。私は統計の専門家ではなく、私の質問が明確であったかどうかわからないため、別の質問をします。
サラ
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.