Startvärden: Varianser och korrelationer · 2007-01-20 19:07 (1 komm.)
Fortsättning på föregående artikel om fleregenskapsanalyser med Asreml, startvärden: Hur (i vilken ordning) skriver man startvärdena i kommandofilen när skattningarna ska baseras på korrelationer i stället för kovarianser?
Som regel är det lättare att hitta bra startvärden ifall man använder en korrelationsmodell (vi vet till exempel att korrelationsvärden håller sig mellan 1.0 och –1.0). Vårt exempel (från Avsevärda arvbarheter) har två rader med startvärden:
...
(1) Trait 0 US 5 0 10 0 0 .25
...
(2) Trait 0 CORGH 0.8 0.5 0.5 .5 .5 .5
...
Rad (2)
avser den genetiska delen i analysen. CORGH
talar om att varians/kovarianskomponenterna för de tre egenskaperna ska skattas med ostrukturerade varianser och korrelationer i stället för med varianser och kovarianser (som i rad (1)
med US
). Egentligen är vi ute efter en likadan 3×3-varians/kovariansmatrismatris, men korrelationer är lättare att arbeta med och de kan enkelt räknas om till kovarianser. Korrelationsmatrisen
1 | r12 | r13 |
r21 | 1 | r23 |
r31 | r32 | 1 |
(r12 = r21, r13 = r31, r23 = r32) med tillhörande varianser [ v1 v2 v3 ] räknas lätt om till en varians/kovariansmatris (som den i föregående artikel).
Det som räknas upp med 0.8 0.5 0.5 .5 .5 .5
är nu korrelationsmatrisens nedre triangel förutom ettorna samt varianserna därunder:
r21 | ||
r31 | r32 | |
v1 | v2 | v3 |
i den ordning de står: r21 r31 r32 v1 v2 v3
. Med matrisuppställning blir det
0.8 | ||
0.5 | 0.5 | |
.5 | .5 | .5 |
På samma sätt som tidigare kan vi byta ut rad (2)
mot de fyra raderna
...
Trait 0 CORGH !+6
0.8
0.5 0.5
.5 .5 .5
...
Och ordningsföljden mellan startvärdena bibehålls i resultatfilens ”components”-tabell.
kommentarer
-
Av Tore 2008-06-27 20:27 #
I och med Asreml 2 (2008) behöver man inte längre ha med antal-värden-på-nästa-rad-specifikarionen (här
!+6
), eftersom de värden som behövs nu automatiskt hittas på efterföljande rader.
kommentera: