Include Cov(\sigma^{2}) in output for variance component model (w/ sum stats)

genetics-statistics / GEMMA

Genome-wide Efficient Mixed Model Association

https://github.com/genetics-statistics/GEMMA

GNU General Public License v3.0

333 stars 125 forks source link

Include Cov(\sigma^{2}) in output for variance component model (w/ sum stats) #209

Open julirsch opened 5 years ago

julirsch commented 5 years ago

Hi,

Would it be possible to include Cov(\sigma^{2}) (and not just the variances) from the sum stats version of the variance component model in the output. Alternatively, I'd be happy to fork and implement myself if you could point me to where in the code this is calculated.

Best, Jacob

pjotrp commented 5 years ago

Sure, feel free. The code is not hard to find, just search for the other columns.

julirsch commented 5 years ago

I looked through and struggled -- could you point me to which specific variable this it?

pjotrp commented 5 years ago

Look for v_sigma2 in gemma.cpp