我非常感谢有关我matlines()的错误的任何建议。这是一些模拟数据
al<-seq(1,65)
d<-seq(1,65)
r<-c("a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","c","c","c","c","c","c","c","c","c","c","c","c","c","c","c","c","c","c","c","c","c","c","c","c","c",)
v<-c("a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","a","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","b","c","c","c","c","c")
lq<-seq(1,65)
data11<-data.frame(al,d,r,v,lq)
这是模型和提取的预测值。
m11<-glm(formula = al ~ r+log(d) + v +offset(lq), family = poisson, data = data11);summary(m11)
p22<-predict(m11,newdata=data11,se.fit=TRUE,type="response")
data11<- data.frame(data11,p22)
然后尝试绘制它。
plot(al~d, data=data11, pch=as.character(v))
y <- unstack(mydata, mydata$fit~mydata$vegetation)
x <- unstack(data11, data11$d~data11$v)
matlines(x,y)
但是我收到了这个错误：
xy.coords中的错误（x，y，xlabel，ylabel，log = log）： （list）对象不能被强制输入'double'