Я работаю над проблемой моделирования с использованием функции gee в R. Пример:
gee(Y~X, family="binomial", data=mdata, id=id, constr="exchangable")
Это парная двоичная проблема. mdata моделируется n=100 раз. Пример набора данных:
Y X id
1 1 1
0 0 1
1 1 2
1 0 2
1 1 3
0 0 3
1 0 4
0 0 4
Истинные коэффициенты для пересечения и X равны 0,4 и 1. Для большей части моделирования функция gee дает довольно хорошие оценки, например,
(intercept) X
0.44143151 1.05925
Однако есть одно моделирование, которое приводит к странной оценке
(intercept) X
-19.320225 20.441242
И в то же время мой R застревает в этом моделировании.
Мне интересно если это из-за неконвергентности функции. Если так, могу ли я пропустить это, когда это произойдет?
Надеюсь на некоторую помощь!
Спасибо,
Аллен