[GenABEL-dev] faster polygenic

[Xia Shen]

This looks weird... I need to dig into the code to figure out why hglm did wrong. Will come back to you when having some answers.

Xia

On Feb 28, 2011, at 12:24 AM, Yurii Aulchenko wrote:

> some extra info:
> 
> was running the checks cf. simulated heritability
> 
> t-test
> 
> difference (h2simulated - h2estimated-old-polygenic) ns from 0, p=0.36
> 
> difference (h2simulated - h2estimated-hglm) SIG from 0, p=0.006
> 
> difference (h2simulated - h2estimated-fmm) ns from 0, p=0.30
> 
> So, my feeling is that hglm does something wrong; old and fmm do very
> well, but old does better :)
> 
> On Sun, Feb 27, 2011 at 10:51 PM, Yurii Aulchenko
> <yurii.aulchenko at gmail.com> wrote:
>> Dear All,
>> 
>> Some of you may know, that MixABEL::FastMixedModel by William Astle
>> first does polygenic model estimation (e.g. h2-estimate is given by
>> FastMixedModel(...)$null.herit), though it is build 'within' GWA and
>> it is not yet possible to get hold on all necessary 'null'
>> heritability estimates.
>> 
> _______________________________________________
> genabel-devel mailing list
> genabel-devel at lists.r-forge.r-project.org
> https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/genabel-devel