R-код для GEE занимает слишком много времени или не будет работать
Я пытаюсь запустить многомерную порядковую модель GEE с использованием R. Поскольку я работаю с порядковыми данными (по порядку 1, 2, 3), я рассмотрел два пакета, которые нужно использовать: 1) ordgee из geepack и 2) ordLORgee от мульт. Мой результат также коррелирует, поэтому я использую ссылку logit.
Мой набор данных включает более 40 000 продольных наблюдений. Мой код из Multgee и Geepack работает, если я работаю с меньшим набором данных (например, 500, 1000), но когда я пытаюсь запустить эти функции с использованием полного набора данных, требуется более 3 часов, чтобы получить сообщение об ошибке, что мое соединение R потерял.
Из мультя, я пытаюсь соответствовать:
ordLORgee(output ~ day + pre2 + pre3 + pre4, id = person, link = 'logit', data = sampledf)
Из geepack, я пытаюсь соответствовать:
ordgee(ordered(output) ~ day + pre2 + pre3 + pre4, id = person, data = final)
Может кто-нибудь сказать мне, откуда эта проблема? В идеале я хочу использовать пакет multgee, поскольку этот метод package/GEE является относительно новым, и я не мог получить большую помощь в Интернете. Спасибо.