Question

我正在研究具有二项式链接的广义线性混合模型我的因变量是107个受试者中的每一个的从0到12的分数（63个受试者是对照，44个受试者具有疾病A）。导致得分的测试用不同的版本重复2次我想测试是否存在组间差异（控制VS疾病A），还有版本差异和组与版本之间的相互作用。

我使用glmer中的lme4：

glmer(cbind(score, 12-score) ~ gender + age + group + version + group:version  + (1|id), data = data, family="binomial")



Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: cbind(score, 12-score) ~ gender + age + group + version + group:version + (1 | id)
   Data: data
     AIC      BIC   logLik deviance df.resid 
   764.7    788.2   -375.4    750.7      206 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-6.1421 -0.6240  0.3693  0.7269  3.4653 

Random effects:
 Groups Name        Variance Std.Dev.
 id     (Intercept) 0.3852   0.6207  
Number of obs: 213, groups:  id, 107

Fixed effects:
                              Estimate Std. Error z value Pr(>|z|)    
(Intercept)                  3.4707862  0.0017104  2029.2   <2e-16 ***
genderH                      0.3402744  0.0017093   199.1   <2e-16 ***
age                         -0.0152378  0.0009988   -15.3   <2e-16 ***
groupCTRL                    0.9554189  0.0017101   558.7   <2e-16 ***
versionunknown              -2.0853952  0.0017089 -1220.3   <2e-16 ***
groupCTRL:versionunknown     0.1156636  0.0017092    67.7   <2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) SexeH  age    gropNC fmlrtn
SexeH        0.000                            
age         -0.016 -0.006                     
groupeNC     0.001  0.000 -0.008              
familrtncnn  0.000  0.000 -0.013  0.000       
grpNC:fmlrt  0.000  0.000 -0.007  0.000  0.000
convergence code: 0
Model failed to converge with max|grad| = 0.0420639 (tol = 0.001, component 1)
Model is nearly unidentifiable: very large eigenvalue
 - Rescale variables?

我按照此链接中的不同步骤操作：https://rstudio-pubs-static.s3.amazonaws.com/33653_57fc7b8e5d484c909b615d8633c01d51.html以处理收敛失败这些步骤都不能解决这个问题。

此外，除年龄系数外，估计固定系数的标准误差几乎相同。所以我可以看到这里有一个真正的麻烦。

有人可以帮助我吗？

来自lme4包的glmer中的收敛警告+估计固定系数的奇怪标准误差

0 个答案: