0

我正在使用这些数据:

'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);

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

知道有什么问题吗?

4

1 回答 1

1

根据 Gavin 的见解,这里有一个建议:

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

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

于 2012-07-05T15:33:47.310 回答