error with glmer
Hi Hank,
I hope to be more clear now.
I am testing the influence of some dental restoration materials
Data, which are positively skewed and clustered, fit the assumption of
being gamma distributed.
I tried glmer()
>mod.glmer=glmer(mis~X_Igruppo_2+X_Igruppo_3+X_Igruppo_4+(1|mon),data=p4,family=Gamma(link='log'))
However,
>summary(mod.glmer)
gives me back the following error message:
Error in asMethod(object) : matrix is not symmetric [1,2]
The same happens when I use lmer() with gamma family. It seems both
glmer() and lmer() don't work when using gamma family. Did someone else
experience this problem?
>sessionInfo()
R version 2.7.2 (2008-08-25)
i386-pc-mingw32
locale:
LC_COLLATE=Italian_Italy.1252;LC_CTYPE=Italian_Italy.1252;LC_MONETARY=Italian_Italy.1252;LC_NUMERIC=C;LC_TIME=Italian_Italy.1252
attached base packages:
[1] stats graphics grDevices datasets tcltk utils
methods base
other attached packages:
[1] lme4_0.999375-20 Matrix_0.999375-10 lattice_0.17-13
xlsReadWrite_1.3.2 svSocket_0.9-42 svMisc_0.9-5
TinnR_1.0.2 R2HTML_1.59
[9] Hmisc_3.4-3
loaded via a namespace (and not attached):
[1] cluster_1.11.11 grid_2.7.2 tools_2.7.2
> 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(mod.glmer)
1: summary(mod.glmer)
Hank Stevens wrote:
Hi Paola, I think you might get more helpful info if you *provide* more info.
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Hank On Dec 29, 2008, at 3:19 AM, paola berchialla wrote:
Dear all, I have some troubles in using lme4 package, glmer function. Here it is the code I am trying to ran unsuccessfully: mod.glmer=glmer(Cmis~X_Igruppo_2+X_Igruppo_3+X_Igruppo_4+(1|mon),data=p4,family=Gamma(link='log')) summary(mod.glmer) The problem is I get the following error message: Error in asMethod(object) : matrix is not symmetric [1,2] Could someone please explain this error to me? I am using R version 2.7.2
_______________________________________________ R-sig-mixed-models at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
Dr. Hank Stevens, Associate Professor 338 Pearson Hall Botany Department Miami University Oxford, OH 45056 Office: (513) 529-4206 Lab: (513) 529-4262 FAX: (513) 529-4243 http://www.cas.muohio.edu/~stevenmh/ http://www.muohio.edu/botany/ http://www.cas.muohio.edu/ecology/ http://www.muohio.edu/ecoreu/ "E Pluribus Unum"
------------------------------------------------------------------------------------ Paola Berchialla, PhD Dept. of Public Health and Microbiology - University of Torino Mobile: +39 3478768532 E-mail: paola.berchialla at unito.it