[GenABEL-dev] faster polygenic

[Yurii Aulchenko]

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.
>