Söktext
 Upp   Ned   Mittemellan 

CORGH  · 2007-02-08 17:37   (2 komm.)

Idag upptäckte jag alldeles oväntat en del om varför det i hela friden ska stå H i CORGH – sådär bara – här är den historien (endast för Asrumlare, tillämpningen är genetisk fleregenskapsanalys):

En bit ned på sidan 120 i Asreml 2.0 user guide
står det vid

CORG   general correlation   CORGH=US

bland annat

number of parameters
homo’s variance:
w(w – 1)/2 + 1

och

number of parameters
hetero’s variance:
w(w – 1)/2 + w

[ w står för antalet variabler (egenskaper) ]

Denna kryptiska framställning kan man kanske smälta förutsatt att man först läst manualen ordentligt innan man börjar läsa den. Troligtvis hänger det ihop så här:

CORG med tillbehöret V, det vill säga CORGV, ska användas då man vill att Asreml ska skatta korrelationerna mellan variablerna samt en gemensam varians för samtliga variabler. Tillägget V beskrivs på sidan 124 (här med EXPV) och hör ihop med homo’s variance ovan och med parameteruppställning w(w – 1)/2 stycken korrelationer + 1 varians (den gemensamma).

CORG med tillbehöret H, det vill säga CORGH, ska användas då man vill att Asreml ska skatta korrelationerna mellan variablerna samt olika varianser för variablerna. Tillägget H beskrivs också på sidan 124 (här med CORUH) och hör ihop med hetero’s variance ovan och med parameteruppställning w(w – 1)/2 stycken korrelationer + w olika varianser.

Den fantasifulle kan kanske ana att det finns fall då man inte vill skatta några varianser alls, vilket i så fall skulle motsvara CORG utan tillägg där de uteblivna variansskattningarna motsvarar korrelationsmatrisens diagonal som är förutbestämd att bestå av ettor.

OK?


 
avdelning: Datorbruk   ämne:  

Startvärden: Varianser och korrelationer Skattning av genetiska parametrar med Asreml


kommentarer

  1. Av  Tore  2007-02-09 16:33    #

    Man kan kanske tänka så här om CORxy:

    x är för korrelationsstrukturen och y är för variansstrukturen. Korrelationerna först och varianserna sist. G = generell korrelationsmatris utan restriktioner och H = ostrukturerade varianser utan restriktioner.

    (Därför kommentaren CORGH=US, det vill säga hela matrisen utan restriktioner på samma sätt som med US för en varians/kovariansmatris)

  2. Av  Tore  2008-10-18 16:46    #

    Hade nyss anledning att testa CORGV (i stället för CORGH) men i det speciella fallet fungerade det inte, men CORGH och !=00...AA... gick bra i stället (nollor betyder ingen restriktion för korrelationerna som kommer först och A:n för alla varianserna lika därefter)


kommentera:

 
     
 
 
 
Upp
Om du misstänker att jag inte fått din kommentar: Skicka samma kommentar igen! Jag vill gärna visa din kommentar om den rör artikeln men har inte själv kontroll över att alla servrar fungerar... / Tore