我在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
更改为要分解的部分,从而使代码得以运行。很奇怪。其他三列绝对是数字。