[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
Team Biometrie & Kwaliteitszorg / Team Biometrics & Quality Assurance 
thierry.onkelinx using inbo.be
Havenlaan 88 bus 73, 1000 Brussel

///////////////////////////////////////////////////////////////////////////////////////////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>:

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

	[[alternative HTML version deleted]]

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