为什么我的二项式 GLM 会得到多条预测线?
Why do I get multiple predictions lines for my binomial GLM?
我的数据是以千克为单位的后代质量,以及一列 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)
当我画这个的时候,我的图上有多条线,有没有办法把它改成一条线?
任何帮助将不胜感激:)
在绘图之前对您的预测变量值进行排序。
使用带有一些变化的鸢尾花数据集:
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")
我的数据是以千克为单位的后代质量,以及一列 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)
当我画这个的时候,我的图上有多条线,有没有办法把它改成一条线?
任何帮助将不胜感激:)
在绘图之前对您的预测变量值进行排序。
使用带有一些变化的鸢尾花数据集:
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")