People
 :bowtie: |
😄 :smile: |
😆 :laughing: |
|---|---|---|
😊 :blush: |
😃 :smiley: |
:relaxed: |
😏 :smirk: |
😍 :heart_eyes: |
😘 :kissing_heart: |
😚 :kissing_closed_eyes: |
😳 :flushed: |
😌 :relieved: |
😆 :satisfied: |
😁 :grin: |
😉 :wink: |
😜 :stuck_out_tongue_winking_eye: |
😝 :stuck_out_tongue_closed_eyes: |
😀 :grinning: |
😗 :kissing: |
😙 :kissing_smiling_eyes: |
😛 :stuck_out_tongue: |
| ## Model: a specific kind of linear mixed model known as "animal model" by geneticists | |
| ## y = mu 1_N + X b + Z u + e = W a + Z u + e | |
| ## y is N x 1; X is N x P; Z is N x Q; W is N x (P+1) | |
| ## u ~ Norm_Q(0, sigma_u^2 A); e ~ Norm_N(0, sigma^2 I_N) | |
| ## Goal of this document: estimate the variance components sigma_u^2 and sigma^2 | |
| ## 1) simulate some data | |
| ## 2) fit the model above using the package rrBLUP (v4.3, on CRAN) | |
| ## 3) fit the model above using lme4 (v1.7, on CRAN) | |
| ## Author: Timothée Flutre (INRA) |
