[GenABEL-dev] Interesting idea by Gulnara Svischeva

[Lars Rönnegård]

Hi,
I attach an implementation of an idea proposed by Gulnara. It uses weighted least squares to compute the REML heritability estimate.
The function is found in "REMLrotateFun.R" and some example code is found in "ExampleCode_REMLrotate.R" that uses the data in "RemlEx.RData".

If someone would be interested in trying it (and possibly also suggest improvements), I would appreciate it.
It works for linear mixed models having a polygenic random effect and iid residuals. It seems to be very fast. The example I attach contains 680 observations.
Is this a new algorithm or has it been done already?

Best wishes,
Lars Rönnegård
Dalarna University
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/genabel-devel/attachments/20110518/f7727348/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: REMLrotateFun.R
Type: application/octet-stream
Size: 2350 bytes
Desc: REMLrotateFun.R
URL: <http://lists.r-forge.r-project.org/pipermail/genabel-devel/attachments/20110518/f7727348/attachment-0003.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ExampleCode_REMLrotate.R
Type: application/octet-stream
Size: 293 bytes
Desc: ExampleCode_REMLrotate.R
URL: <http://lists.r-forge.r-project.org/pipermail/genabel-devel/attachments/20110518/f7727348/attachment-0004.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: RemlEx.RData
Type: application/octet-stream
Size: 37570 bytes
Desc: RemlEx.RData
URL: <http://lists.r-forge.r-project.org/pipermail/genabel-devel/attachments/20110518/f7727348/attachment-0005.obj>