[R-sig-ME] Error in asMethod(object) : matrix is not symmetric [1, 2]

Chaudhari, Monica mchaudhari at deltadentalwa.com
Wed Dec 1 00:57:10 CET 2010


Hello All,

I am receiving this error while trying to model follow-up  restorative
costs on  patients' teeth using Gamma family mixed model. 

> lmer9<-lmer(cumappAmtCPIAdj ~ toothNo  + yrFromBase +  (yrFromBase |
ptntId/quad), data = datsample[procClass == "Restorative" & ptntId %in%
levels(ptntId)[1:10],],family = Gamma(link="log"))
The model appears to fit until you give the bellow command:
> summary(lmer9)
Error in asMethod(object) : matrix is not symmetric [1,2]

When tracing back for error, it gives:
> traceback()
15: .Call(dense_to_symmetric, from, "U", TRUE)
14: asMethod(object)
13: as(from, "symmetricMatrix")
12: .class1(object)
11: as(as(from, "symmetricMatrix"), "dMatrix")
10: .class1(object)
9: as(as(as(from, "symmetricMatrix"), "dMatrix"), "denseMatrix")
8: .class1(object)
7: as(as(as(as(from, "symmetricMatrix"), "dMatrix"), "denseMatrix"), 
       "dpoMatrix")
6: asMethod(object)
5: as(sigma(object)^2 * chol2inv(object at RX, size = object at dims["p"]), 
       "dpoMatrix")
4: vcov(object)
3: vcov(object)
2: summary(lmer9)
1: summary(lmer9)

Below is additional information about the session:
> sessionInfo()
R version 2.11.1 (2010-05-31) 
i386-pc-mingw32 

locale:
[1] LC_COLLATE=English_United States.1252  LC_CTYPE=English_United
States.1252    LC_MONETARY=English_United States.1252 LC_NUMERIC=C
LC_TIME=English_United States.1252    

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base


other attached packages:
[1] lme4_0.999375-37   Matrix_0.999375-45 lattice_0.19-13   

loaded via a namespace (and not attached):
[1] grid_2.11.1   nlme_3.1-97   stats4_2.11.1

I have not had success fitting Gamma(link="log") even a single time
using lme4. It returns the above error always. Could you please point
any resolutions?

Thanks,
Monica

______________________________________________________________________
The information contained in this e-mail and subsequent ...{{dropped:9}}




More information about the R-sig-mixed-models mailing list