我无法拟合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模型适合我的数据吗?