微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

R 未标记动态模型 - 对比错误

如何解决R 未标记动态模型 - 对比错误

第一次发布一个问题,如果做得不好,请不要问我。在运行与季节相关的殖民化和灭绝以及未标记的完全时间相关模型时,我收到“对比”消息中的错误。我用以前的数据集运行了一个非常相似的脚本并且它起作用了,我的变量都没有少于这两个因素,所以我不知所措。 见示例:

dput(sitecovariates)

structure(list(site = c("c1","c2","c3","c4","c5","c6","c7","c8","c9","c10"),habitat = c("forest","forest","intervened","intervened"),z_dis_area_cons = c(0.01417616,1.293259454,-1.239008999,1.450977679,0.996786633,-1.386719679,-0.345627422,0.155103464,-0.369931339,1.402797566)),row.names = c(NA,10L),class = "data.frame")

dput(HistDet)

structure(list(s1v1 = c(0L,0L,1L,1L
),s1v2 = c(1L,NA,0L),s1v3 = c(0L,s1v4 = c(0L,s1v5 = c(0L,s1v6 = c(0L,s2v1 = c(0L,1L),s2v2 = c(0L,s2v3 = c(0L,s2v4 = c(0L,s2v5 = c(1L,s2v6 = c(0L,0L)),class = "data.frame")
season<-as.character(1:2)
season<-matrix(season,nrow(HistDet),2,byrow=TRUE)
Datos <- unmarkedMultFrame(y = HistDet,siteCovs = sitecovariates,yearlySiteCovs=list(s=season),numPrimary=2)
mo0<-colext(~1,~1,Datos)
mo1<-colext(~1,~s,Datos)

#季节依赖的殖民和灭绝

mo2<-colext(~1,~s-1,Datos)

#完全依赖时间

mo3<-colext(~1,Datos)

我很感激任何建议。

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。