[GenABEL-dev] faster polygenic

[Yurii Aulchenko]

Thanks, will update the code, give it a try in larger # of
simulations, and tell the results!

Yurii

On Mon, Mar 7, 2011 at 4:52 PM, Xia Shen <xia.shen at lcb.uu.se> wrote:
> Hi Yurii,
>
> I just 'sudo-randomly' tried some different argument value in the hglm function call, and happened to have solved the bias issue (at least it seems to be so). Please try putting the following arguments in your hglm wrapper:
>
> hglm(y = ph, X = desmat, Z = L, conv = 1e-6, maxit = 40)
>
> I re-did 100 simulations and the bias is not sig from 0 with p-value = 0.77!
>
> I guess the problem here is just numerical. Since IBS matrix is generally sparse (close to identity), the estimation of variance components will be sensitive numerically. Just make the convergence criteria more strict, and it works.
>
> Note that this might slow down hglm a bit (not very much I guess), but I think it should be more than 20x faster than the old polygenic. Not familiar with fmm, but: 1. our hglm can deal with any GLM family of the trait distribution, even non-normal genetic effects. 2. we produce standard errors for everything - every part is checkable using h-likelihood :)
>
> Cheers
>
> Xia
>