观察值与拟合值的关系图

Question

我想绘制观察到的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缺失。我做错了什么？

Answer 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")