尝试创建depmix模型时出现“模型响应对于二项式模型无效”错误

时间:2019-09-13 04:49:18

标签: r

我在R中使用depmixS4包,我想创建一个将响应建模为二项式函数的depmix模型。我的代码目前看起来像

mod <- depmix(list(y~1),
          data=train, 
          nstate=18,
          family=list(binomial()),
          ntimes=rep(28, 100))

y是数据框中的数字列,仅由0和1组成。当我运行上面的代码时,R引发错误

  

.locl(公式,数据,族,pstart,固定,概率,...)中的错误:模型响应对于二项式模型无效

我在https://github.com/cran/depmixS4/blob/master/R/responseGLM.R处查看了源代码,但无法弄清楚为什么会引发此错误。我还有其他三个二进制数值列,我们称它们为y1,y2,y3,如果我将上述代码中的y替换为这三列中的任何一列,则代码运行良好。

有人可以建议我如何解决该问题或如何调查此问题?

谢谢

编辑:

train <- data.frame(y = c(0, 0, 0, 0, 0, 0), 
                   y1 = c(0, 0, 0, 0, 0, 0), 
                   y2 = c(0, 0, 0, 0, 0, 0), 
                   y3 = c(0, 0, 0, 0, 0, 1))

我最终将列y更改为要分解的部分,从而使代码得以运行。很奇怪。其他三列绝对是数字。

0 个答案:

没有答案