0

ロジット モデルを実行し、多くの重要な予測変数の確率曲線をプロットしようとしています。Effectsパッケージでこれをやろうとしています。

df=data.frame(income=c(5,5,3,3,6,5),
              won=c(0,0,1,1,1,0),
              age=c(18,18,23,50,19,39),
              home=c(0,0,1,0,0,1))
str(df)

md1 = glm(factor(won) ~ income + age + home, 
          data=df, family=binomial(link="logit"))

summary(md1)

plot(effect("income", md1), grid=TRUE)

しかし、収入 (またはその他の予測因子) に基づいて勝った確率 (応答) を示すようにグラフをプロットする方法を知りたいです。

ただし、私がやりたいのは、同じプロットを生成することです。y 軸にウォン、x 軸に収入を示しますが、曲線は年齢と家の確率を (別々のプロットで) 示しています。

効果のドキュメントでこれを行う方法がわかりません。ヘルプ！

ありがとう。

4

1 に答える 1

1

私の理解が正しければ、相互作用項をモデルに導入する必要があります。例えば、

md1 <- glm(factor(won) ~ income*age + home,
  data=df, family=binomial(link="logit"))

with(df, mean(age))
#[1] 27.83333
with(df, mean(age)-sd(age))
#[1] 14.33642
with(df, mean(age)+sd(age))
#[1] 41.33025

plot(effect("income*age", md1, 
  xlevels=list(income=0:10, age=c(14.3, 27.8, 41.3))), multiline=TRUE)

年齢のさまざまな値で獲得した収入の影響を示します。

于 2012-08-13T21:26:18.067 に答える