[R-sig-ME] error with glmer
paola berchialla
paola.berchialla at unito.it
Mon Dec 29 22:49:18 CET 2008
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
More information about the R-sig-mixed-models
mailing list