观测值与拟合值的图答案

【问题标题】：Plot of observed values vs fitted values观测值与拟合值的图
【发布时间】：2014-11-17 12:51:13
【问题描述】：

我想绘制观察到的 X 与观察和拟合的 Y。这是我的代码：

conc <- c(1.6907, 1.7242, 1.7552, 1.7842, 1.8113, 1.8369, 1.8610,1.8839)
number <- c(59,60,62,56,63,59,62,60)
dead <- c(6,13,18,28,52,53,61,60)
y <- dead/number

minimal.model <- glm(y ~1, family = binomial(link=logit), weights=number)
plot(conc, minimal.model$fitted)
points(conc,y,pch=4)

这是我需要得到的

这就是我得到的

请注意，一些观察到的 Y 缺失。我做错了什么？

【问题讨论】：

标签： r plot logistic-regression

【解决方案1】：

我不知道 x 轴标签是怎么回事：您可能想尝试打开一个新的图形窗口。

对于点，需要适当设置y轴范围

par(las=1,bty="l") ## cosmetic
plot(conc, minimal.model$fitted,ylim=c(0,1))
points(conc,y,pch=4)

您也可以使用matplot() 同时绘制两者：

matplot(conc,cbind(minimal.model$fitted,y),
        col=1,pch=c(1,4),ylab="y")

【讨论】：

谢谢本。我喜欢 matplot。