在vglm

时间:2019-04-30 15:07:17

标签: r regression linear-regression glm vgam

我无法拟合Tobit模型,因为显然需要etastart,mustart,coefstart的起始值,尽管我对大约40个协变量有大约8000个观察值。

我的数据包括许多0,我对研究与我的因变量的一些协变量之间的关联感兴趣。

我正在使用vglm软件包中的Tobit规范(不同的软件包出现类似的问题)。

如果我只用几个协变量运行下面的代码,它就可以正常工作。添加更多内容后,我将显示以下错误消息。

m1 <- vglm(y ~ sexe  + age + nb_children + engaged, 
                   tobit(Lower  = 0, Upper = 23),
                   coefstart = eta1,
                   data = full_data_18)
  

vglm.fitter中的错误(x = x,y = y,w = w,offset = offset,Xm2 = Xm2,:    无法获得有效的初始值。尝试使用“ etastart”,“ coefstart”>或“ mustart”,否则使用家族特定的参数,例如“ imethod”。

然后我尝试指定coefstart值(所有协变量+两个截距)a

coef1<- c(0,0,0,0,0,0)

m1 <- vglm(y ~ sexe  + age + nb_children + engaged, 
                   tobit(Lower  = 0, Upper = 23),
                   coefstart = coef1,
                   data = full_data_18)

并收到错误消息:

  

vglm.fitter中的错误(x = x,y = y,w = w,offset = offset,Xm2 = Xm2,:    无法获得有效的初始值。尝试使用“ etastart”,“ coefstart”或“ mustart”,否则使用家族特定的参数,例如“ imethod”。

我不确定如何准确地定义etastart和/或mustart,并且我缺乏理解如何正确定义这些起始值的直觉?我本着“保守先验”的精神在这里选择0。

除了这个起始值问题之外,我还应该在vglm函数中指定其他内容,以确保它可以将Tobit模型适合我的数据吗?

0 个答案:

没有答案