woma c .pql <-glmmPQL (тип ~ боль + жесткая + разница, случайный = ~ 1 | id, данные = mydata, семейство = бином, niter = 10) </p>
Value Std.Error DF t-value p-value
(Intercept) 1.3009843 4.079884 15 0.3 0.7542
pain 0.0415517 0.000000 15 885745.6 0.0000
stiff -0.0183675 0.000000 15 -249125.8 0.0000
diff -0.0077214 0.000000 15 -519945.0 0.0000
Мне трудно принять этот результат, поэтому я попытался увеличить число niter.however, когда я увеличиваю значение для niter = 30, я получаю эту ошибку, которая, как мне кажется, из-за сбоя моей модели. как я уже пробовал использовать другую функцию для glmm, такую как glmer в lme4, bglmer в blme и brm в brms. из всего пакета, с которым я столкнулся с проблемой сближения.
Error in solve.default(estimates[dimE[1] - (p:1), dimE[2] - (p:1), drop = FALSE]) : system is computationally singular: reciprocal condition number = 2.06865e-18
итак, результат, который я получил с niter = 10, был достаточно хорошим или неправильным, чтобы принять его?