Estimating Heritability With Sire-Offspring Regression In R
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
• 3,364 views
•
link
updated
by
Biostar
•
written
by
Recology_ •
•
written
by
Recology_ •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
Hi Scott, Can you please add the link to the book you have mentione in this question ?
done, there's no open source version that I know of though...
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.