「交互作用項を入れたときの主効果には意味がない!」などという人がいますが、意味がないことはありません。意味がないと言いたくなる人の気持ちも、場合によってはわかります。年齢と喫煙の交互作用項をモデルに入れる場合、年齢の主効果は喫煙=0の時の効果、喫煙の主効果は年齢=0の時の効果になります。0歳の時の喫煙の効果を調べることに意味があるかと言えば、0歳で喫煙している人はいないでしょうから、「意味はありません」。それが何にも意味がないかと言われれば、今言ったような内容が結果の解釈となるわけです。1逆に、解釈する意味のある値にしたいのであれば、平均年齢が0になるように変数を変換するなどの方法が考えられます。Stataで交互作用項を入れる場合には「♯」を使います。i.A#i.Bのいずれかが一般的でしょう。「i.」はその変数が離散変数20/1とか。variableという。あることを宣言する接頭辞です。上記二つの方法に加えて、AとBの交互作用項を事前に作るという方法もあります。Int交互作用を検討する変数が連続変数である場合には、「c.」を変数名の前につけましょう。c.C##c.D上記例だと一つ目が離散変数Aと連続変数Cの交互作用、二つ目は連続変数Cと連続変数Dの交互作用になります。ちなみに、離散変数どうしの交互作用項を見るときに、♯ひとつだけをつかうと、交互作用を検討している変数が取りうる値の組み合わせごとの推定値を返します。B=1、それぞれの場合の値を算出します。groupと比較した場合の数値です。それらの値単独では交互作用があったかどうかは分かりませんが、それぞれの組み合わせの場合にどのような値を取ることになるのか確認するという意味では便利です。結果の説明をしやすいように年齢階級を減らしました。II・III度高血圧を収縮期血圧のみで定義しました。3かなり雑な定義です。♯ひとつをつかって交互作用を定義した場合です。ロジスティック回帰分析にしています。estimatorなどが望ましいです。♯ふたつをつかって交互作用を定義した場合です。marginを使用しています。sexが0と1、agegrpが1と3の場合に取りうる組み合わせごとに予測値を計算しろという指令です。6行目の後にやっても、5行目の後にやっても同じ結果を返します。marginで計算した値をグラフにするコマンドです。mmHg以上になるオッズ比を計算しています。.9330996------------------------------------------------------------------------------8行目の結果:sex#agegrpが有意ではないので、交互作用があるとは言えないという結論になります。.9330996------------------------------------------------------------------------------これを図にするとこうなります。interactionがあります。interaction(具体例2)になります。interactionを計算することもできます。interactionを計算するのかという理由について深く知りたい人は、以下のリンクを参考にするといいと思います。interactionをどのように計算するか示しておくだけにします。交互作用項を別途新規の変数として定義する必要があります。scaleでも交互作用があるとは言えないということです。コメントをキャンセルメールアドレスが公開されることはありません。センターで研究員をしている井上陽介の個人的な情報・記録の置き場所です。中国海南島の調査からアメリカ・ノースカロライナでのポスドクを経て、東京にたどり着きました。ご質問がある方はメールにてご連絡ください!yoinoue[at]hosp.ncgm.go.jp※[at]を@に変えてください。少数民族居住地域で、急速な経済発展にともなう生活環境の変化が、健康にどのような影響を与えるか評価するために、フィールドワークをしています。仕事は村人と酒を飲むことです。2016年7月からアメリカ・ノースカロライナに滞在中。