Re: Tulosmatriiseista

[vastaus aiempaan viestiin]

Edellisessä viestissäni on virhe, josta Simo Puntanen on ystävällisesti
huomauttanut.
Esitin regressiokertoimien kovarianssien laskemiseksi komentoa
MAT REGCOV=INV(RG.M(MSE))*DV(REGS.M(*,2))*REG_CORR.M*DV(REGS.M(*,2))
Siinä ensimmäinen tekijä INV(RG.M(MSE)) on täysin tarpeeton eli oikea
lauseke on luonnollisesti

MAT REGCOV=DV(REGS.M(*,2))*REG_CORR.M*DV(REGS.M(*,2))

koska ihan yleisestikin korrelaatioista saadaan kovarianssi
kertomalla hajontojen tulolla.
Tällöin REGCOV voidaan laskea (samaan tapaan kuin muuttujien
kovarianssimatriisi CORR-operaation jälkeen) myös sukrokomennolla

/COV REGCOV,REG_CORR.M,REGS.M

ja kovarianssimatriisiksi tulee

LOADM REGCOV,123456.123456789,CUR+1
DV(REGS.M(*,2))*regr(DECA)*DV(REGS.M(*,2))
                 Constant             100m           L_jump
Constant 266579.482774176    -43.534535572    -54.783468026
100m        -43.534535572      0.093856508     -0.013747811
L_jump      -54.783468026     -0.013747811      0.091492077
Shot_put    -57.470517435     -0.012697817      0.001590386
Hi_jump     -76.551077502      0.026277435     -0.007146828
400m        -98.357531812     -0.042305004     -0.007523144

                 Shot_put          Hi_jump             400m
Constant    -57.470517435    -76.551077502    -98.357531812
100m         -0.012697817      0.026277435     -0.042305004
L_jump        0.001590386     -0.007146828     -0.007523144
Shot_put      0.067050560     -0.007744451      0.028546119
Hi_jump      -0.007744451      0.068460798      0.013965032
400m          0.028546119      0.013965032      0.131917574

-Seppo

Luettelo keskustelunaiheista

Survo-keskustelupalstan (2001-2013) viestit arkistoitiin aika ajoin sukrolla, joka automaattisesti rakensi viesteistä (yli 1600 kpl) HTML-muotoisen sivukokonaisuuden. Vuoden 2013 alusta Survo-keskustelua on jatkettu entistäkin aktiivisemmin osoitteessa forum.survo.fi. Tervetuloa mukaan!