[R-sig-ME] Force the random effects variances to be the same
piuo244 m@iii@g oii @ucki@@du@i@@c@@z
piuo244 m@iii@g oii @ucki@@du@i@@c@@z
Fri May 13 11:17:04 CEST 2022
Hi Thierry,
Because I would like to fit a heritability model, where the genetic variance should be the same for mother and father. I was doing this with lme4qtl which allows correlated random effects. In my case, the correlation matrix is a block diagonal matrix, and a block looks like
1 0 1/2 1/2 1/20 1 1/2 1/2 1/21/2 1/2 1 1/2 1/21/2 1/2 1/2 1 1/2
1/2 1/2 1/2 1/2 1
But the performance of the Laplace approximation is quite poor for binary data with small cluster sizes. I think letting nAGQ=3 might help, but lme4qtl does not allow me to do so.
Cheers,Zoe On Friday, May 13, 2022, 07:37:07 PM GMT+12, Thierry Onkelinx <thierry.onkelinx using inbo.be> wrote:
Why would you want to do that?
ir. Thierry Onkelinx
Statisticus / Statistician
Vlaamse Overheid / Government of Flanders
INSTITUUT VOOR NATUUR- EN BOSONDERZOEK / RESEARCH INSTITUTE FOR NATURE AND FOREST
Team Biometrie & Kwaliteitszorg / Team Biometrics & Quality Assurance
thierry.onkelinx using inbo.be
Havenlaan 88 bus 73, 1000 Brussel
www.inbo.be
///////////////////////////////////////////////////////////////////////////////////////////To call in the statistician after the experiment is done may be no more than asking him to perform a post-mortem examination: he may be able to say what the experiment died of. ~ Sir Ronald Aylmer Fisher
The plural of anecdote is not data. ~ Roger Brinner
The combination of some data and an aching desire for an answer does not ensure that a reasonable answer can be extracted from a given body of data. ~ John Tukey
///////////////////////////////////////////////////////////////////////////////////////////
Op do 12 mei 2022 om 23:29 schreef pluo244 using aucklanduni.ac.nz <pluo244 using aucklanduni.ac.nz>:
Hi,
I would like to use lme4::glmer to fit a model of the form
y ~ x + (0+m|group) + (0+f|group)
where a group looks like
m f1 00 11/2 1/21/2 1/2
1/2 1/2
I was wondering if it is possible to force the variances of the random effects to be the same?
Many thanks,Zoe
[[alternative HTML version deleted]]
_______________________________________________
R-sig-mixed-models using r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
[[alternative HTML version deleted]]
More information about the R-sig-mixed-models
mailing list