[GenABEL-dev] Interesting idea by Gulnara Svischeva

[L.C. Karssen]

Wow! This looks great indeed, good news to start the weekend with :-).


Lennart.


On Fri, 2011-05-20 at 21:16 +0200, Yurii Aulchenko wrote:
> In attached PDF I show the timing for running 'polygenic' (red line;
> ML operating on Gulya's idea), polygenic_hglm (black; Xia/Lars), and
> REML.rotate (green; Lars's implementation/extension of Gulya's idea)
> as a function of sample size (X axes).
> 
> While polygenic and hglm show exponential trend, REMLrotate line is
> not only the lowest, but also looks almost linear to me (also looking
> at the code, it should be roughly linear on time with no of subjects).
> So, with REMLrotate we can get arbitrary big speed-up (x100, x1000,
> x10000, ... you name it!) cf other methods, just by adding more people
> to the data :) Not a bad result!
> 
> Lars and Gulya, you are the champions!
> 
> best wishes,
> Yurii
> 
> 2011/5/19 Yurii Aulchenko <yurii.aulchenko at gmail.com>:
> > Correction: as Lars spotted out, I did placed my 'tags' wrongly
> > concerning time performance results, so actually REML.rotate is the
> > fastest.
> >
> > Here are updated results:
> >
> >> With conv=1e-8 mean time for estimation on 300 IDs was:
> >>
> > 17.81 sec for polygenic_hglm
> >> 2.27 sec for 'polygenic'-ML
> > 1.33 sec for REML.rotate
> >
> > best wishes, and sorry for inconvenience,
> > Yurii
> >
> _______________________________________________
> 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

-- 
-----------------------------------------------
L.C. Karssen
Erasmus MC
Department of Epidemiology
Room Ee 2224

Postbus 2040
3000 CA Rotterdam
The Netherlands

phone: +31-10-7044217
fax: +31-10-7044657
e-mail: l.karssen at erasmusmc.nl
GPG key ID: 0E1D39E3
-----------------------------------------------
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 198 bytes
Desc: This is a digitally signed message part
URL: <http://lists.r-forge.r-project.org/pipermail/genabel-devel/attachments/20110520/58acc36f/attachment.pgp>