Temat: Z Fishera dla różnic między korelacjami r Pearsona -...

Witam Wszystkich!

Czy mógłby mi ktoś wskazać sposób na przeprowadzenie Z-Fishera w celu porównania korelacji r Perasona dla dwóch grup w SPSS?

Z góry dziękuję za odpowiedź!

Temat: Z Fishera dla różnic między korelacjami r Pearsona -...

Witam,

tak na szybko znalazłem na necie kod:

* Independent populations: Fisher's Z test *. 

DEFINE TWOCORRS (r1=!TOKENS(1)/n1=!TOKENS(1)/r2=!TOKENS(1)/n2=!TOKENS(1)). 

MATRIX. 

PRINT /TITLE='COMPARING TWO INDEPENDENT CORRELATION COEFFICIENTS'. 

COMPUTE r1=!r1. 

COMPUTE n1=!n1. 

COMPUTE r2=!r2. 

COMPUTE n2=!n2. 

PRINT {r1,r2} 

 /FORMAT='F8.2' 

 /CLABELS='R1' 'R2' 

 /TITLE='Correlation values being compared:'. 

PRINT {n1,n2} 

 /FORMAT='F8.0' 

 /CLABELS='N1','N2' 

 /TITLE='Sample sizes:'. 

COMPUTE zr1=.5*(LN((1+r1)/(1-r1))). 

COMPUTE varzr1=1/(n1-3). 

COMPUTE zr2=.5*(LN((1+r2)/(1-r2))). 

COMPUTE varzr2=1/(n2-3). 

COMPUTE z=ABS(zr1-zr2)/SQRT(varzr1+varzr2). 

COMPUTE z_sig=2*(1-CDFNORM(z)). 

PRINT {z,z_sig,(z_sig/2)} 

 /FORMAT='F8.3' 

 /CLABELS='Z' '2-tails' '1-tail' 

 /TITLE='Z value and significance (2-tailed & 1-tailed)'. 

END MATRIX. 

!ENDDEFINE. 



*odpalasz przykład podając wartości współczynników korelacji i liczebności

TWOCORRS r1=0.397 n1=76 r2=0.819 n2=64. 

źródło: http://spssx-discussion.1045642.n5.nabble.com/comparin...

edit: proszę sprawdź poprawność.Maciej B. edytował(a) ten post dnia 27.04.12 o godzinie 13:39

Temat: Z Fishera dla różnic między korelacjami r Pearsona -...

nie wiem, jak to liczy Statistica 6, ale wyniki nie są spójne.

Natomiast procedura ta jest spójna z tym, co pisze Góralski, A. (1987). Metody opisu i wnioskowania statystycznego w psychologii i pedagogice. Warszawa: Państwowe Wydawnictwo Naukowe, ss. 221-222.

Temat: Z Fishera dla różnic między korelacjami r Pearsona -...

Według helpa Statistici (Inne testy istotności) wzór na statystykę testową jest następujący:

t=(r'1-r'2)/  sqrt((n1+n2-6)/ ((n1-3)*(n2-3)))

i ma rozkład t-studenta z n1+n2-4 stopniami swobody. Gdzie r' to współczynnik korelacji przekształcony według wzoru (taki sam jak w synatxie):

r'=0.5*[ln(1+r) - ln(1-r)]

a w kodzie, który zalazłem (i Pan Andrzej Wichrowski potwierdził w literaturze) jest następujący:

z=|r'1-r'2| / sqrt( 1/(n1-3)+1/(n2-3))=|r'1-r'2| / sqrt((n1+n2-6)/ ((n1-3)*(n2-3)))

i ma rozkład normalny.

Stąd (inny asymptotyczny rozkład) mogą wynikać rozbieżności w wynikach ale decyzja wciąż jest taka sama.Maciej B. edytował(a) ten post dnia 30.04.12 o godzinie 17:58