R语言unmarked包报错：obsCovs数据框行数不正确求助

核心问题定位 #

unmarkedFrameOccu对obsCovs的行数要求是站点总数 × 每个站点的重复调查次数，而非与站点数或检测矩阵行数一致。新手常误将obsCovs设为和站点数相同的行数，这是触发该错误的最常见原因。

分步排查步骤 #

• 确认检测矩阵结构：先检查你的检测数据y是否为矩阵格式，执行dim(y)查看结果，格式应为[站点数, 调查次数]。
• 计算obsCovs所需行数：用公式nrow(y) * ncol(y)算出obsCovs必须的行数，对比你当前obsCovs的nrow(obsCovs)结果，若不相等则直接定位问题。
• 检查obsCovs数据格式：确保obsCovs是长格式——每个站点的每一次调查对应一行数据，而非每个站点一行包含所有调查的协变量（比如wind1、wind2这类宽格式列需要转换）。

正确格式示例 #

假设2个站点，每个站点调查2次：

# 检测数据（2站点×2调查）
y <- matrix(c(1,0,1,1), nrow=2, ncol=2)
# 站点协变量（2行，对应2个站点）
siteCovs <- data.frame(habitat=c("forest", "grass"))
# 观测协变量（4行=2×2，对应每一次站点调查）
obsCovs <- data.frame(wind=c(1,3,2,1))
# 创建unmarkedFrameOccu对象
umf <- unmarkedFrameOccu(y=y, siteCovs=siteCovs, obsCovs=obsCovs)

宽转长格式处理 #

如果你的原始观测协变量是宽格式（每个站点一行，含多列调查数据），用tidyr包转成长格式：

library(tidyr)
# 假设原始数据中wind1、wind2是两次调查的风速
obsCovs_long <- pivot_longer(你的原始数据框, 
                             cols = starts_with("wind"), # 匹配所有调查协变量列
                             names_to = "survey", 
                             values_to = "wind")
# 保留需要的观测协变量列，确保行数符合要求
obsCovs <- obsCovs_long[, c("wind")]

内容的提问来源于stack exchange，提问作者Leonardo Villa