Re: Selittelypuheenvuoro

[vastaus aiempaan viestiin]

Esko,

Nyt esimerkissäsi on jostain syystä pielessä varsinainen
korrelaatiomatriisi R (monet luvuista ovat yli ykkösen)!
RELIAB-ohjelma ei tätä tarkista; se on parasta tehdä itse
esim. ominaisarvohajotelman avulla. Tässä tapauksessa yksi
om.arvoista on negatiivinen, joten R ei ole korr.matriisi.

Tällaisissa asetelmissa lienee yksinkertaisinta tallettaa
pelkkä latausmatriisi (tässä L) ja laskea siitä vastaava R
seuraavasti:

MAT RLABELS X TO L / muuttujille nimet
MAT R=L*L'-DIAG(L*L')+IDN(p,p) / p=10
MAT LOAD R #.### CUR+2 / katsotaan:

MATRIX R
L*L'-DIAG(L*L')+IDN
///         X1    X2    X3    X4    X5    X6    X7    X8    X9   X10
X1       1.000 0.560 0.480 0.400 0.320 0.000 0.000 0.000 0.000 0.000
X2       0.560 1.000 0.420 0.350 0.280 0.000 0.000 0.000 0.000 0.000
X3       0.480 0.420 1.000 0.300 0.240 0.000 0.000 0.000 0.000 0.000
X4       0.400 0.350 0.300 1.000 0.200 0.000 0.000 0.000 0.000 0.000
X5       0.320 0.280 0.240 0.200 1.000 0.000 0.000 0.000 0.000 0.000
X6       0.000 0.000 0.000 0.000 0.000 1.000 0.560 0.480 0.400 0.320
X7       0.000 0.000 0.000 0.000 0.000 0.560 1.000 0.420 0.350 0.280
X8       0.000 0.000 0.000 0.000 0.000 0.480 0.420 1.000 0.300 0.240
X9       0.000 0.000 0.000 0.000 0.000 0.400 0.350 0.300 1.000 0.200
X10      0.000 0.000 0.000 0.000 0.000 0.320 0.280 0.240 0.200 1.000

Oletan seuraavassa, että tarkoitit tällaista asetelmaa.

Asetelma on täysin yksinkertaisen rakenteen (simple structure)
mukainen, ts. muuttujaryhmien (X1--X5) ja (X6--X10) välillä ei
ole lainkaan korrelaatiota. Tilanne vastaa siis täydelleen kahta
ortogonaalista faktoria (jopa keskenään identtistä, kun latauksetkin
ovat aivan samat). Tilanne ei tästä muutu, vaikka faktoreiden
korrelaatiomatriisia muutettaisiin. Eroja alkaa tulla näkyviin vasta,
kun faktoreiden korrelaatiot heijastuvat edes jollain tavoin myös
muuttujien väliseen korrelaatiomatriisiin.

- Kimmo

PS. Ohessa samat reliabiliteettilaskelmat suoraan matriisitulkilla:

MAT RHO=INV(IDN(2,2)+DIAG(W'*(R-L*P*L')*W)*INV(DIAG(W'*L*P*L'*W)))
MAT LOAD RHO CUR+2

MATRIX RHO
INV(IDN+DIAG(W'*(L*L'-DIAG(L*L')+IDN-L*P*L')*W)*INV(DIAG(W'*L*P*L'*W)))
///             1        2
  1      0.743802 0.000000
  2      0.000000 0.743802

• Re: Selittelypuheenvuoro (Esko Kaukonen 22.3.2009 17:49)

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!