为什么我的二项式GLM有多个预测行？

Question

我的数据是以千克为单位的后代数量，以及一列1和0来表示母亲是否已步入晚年。

Chick Mass        Terminal Effect
   3.4                   0
   3.1                   1
   2.4                   1
   3.6                   0

等等。

因此，我有一个模型来评估质量(以公斤为单位)是否对死亡率(二项式)有影响。

m10 <-glm(Terminal_Effect~chick_mass, data = cranesData, family = binomial(link="logit"))
summary(m10)

plot(cranesData$Terminal_Effect~cranesData$chick_mass, xlab = "Chick Mass (kg)", ylab = "Probability of Mother Death", pch = 19)

当我绘制这个图时，我的图上有多条线，有没有办法把它改成一条线？

​

​

如有任何帮助，我们将不胜感激:)

Answer 1

在绘图之前，对您的预测值进行排序。

使用带有一些变体的虹膜数据集：

set.seed(111)
dat = iris
dat$Species = as.numeric(dat$Species=="setosa")
dat$Petal.Width = dat$Petal.Width + rnorm(nrow(dat))

对预测器进行排序，在本例中为petal.width：

dat = dat[order(dat$Petal.Width),]

fit = glm(Species ~ Petal.Width,data=dat,family="binomial")

plot(dat$Species ~ dat$Petal.Width)
lines(dat$Petal.Width,fit$fitted.values,col="blue")

​

​