Re: Korr.matriisin tallentaminen ja ihan muuta

[vastaus aiempaan viestiin]

Kirjoittaja: Seppo Mustonen
Sähköposti:    -
Päiväys: 16.12.2004 16:35

Kiitos Kimmolle korrelaatiomatriisista! Säästyy liioilta vaivoilta,
kun tallettaa aluksi vain alakolmiomatriisin ja täydentää sen sitten
matriisikäskyillä lopulliseen muotoon, kuten Kimmo osoitti.

Reijon vihjeiden perusteella sain selville Thurstonen alkuperäiset
20 laatikkoa. Niiden dimensiot x,y,z on taulukoitu artikkelissa
Kaiser&Horst (1975): A score matrix for Thurstone's box problem.
Multivariate Behavioral Research, sivut 17-25.

x y z       x y z
3 2 1       4 4 1
3 2 2       4 4 2
3 3 1       4 4 3
3 3 2       5 2 1
3 3 3       5 2 2
4 2 1       5 3 2
4 2 2       5 3 3
4 3 1       5 4 1
4 3 2       5 4 2
4 3 3       5 4 3

26 muuttujan luettelo taas löytyy mm. toistaiseksi julkaisemattomasta
Jennrichin artikkelista "Rotation to simple loadings using component
loss functions: The oblique case" ja ne ovat

 1. x           11. y2z=y^2*z      21. ypz=2*(y+z)
 2. y           12. yz2=y*z^2      22. sxy=sqrt(x^2+y^2)
 3. z           13. xdy=x/y        23. sxz=sqrt(x^2+z^2)
 4. xy=x*y      14. ydx=y/x        24. syz=sqrt(y^2+z^2)
 5. xz=x*z      15. xdz=x/z        25. xyz=x*y*z
 6. yz=y*z      16. zdx=z/x        26. sxyz=sqrt(x^2+y^2+z^2)
 7. x2y=x^2*y   17. ydz=y/z
 8. xy2=x*y^2   18. zdy=z/y
 9. x2z=x^2*z   19. xpy=2*(x+y)
10. xz2=x*z^2   20. xpz=2*(x+z)

Laskin näiden tietojen pohjalta korrelaatiomatriisin R26, jonka
pitäisi olla sama kuin Thurstonen kirjassa esitetty. Tulee ottaa
kuitenkin huomioon, että 1940-luvulla po. laskutoimitukset olivat
aikamoisen työn takana. Virheitä tuli helpommin, joten pienet
poikkeamat ovat mahdollisia.

Nyt erot ovat kuitenkin melkoisia jo alkuperäisten muuttujien x,y,z
(taulukko yllä) kohdalla:

MAT LOAD P370(1:3,1:3)
MATRIX P370
///             x        y        z
x         1.00000  0.24000  0.26000
y         0.24000  1.00000  0.35000
z         0.26000  0.35000  1.00000

MAT LOAD R26(1:3,1:3)
MATRIX R26
R(BOXES)
///             x        y        z
x         1.00000  0.25211  0.10169
y         0.25211  1.00000  0.25211
z         0.10169  0.25211  1.00000

Vaikuttaa siltä, että jokin mättää erityisesti z-muuttujassa.
Lisäksi matriisista P370 näkyy, että esim. muuttujien x/y ja y/x
korrelaatiokerroin on tasan -1 (eikä luokkaa -0.95, mitä se on
matriisissa R26). Tämä viittaa siihen, että muuttujien x ja y suhteiden
asemasta kyseessä onkin ao. muuttujien erotukset x-y ja y-x.
Tämä piirre kuvastuu kaikissa näkemissäni po. aineiston eri tavoin
rotatoiduissa faktorimatriiseissa niin, että esim. muuttujilla x/y ja
y/x ovat täsmälleen samat, mutta vastakkaismerkkiset lataukset.
(Sama koskee yhdistelmiä x/z,z/x ja y/z,z/y.)

x/y vs. x-y -sekaannus näkyy hyvin myös seuraavista laskelmista:

Tapaus x/y:
MAT E=P370-R26  / Korrelaatiomatriisien erotus

MAT S=SUM(SUM(E,2)')
MAT_S(1,1)=22.448726718193 / Erotusten neliösumma

MAT S2=SUM(SUM(P370,2)')
MAT_S2(1,1)-26=247.51 / korrelaatiokertoimien neliöiden summa

Matriisit ovat siis karkesti ottaen samat, mutta poikkeamat silti
kohtuuttomia.

Tapaus x-y:
MAT S=SUM(SUM(E,2)')
MAT_S(1,1)=12.467262207771 / Erotusten neliösumma (selvästi pienempi)
MAT_S(1,1)/(26*25)=0.01918040339657 / ero alkiota kohti

Yhteensopivuus on selvästi parempi, mutta jäljelle jää epäilys
siitä, että joko alkuperäinen laatikosto on erilainen kuin mitä
lähteet kertovat tai sitten Thurstone tai hänen apulaisensa
on tehnyt pieniä laskuvirheitä.

Olen saamassa postitse vasta huomenna Thurstonen kirjan. Toivon, että
se antaa lisävalaistusta.

                     - - -

Kimmon viestin kuittaus
- "html kiva moka hei"
lienee seurausta viime seminaaritilaisuudessa käydystä keskustelusta.
LASKE-sukroa koskevan monologini viimeisen viestissä
taisin verrata po. sukron käyttöä taikatemppuihin ja yllytin
survoilijoita hämmästyttämään ystäviään. Viestini kuittaus oli
- On se temppus', on!

Kysyin seminaarissa (joiltakin jo sitä ennen), mitä mahtaa moinen
signeeraus oikein tarkoittaa. Kukaan (???!!!) ei ymmärtääkseni
ilman vihjeitä tajunnut, että kyseessä oli anagrammi. Osasinpas hämätä!

En muistaakseni ole koskaan aikaisemmin tehnyt anagrammeja.
Pari viikkoa sitten keskustellessani puhelimessa musiikista
kiinnostuneen ystäväni kanssa tuli mieleeni mainita, että
säveltäjä Kalevi Aho (Ihka Ovela) on ollut etevä tekemään
kollegoistaan anagrammeja tyyliin

Ei luja bisnes - Jean Sibelius
Kopean salskeena - Esa-Pekka Salonen
Piinajuonne kai - Jouni Kaipainen
Komein köhii - Mikko Heiniö
Aasin löpinät - Tapani Länsiö
Lihakammoko? - Kimmo Hakola

ja eräästä nuoremmasta veikosta hän on laatinut jopa runonpätkän
Melun, loiston/ tunnelmis' olo/ solot unelmin

Koska lupasin lähettää ystävälleni po. jutut sähköpostitse, piti
väsätä kuittaus samassa hengessä ja noin viidentoista minuutin
Survo-avusteisen ähellyksen jälkeen kehkeytyi
- On se temppus', on!

Siitä miettimään, josko kannattaisi Survoon liittää joitain
asiaan liittyviä toimia, mutta onneksi tyttäreni Elina (Lintusen omena
- Suomen tilanne - Onnesi unelmat) tiesi, että verkosta löytyy valmis
Arrak-anagrammi-palvelu
http://www.arrak.fi/ag/ 
joka erittäin nokkelan ohjelman (laatija teekkari A. Gustafsson 1992)
avulla tuottaa nopeasti anagrammeja mielin määrin mistä sanoista
tahansa ja jopa usealla kielellä.
Ohjelma on alunperin C-kielinen ja sen lähdekoodi (vain noin 1500
merkkiä) löytyy mainitun sivun linkkien välityksellä. Se käyttää
valmista esim. verkosta suomenkielisiä sivuja poimimalla saatua ja
tekstitiedostoksi talletettua sanastoa.
Ohjelma on palkittu ("most humourous output") vuoden 1992 kilpailussa
"International Obfuscated C Code Contest".
Todettakoon lopuksi, että anagrammia "On se temppus', on!" tuo ohjelma
ei keksi, sillä sen sanavarastoon ei kuulu runokielen "temppus" eikä
se anna vastauksia, joissa sama sana esiintyy useamman kerran.

- Se pop

Vastaukset:

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!

Etusivu  |  Keskustelu
Copyright © Survo Systems 2001-2013. All rights reserved.
Updated 2013-06-15.