У меня проблема с запуском glm Модель, мои данные имеют только 9 субъектов с 3 периодами записи измерений, до, немедленно, после. это часть кода, который я пробовал.
m1<-glmer(sleeve~pain+stiff+diff+(1|patient)+(1|time),data = oa3,
family = binomial(link="logit"))
Обычно, когда я запускаю свой код, я получаю сообщение об ошибке или предупреждение: модель не сходится, проблема сингулярности, и если я помещаю все переменные, я получаю предупреждение на
невозможно оценить масштабированный градиент. Гессиан численно сингулярен: параметры не определены однозначно
проблема с моими данными в том, что мой уровень группы случайных эффектов <5 и мой Размер выборки мал со многими переменными. Я читал, что люди предлагают использовать <a href="/questions/tagged/bayesian" class="post-tag" title="show questions tagged 'bayesian'" rel="tag"> байесовский (предварительная информация) для решения этой проблемы, но я не знаю об этом и не могу понять, как выбрать предыдущую информацию. это я имею в виду blme, который добавляет байесовский пакет в lme4.
Это мои данные https://drive.google.com/file/d/19iTyWM4u11N9bNlTcU2QHZOMzkZ9IEdM/view?usp=sharing
Итак, мой вопрос: можно использовать надежный и glmm для решить мою проблему? есть ли пакет для надежного glmm?