r - 无法在 R 中绘制预测线

Question

我正在使用这些数据：

'data.frame':   1584 obs. of  3 variables:
$ Individual: Factor w/ 3 levels "AG201","AG202",..: 1 1 1 1 1 1 1 1 1 1 ...
$ Used      : Factor w/ 2 levels "no","yes": 2 2 2 2 2 2 2 2 2 2 ...
$ NDVI      : int  4724 4576 4894 4297 4670 4932 4346 3810 3481 4058 ...

我正在使用“NDVI”作为连续解释变量进行 glm，然后通过数据的散点图绘制模型（我正在复制与 R 书中相同的脚本，Crawley，p.596）

model<-glm(Used~NDVI,binomial);
xv<-seq(0,10000,0.2);
yv<-predict(model,list(NDVI=xv),type="response");
plot(NDVI,Used);
lines(xv,yv);

我的问题是我的图表上没有出现任何线条......

知道有什么问题吗？

Answer 1

根据 Gavin 的见解，这里有一个建议：

plot(NDVI, as.numeric(Used)-1 )
lines(xv,yv)

因子是整数向量1L，默认情况下以标签的 alpha 顺序从分配开始。因此，您应该可以接受“否”<“是”，导致“否”为 1，“是”为 2，然后向下移动到正确的比例 [0,1]。您可能还需要查看str(yv)