我正在根据该模型的估计值绘制一个线性比例图,该密度图是真实比例下密度与上半正态模型(hn.l.s)中协变量的关系。如何摆脱2条错误消息并绘制线图?
`SELECT list is not in GROUP BY clause and contains nonaggregated column '.id' which is not functionally dependent on columns in GROUP BY clause; this is incompatible with sql_mode=only_full_group_by`
.local(obj,...)中的错误:
不能直接向后转换长度大于1的无标记估计。
dist <- seq(0, 60, 15)
obs <- as.matrix(dat[, 1:4])
site_covs <- data.frame(shrub = dat[, c("d.shrub")], light=dat[, c("light")])
umf<-unmarkedFrameDS(y = obs, siteCovs = site_covs, dist.breaks = c(dist),
survey = "point", unitsIn = "m")
hn.l.s <- distsamp(~light~shrub, umf, keyfun = "halfnorm", output = "density",
unitsOut = "ha")
pred.hn.l.s <- predict(hn.l.s, "state", se.fit = TRUE)
backTransform(hn.l.s, type = "state")
xy.coords(x,y,xlabel,ylabel,log)中的错误:
“ x”和“ y”的长度不同