can';在R中绘制预测线
我正在使用这些数据:can';在R中绘制预测线,r,plot,R,Plot,我正在使用这些数据: '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 4
'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,第596页)
model根据Gavin的见解,这里有一个建议:
plot(NDVI, as.numeric(Used)-1 )
lines(xv,yv)
因子是从1L
开始的整数向量,默认情况下按标签的字母顺序赋值。所以你应该同意“否”<“是”导致否为1,是为2,然后向下移动到正确的刻度[0,1]。您可能还需要查看str(yv)
您查看过yv吗?范围(xv)和范围(yv)是否在轴限制内?猜测一下,y轴有一个系数(请参见str()
输出中的$Used
)。这并没有画出你认为的情节。绘制绘图后,请使用par(“usr”)
查看绘图的用户坐标,然后查看range(yv)
,希望这会对您有所启发。您需要一个覆盖范围为0,…,1的连续绘图,当您在y轴上绘制一个因子时,您可能无法获得该范围。您显然也执行了attach(dfrm)
,这通常是不可取的,但如果您要使用attach
,则应将其包含在发布的代码中。
plot(NDVI, as.numeric(Used)-1 )
lines(xv,yv)