0

次のように、R の glm 関数を使用してロジスティック回帰モデル (「mylogit」と呼ばれる) を構築しました。

mylogit <- glm(answer ~ as.factor(gender) + age, data = mydata, family = "binomial")

ここで、年齢は数値、性別はカテゴリ (男性と女性) です。

次に、構築したモデルを使用して予測を行いました。

pred <- predict(mylogit, type = "response")

次のようにして、予測の時系列プロットを簡単に作成できます。

plot.ts(ts(pred))

次のようなプロットを提供します。

予測に対する時間のプロット

これにより、予測のプロットが得られます。

私の質問はこれです:glmで指定された性別(男性または女性)に従ってx軸をセグメントに入れることは可能ですか? つまり、y 軸に予測、x 軸に性別 (男性と女性に分けられます) を設定できますか?

プロットしたいデータのサンプルは次のようになります。

やった：

bind = cbind(mydata, pred)

「バインド」は次のようになります。

pred          age        gender
0.9461198     32          male
0.9463577     45         female
0.9461198     45         female
0.9461198     37         female
0.9477645     40          male
0.8304513     32         female
4

2 に答える 2

0

あなたが持っているデータは時系列ではないので、使用する必要はないtsと思いますよね？プロットする前にplot.ts並べ替えるだけです。pred

# Get data
str <- "pred,age,gender
0.9461198,32,male
0.9463577,45,female
0.9461198,45,female
0.9461198,37,female
0.9477645,40,male
0.8304513,32,female"
bind <- read.csv(textConnection(str))

# Plot
bind <- bind[order(bind$gender),]
plot(bind$pred, col = bind$gender)

library(ggplot2)
ggplot(bind, aes(x = gender, y = pred)) +
  geom_point(position = position_jitter(width = .3))

または、作成せずにbind行うことができますplot(pred[order(mydata$gender)])。

于 2016-02-25T09:05:42.803 に答える