我刚刚得知 lmeNB 包的存在，并且考虑到我有过度分散的正计数数据数据，它似乎就在我的袖子里。但是，我无法让它工作。
包的语法可以在这里找到
我正在尝试运行一个包含 37 个人和每个个人 16 年的完美平衡的面板（我的数据框按字母顺序排列，并且随着年份的增加，我认为这与语法中的描述一致。
library(lmeNB)
mle.ar1.fun(abs_pb_t ~ RI1 + RA1 + abs_pb + cddom + cddom2, data=s.data, ID = p.data$ID, Vcode = p.data$Year
现在我收到以下错误消息：
Error in[<-.factor (tmp, upID == uniID[i], value = c("1--1", "1--2", :
NAs are not allowed in subscripted assignments
In addition: Warning messages:
1: In [<-.factor (tmp, upID == uniID[i], value = c("1--1", "1--2", :
invalid factor level, NA generated
2: In [<-.factor (tmp, upID == uniID[i], value = c("1--1", "1--2", :
invalid factor level, NA generated
我认为警告中确定的问题会导致错误，但我的 data.frame 中没有单个缺失值，也没有完美的分离。
有谁之前经历过这个吗？
西蒙