線形回帰+交絡因子
交絡因子Zによって調整された結果Yと変数Xの間の効果サイズと有意性にアクセスしたいとします。 私の質問は、以下のシナリオ間でXの効果サイズと有意性を決定するための違いがある場合です。 変数と交絡因子を線形回帰モデルにまとめます。これは、Y〜X + Zの回帰モデルをちょうど当てはめて、係数とXの p値を計算することを意味します。 残留、取得RからY〜Zを、そしてその後の適合回帰モデルR〜Xは、その後、係数のそのp値を算出X(R〜Xからの)。 ここから交絡者を学びます。 編集----- @Gordon Smythの回答に感謝します。ただし、シミュレーション調査(以下のコード)から、メソッド1、メソッド2、およびゴードンスミスの回答のメソッド3の誤検出率を比較したところ、驚くことに、メソッド2の誤検出率がかなり低いことがわかりました。 方法1は「教科書」として正しいことを理解しています。論理的にmethod2の何が問題なのでしょうか?その上、「すべてのモデルが間違っていますが、一部は便利です」。 p1 = p2 = p3 = c() i=0 while(i<10000){ y = rnorm(10) x = rnorm(10) c = rnorm(10) # method 1 p1[i] = summary(lm(y~x + c))$coefficients[2,4] # method 2 p2[i] = summary(lm(lm(y ~ c)$res ~ x))$coefficients[2,4] # method …