Estimating Heritability With Sire-Offspring Regression In R

2

Entering edit mode

12.1 years ago

Recology_ ▴ 100

Has anyone written code in R for estimating heritability with sire offspring regression following the Lynch and Walsh chapter (Chapter 17) using iteratively re-weighted least squares regression?

There is no open source version of the Lynch and Walsh book, but here is a link to it on Amazon: http://www.amazon.com/Genetics-Analysis-Quantitative-Traits-Michael/dp/0878934812

The rlm() function in the MASS R package does iteratively re-weighted least squares regression, but I am not sure how to specify the correct weight.

Thanks, Scott Chamberlain

r • 2.8k views

ADD COMMENT • link updated 4.1 years ago by Biostar 20 • written 12.1 years ago by Recology_ ▴ 100

0

Entering edit mode

Hi Scott, Can you please add the link to the book you have mentione in this question ?

ADD REPLY • link 12.1 years ago by Khader Shameer 18k

0

Entering edit mode

done, there's no open source version that I know of though...

ADD REPLY • link 12.1 years ago by Recology_ ▴ 100

0

Entering edit mode

Hi, 'lm' and 'glm' functions have the argument 'weights'... You can add a vector with the sample size (n) of each 'offspring family', and perform the weighted regression following the Lynch/Walsh recommendation! All the best, Anselmo.

ADD REPLY • link 10.6 years ago by anselmobiologia • 0

Login before adding your answer.