Я не могу подогнать модель Tobit, так как, по-видимому, требуются начальные значения для etastart, mustart, coefstart, хотя у меня есть около 8000 наблюдений для примерно 40 ковариат, которые меня интересуют.
Мои данные включают в себя множество нулейи я заинтересован в изучении связи некоторых ковариат с моей зависимой переменной.
Я использую спецификацию Tobit из пакета vglm (аналогичная проблема возникла с различными пакетами).
Если я запускаю приведенный ниже код только с несколькими ковариатами, он работает просто отлично.Как только я добавляю больше, я получаю сообщение об ошибке, показанное ниже.
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 и / илидолжен точно, и мне не хватает интуиции, чтобы понять, как правильно определить эти начальные значения? Я выбрал 0 здесь в духе «консервативного априора».
Помимо этой проблемы начальных значений, есть ли что-то еще, что я долженуказать в функции vglm, чтобы она соответствовала модели Tobit моим данным?