Heritability (h2) estimation

Open YinLiLin opened this issue 2 years ago • 3 comments

Hi,

Thank you for developing such an useful tool! I am wondering how to obtain the estimated heritability of a disease from the return of the function glmmkin. I tried the command as follows:

> model0$theta
dispersion      kins1 
 1.0000000  0.3377331

Is the value 0.3377331 an estimated heritability of the disease in analysis? or the genetic variance of the disease?

With regards, LL

Jun 19 '23 04:06 YinLiLin

Hi LL,

This is not the estimated heritability as you are using a logit link function.

Best, Han

Jun 30 '23 19:06 hanchenphd

Thanks, so is it possible to derive the estimated heritability from the returns?

Jul 01 '23 01:07 YinLiLin

Unfortunately, it is unclear to me how to derive that based on the logit link function.

Jul 03 '23 15:07 hanchenphd