[R] what does this warnings mean? and what should I do?

ronggui 042045003 at fudan.edu.cn
Tue Dec 13 21:20:34 CET 2005


I use lmer to fit a mixed effect model.It give some warnings.what does this warnings mean? and what should I do?

> (fm2.mlm <- lmer(qd ~ edu + jiankang + peixun +hunyin + cadcj + age + age2 + sex + dangyuan + Comp.1 + Comp.2+trust.cz1 +(trust.cz1|commid), data = individual,na.action = "na.exclude",family="quasibinomial")) 
Generalized linear mixed model fit using PQL 
Formula: qd ~ edu + jiankang + peixun + hunyin + cadcj + age + age2 +      sex + dangyuan + Comp.1 + Comp.2 + trust.cz1 + (trust.cz1 |      commid) 
   Data: individual 
 Family: quasibinomial(logit link)
      AIC      BIC    logLik deviance
 736.7059 821.8267 -349.3529 698.7059
Random effects:
 Groups   Name        Variance Std.Dev. Corr  
 commid   (Intercept) 1.56413  1.25065        
          trust.cz1   0.17922  0.42334  1.000 
 Residual             0.89728  0.94725        
# of obs: 652, groups: commid, 39

Fixed effects:
                  Estimate  Std. Error  DF t value Pr(>|t|)  
(Intercept)    -1.6115e-01  6.7997e-01 637 -0.2370  0.81274  
edu            -5.2585e-02  4.1048e-02 637 -1.2810  0.20064  
jiankang       -9.8243e-01  4.4645e-01 637 -2.2005  0.02813 *
peixun         -4.6307e-01  2.6397e-01 637 -1.7542  0.07988 .
hunyin         -1.2255e-02  2.8151e-01 637 -0.0435  0.96529  
hunyin         -2.7726e-01  1.3846e+00 637 -0.2002  0.84136  
hunyin         -2.9759e-01  8.7180e-01 637 -0.3414  0.73295  
cadcj           2.2366e-01  7.6467e-01 637  0.2925  0.77000  
age             9.3626e-02  4.0390e-02 637  2.3180  0.02076 *
age2           -1.3095e-03  5.5104e-04 637 -2.3763  0.01778 *
sex             3.9188e-01  1.9759e-01 637  1.9833  0.04776 *
dangyuan       -5.2558e-01  5.9091e-01 637 -0.8894  0.37410  
Comp.1          5.2463e-02  1.0309e-01 637  0.5089  0.61100  
Comp.2         -1.5048e-01  1.1435e-01 637 -1.3160  0.18863  
trust.cz1      -8.0709e-01  4.4632e-01 637 -1.8083  0.07103 .
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
There were 11 warnings (use warnings() to see them)
> warnings()
Warning messages:
1: optim or nlminb returned message false convergence (8) 
 in: LMEopt(x = mer, value = cv) 
2: optim or nlminb returned message false convergence (8) 
 in: LMEopt(x = mer, value = cv) 
3: optim or nlminb returned message false convergence (8) 
 in: LMEopt(x = mer, value = cv) 
4: optim or nlminb returned message false convergence (8) 
 in: LMEopt(x = mer, value = cv) 
5: optim or nlminb returned message false convergence (8) 
 in: LMEopt(x = mer, value = cv) 
6: optim or nlminb returned message false convergence (8) 
 in: LMEopt(x = mer, value = cv) 
7: optim or nlminb returned message false convergence (8) 
 in: LMEopt(x = mer, value = cv) 
8: optim or nlminb returned message false convergence (8) 
 in: LMEopt(x = mer, value = cv) 
9: optim or nlminb returned message singular convergence (7) 
 in: LMEopt(x = mer, value = cv) 
10: optim or nlminb returned message singular convergence (7) 
 in: LMEopt(x = mer, value = cv) 
11: optim or nlminb returned message singular convergence (7) 
 in: LMEopt(x = mer, value = cv) 

> version
         _              
platform i386-pc-mingw32
arch     i386           
os       mingw32        
system   i386, mingw32  
status                  
major    2              
minor    2.0            
year     2005           
month    10             
day      06             
svn rev  35749          
language R     	



 				


2005-12-14

------
Deparment of Sociology
Fudan University

My new mail addres is ronggui.huang at gmail.com
Blog:http://sociology.yculblog.com




More information about the R-help mailing list