projects / pspp / blobdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | inline | side by side
CROSSTABS: Add tests for Pearson's R, and update documentation.
[pspp] / tests / language / stats / crosstabs.at
diff --git a/tests/language/stats/crosstabs.at b/tests/language/stats/crosstabs.at
index 4492e9a25002be3cbc73a4a0a96bfa5fdb9b501e..4dbd0bdc78e3ea554e45872007b745955d5c14e7 100644 (file)
--- a/tests/language/stats/crosstabs.at
+++ b/tests/language/stats/crosstabs.at
@@ -175,9 +175,9 @@ Total,2.00,2.00,4.00
 
 Table: Chi-square tests.
 Statistic,Value,df,Asymp. Sig. (2-tailed)
-Pearson Chi-Square,2.00,2,.37
-Likelihood Ratio,2.77,2,.25
-Linear-by-Linear Association,.27,1,.60
+Pearson Chi-Square,2.00,2,.368
+Likelihood Ratio,2.77,2,.250
+Linear-by-Linear Association,.27,1,.602
 N of Valid Cases,4,,
 ]])
 AT_CLEANUP
@@ -207,6 +207,7 @@ CROSSTABS
     /STATISTICS=CHISQ
     /CELLS=COUNT ROW COLUMN TOTAL.
 ])
+
 AT_CHECK([pspp -O format=csv crosstabs.sps], [0],
   [[Variable,Value,Label
 v0,a ,
@@ -235,10 +236,10 @@ Total,4.00,2.00,6.00
 
 Table: Chi-square tests.
 Statistic,Value,df,Asymp. Sig. (2-tailed),Exact Sig. (2-tailed),Exact Sig. (1-tailed)
-Pearson Chi-Square,.38,1,.54,,
-Likelihood Ratio,.37,1,.54,,
-Fisher's Exact Test,,,,1.00,.60
-Continuity Correction,.00,1,1.00,,
+Pearson Chi-Square,.38,1,.540,,
+Likelihood Ratio,.37,1,.545,,
+Fisher's Exact Test,,,,1.000,.600
+Continuity Correction,.00,1,1.000,,
 N of Valid Cases,6,,,,
 
 Variable,Value,Label
@@ -268,10 +269,10 @@ Total,1.00,3.00,4.00
 
 Table: Chi-square tests.
 Statistic,Value,df,Asymp. Sig. (2-tailed),Exact Sig. (2-tailed),Exact Sig. (1-tailed)
-Pearson Chi-Square,.44,1,.50,,
-Likelihood Ratio,.68,1,.41,,
-Fisher's Exact Test,,,,1.00,.75
-Continuity Correction,.00,1,1.00,,
+Pearson Chi-Square,.44,1,.505,,
+Likelihood Ratio,.68,1,.410,,
+Fisher's Exact Test,,,,1.000,.750
+Continuity Correction,.00,1,1.000,,
 N of Valid Cases,4,,,,
 ]])
 AT_CLEANUP
@@ -343,13 +344,13 @@ Total,,3.00,1.00,4.00
 
 Table: Chi-square tests.
 z,Statistic,Value,df,Asymp. Sig. (2-tailed)
-1,Pearson Chi-Square,5.00,4,.29
-,Likelihood Ratio,5.00,4,.29
-,Linear-by-Linear Association,.01,1,.94
+1,Pearson Chi-Square,5.00,4,.287
+,Likelihood Ratio,5.00,4,.287
+,Linear-by-Linear Association,.01,1,.938
 ,N of Valid Cases,5,,
-2,Pearson Chi-Square,4.00,3,.26
-,Likelihood Ratio,4.50,3,.21
-,Linear-by-Linear Association,1.58,1,.21
+2,Pearson Chi-Square,4.00,3,.261
+,Likelihood Ratio,4.50,3,.212
+,Linear-by-Linear Association,1.58,1,.209
 ,N of Valid Cases,4,,
 
 Table: Symmetric measures.
@@ -462,3 +463,494 @@ X1 * X2,0,0.0%,1,100.0%,1,100.0%
 crosstabs.sps:8: warning: CROSSTABS: Crosstabulation X1 * X2 contained no non-missing cases.
 ])
 AT_CLEANUP
+
+
+
+dnl This example comes from http://www.ats.ucla.edu/stat/spss/whatstat/whatstat.htm#chisq
+AT_SETUP([CROSSTABS Fisher Exact Test])
+
+AT_DATA([fisher-exact.sps], [dnl
+SET FORMAT F12.3.
+SET DECIMAL DOT.
+
+DATA LIST notable LIST  /schtyp (F9.2) female (F9.2) ses (F9.2) .
+begin data.
+      1.00       .00      1.00 
+      1.00      1.00      2.00 
+      1.00       .00      3.00 
+      1.00       .00      3.00 
+      1.00       .00      2.00 
+      1.00       .00      2.00 
+      1.00       .00      2.00 
+      1.00       .00      2.00 
+      1.00       .00      2.00 
+      1.00       .00      2.00 
+      1.00       .00      2.00 
+      1.00       .00      2.00 
+      1.00       .00      3.00 
+      1.00       .00      3.00 
+      1.00       .00      1.00 
+      1.00       .00      1.00 
+      1.00       .00      3.00 
+      2.00       .00      2.00 
+      1.00       .00      3.00 
+      1.00       .00      2.00 
+      1.00       .00      2.00 
+      1.00       .00      2.00 
+      1.00       .00      2.00 
+      1.00       .00      3.00 
+      1.00       .00      2.00 
+      1.00       .00      2.00 
+      1.00       .00      3.00 
+      2.00       .00      2.00 
+      2.00       .00      3.00 
+      1.00       .00      1.00 
+      1.00       .00      2.00 
+      1.00       .00      3.00 
+      2.00       .00      3.00 
+      1.00       .00      2.00 
+      2.00       .00      3.00 
+      1.00       .00      3.00 
+      2.00       .00      2.00 
+      1.00       .00      3.00 
+      1.00       .00      1.00 
+      1.00       .00      2.00 
+      2.00       .00      2.00 
+      2.00       .00      2.00 
+      1.00       .00      2.00 
+      1.00       .00      1.00 
+      1.00       .00      3.00 
+      1.00       .00      1.00 
+      1.00       .00      3.00 
+      1.00       .00      2.00 
+      2.00       .00      2.00 
+      1.00       .00      2.00 
+      1.00       .00      2.00 
+      1.00       .00      3.00 
+      1.00       .00      2.00 
+      2.00       .00      2.00 
+      1.00       .00      2.00 
+      1.00       .00      3.00 
+      1.00       .00      1.00 
+      1.00       .00      2.00 
+      2.00       .00      2.00 
+      1.00       .00      2.00 
+      2.00       .00      2.00 
+      1.00       .00      3.00 
+      1.00       .00      1.00 
+      1.00       .00      2.00 
+      2.00       .00      3.00 
+      1.00       .00      2.00 
+      1.00       .00      2.00 
+      1.00       .00      1.00 
+      1.00       .00      1.00 
+      1.00       .00      2.00 
+      1.00       .00      2.00 
+      1.00       .00      3.00 
+      1.00       .00      2.00 
+      1.00       .00      2.00 
+      1.00       .00      2.00 
+      1.00       .00      1.00 
+      1.00       .00      3.00 
+      1.00       .00      3.00 
+      1.00       .00      2.00 
+      1.00       .00      3.00 
+      1.00       .00      3.00 
+      1.00       .00      1.00 
+      2.00       .00      2.00 
+      1.00       .00      1.00 
+      1.00       .00      2.00 
+      1.00       .00      3.00 
+      1.00       .00      3.00 
+      1.00       .00      3.00 
+      1.00       .00      2.00 
+      1.00       .00      3.00 
+      1.00       .00      2.00 
+      1.00       .00      1.00 
+      1.00      1.00      3.00 
+      1.00      1.00      1.00 
+      1.00      1.00      1.00 
+      1.00      1.00      1.00 
+      1.00      1.00      2.00 
+      1.00      1.00      3.00 
+      1.00      1.00      1.00 
+      2.00      1.00      3.00 
+      1.00      1.00      3.00 
+      1.00      1.00      3.00 
+      1.00      1.00      1.00 
+      1.00      1.00      3.00 
+      1.00      1.00      2.00 
+      1.00      1.00      2.00 
+      1.00      1.00      3.00 
+      1.00      1.00      1.00 
+      2.00      1.00      1.00 
+      2.00      1.00      3.00 
+      1.00      1.00      2.00 
+      1.00      1.00      1.00 
+      1.00      1.00      3.00 
+      1.00      1.00      1.00 
+      2.00      1.00      3.00 
+      1.00      1.00      2.00 
+      1.00      1.00      3.00 
+      1.00      1.00      3.00 
+      1.00      1.00      1.00 
+      1.00      1.00      1.00 
+      2.00      1.00      1.00 
+      1.00      1.00      2.00 
+      1.00      1.00      2.00 
+      1.00      1.00      2.00 
+      1.00      1.00      1.00 
+      1.00      1.00      3.00 
+      1.00      1.00      2.00 
+      1.00      1.00      2.00 
+      1.00      1.00      3.00 
+      1.00      1.00      1.00 
+      1.00      1.00      2.00 
+      1.00      1.00      1.00 
+      1.00      1.00      2.00 
+      1.00      1.00      2.00 
+      1.00      1.00      1.00 
+      1.00      1.00      3.00 
+      2.00      1.00      2.00 
+      1.00      1.00      2.00 
+      1.00      1.00      2.00 
+      2.00      1.00      2.00 
+      1.00      1.00      1.00 
+      1.00      1.00      3.00 
+      1.00      1.00      2.00 
+      1.00      1.00      2.00 
+      1.00      1.00      2.00 
+      2.00      1.00      3.00 
+      1.00      1.00      2.00 
+      2.00      1.00      2.00 
+      1.00      1.00      1.00 
+      1.00      1.00      1.00 
+      1.00      1.00      1.00 
+      1.00      1.00      3.00 
+      1.00      1.00      2.00 
+      1.00      1.00      2.00 
+      1.00      1.00      2.00 
+      1.00      1.00      2.00 
+      1.00      1.00      2.00 
+      1.00      1.00      2.00 
+      1.00      1.00      2.00 
+      1.00      1.00      3.00 
+      1.00      1.00      1.00 
+      1.00      1.00      2.00 
+      2.00      1.00      3.00 
+      1.00      1.00      1.00 
+      1.00      1.00      2.00 
+      1.00      1.00      1.00 
+      1.00      1.00      2.00 
+      1.00      1.00      1.00 
+      2.00      1.00      2.00 
+      1.00      1.00      1.00 
+      1.00      1.00      1.00 
+      1.00      1.00      2.00 
+      1.00      1.00      3.00 
+      1.00      1.00      3.00 
+      1.00      1.00      1.00 
+      1.00      1.00      1.00 
+      1.00      1.00      2.00 
+      1.00      1.00      2.00 
+      1.00      1.00      3.00 
+      1.00      1.00      1.00 
+      1.00      1.00      2.00 
+      2.00      1.00      2.00 
+      1.00      1.00      3.00 
+      1.00      1.00      2.00 
+      1.00      1.00      3.00 
+      1.00      1.00      1.00 
+      1.00      1.00      2.00 
+      1.00      1.00      2.00 
+      2.00      1.00      3.00 
+      1.00      1.00      1.00 
+      1.00      1.00      1.00 
+      2.00      1.00      3.00 
+      2.00      1.00      2.00 
+      1.00      1.00      3.00 
+      2.00      1.00      2.00 
+      2.00      1.00      2.00 
+      1.00      1.00      2.00 
+      2.00      1.00      2.00 
+      1.00      1.00      2.00 
+      1.00      1.00      3.00 
+end data.
+
+VARIABLE LABEL schtyp 'type of school'.
+ADD VALUE LABELS female 0 male 1 female.
+ADD VALUE LABELS ses 1 low 2 middle 3 high.
+ADD VALUE LABELS schtyp 1 public 2 private.
+
+crosstabs /tables = schtyp by female /statistic = chisq.
+crosstabs /tables = female by ses  /statistic = chisq.
+])
+
+AT_CHECK([pspp -O format=csv fisher-exact.sps], [0], [dnl
+Table: Summary.
+,Cases,,,,,
+,Valid,,Missing,,Total,
+,N,Percent,N,Percent,N,Percent
+type of school * female,200,100.0%,0,0.0%,200,100.0%
+
+Table: type of school * female [[count]].
+,female,,
+type of school,male,female,Total
+public,77.000,91.000,168.000
+private,14.000,18.000,32.000
+Total,91.000,109.000,200.000
+
+Table: Chi-square tests.
+Statistic,Value,df,Asymp. Sig. (2-tailed),Exact Sig. (2-tailed),Exact Sig. (1-tailed)
+Pearson Chi-Square,.047,1,.828,,
+Likelihood Ratio,.047,1,.828,,
+Fisher's Exact Test,,,,.849,.492
+Continuity Correction,.001,1,.981,,
+Linear-by-Linear Association,.047,1,.829,,
+N of Valid Cases,200,,,,
+
+Table: Summary.
+,Cases,,,,,
+,Valid,,Missing,,Total,
+,N,Percent,N,Percent,N,Percent
+female * ses,200,100.0%,0,0.0%,200,100.0%
+
+Table: female * ses [[count]].
+,ses,,,
+female,low,middle,high,Total
+male,15.000,47.000,29.000,91.000
+female,32.000,48.000,29.000,109.000
+Total,47.000,95.000,58.000,200.000
+
+Table: Chi-square tests.
+Statistic,Value,df,Asymp. Sig. (2-tailed)
+Pearson Chi-Square,4.577,2,.101
+Likelihood Ratio,4.679,2,.096
+Linear-by-Linear Association,3.110,1,.078
+N of Valid Cases,200,,
+])
+
+AT_CLEANUP
+
+AT_SETUP([CROSSTABS Pearson's R])
+# Test 1.
+AT_DATA([pearson.sps], [dnl
+* From http://www.statisticslectures.com/topics/pearsonr/.
+DATA LIST FREE/x y.
+BEGIN DATA.
+1 4
+3 6
+5 10
+5 12
+6 13
+END DATA.
+CROSSTABS x BY y/STATISTICS=CORR.
+])
+AT_CHECK([pspp -O format=csv pearson.sps], [0], [dnl
+Table: Summary.
+,Cases,,,,,
+,Valid,,Missing,,Total,
+,N,Percent,N,Percent,N,Percent
+x * y,5,100.0%,0,0.0%,5,100.0%
+
+Table: x * y [[count]].
+,y,,,,,
+x,4.00,6.00,10.00,12.00,13.00,Total
+1.00,1.00,.00,.00,.00,.00,1.00
+3.00,.00,1.00,.00,.00,.00,1.00
+5.00,.00,.00,1.00,1.00,.00,2.00
+6.00,.00,.00,.00,.00,1.00,1.00
+Total,1.00,1.00,1.00,1.00,1.00,5.00
+
+Table: Symmetric measures.
+Category,Statistic,Value,Asymp. Std. Error,Approx. T,Approx. Sig.
+Ordinal by Ordinal,Spearman Correlation,.97,.02,2.84,
+Interval by Interval,Pearson's R,.97,.02,3.36,
+N of Valid Cases,,5,,,
+])
+
+# Test 2.
+AT_DATA([pearson2.sps], [dnl
+* Checked with http://www.socscistatistics.com/tests/pearson/Default2.aspx.
+DATA LIST FREE/x y.
+BEGIN DATA.
+1 1.5
+2 1.5
+3 4
+4 6
+5 5
+6 7
+7 6.5
+8 9
+9 10.5
+10 11
+END DATA.
+CROSSTABS x BY y/STATISTICS=CORR.
+])
+AT_CHECK([pspp -O format=csv pearson2.sps], [0], [dnl
+Table: Summary.
+,Cases,,,,,
+,Valid,,Missing,,Total,
+,N,Percent,N,Percent,N,Percent
+x * y,10,100.0%,0,0.0%,10,100.0%
+
+Table: x * y [[count]].
+,y,,,,,,,,,
+x,1.50,4.00,5.00,6.00,6.50,7.00,9.00,10.50,11.00,Total
+1.00,1.00,.00,.00,.00,.00,.00,.00,.00,.00,1.00
+2.00,1.00,.00,.00,.00,.00,.00,.00,.00,.00,1.00
+3.00,.00,1.00,.00,.00,.00,.00,.00,.00,.00,1.00
+4.00,.00,.00,.00,1.00,.00,.00,.00,.00,.00,1.00
+5.00,.00,.00,1.00,.00,.00,.00,.00,.00,.00,1.00
+6.00,.00,.00,.00,.00,.00,1.00,.00,.00,.00,1.00
+7.00,.00,.00,.00,.00,1.00,.00,.00,.00,.00,1.00
+8.00,.00,.00,.00,.00,.00,.00,1.00,.00,.00,1.00
+9.00,.00,.00,.00,.00,.00,.00,.00,1.00,.00,1.00
+10.00,.00,.00,.00,.00,.00,.00,.00,.00,1.00,1.00
+Total,2.00,1.00,1.00,1.00,1.00,1.00,1.00,1.00,1.00,10.00
+
+Table: Symmetric measures.
+Category,Statistic,Value,Asymp. Std. Error,Approx. T,Approx. Sig.
+Ordinal by Ordinal,Spearman Correlation,.97,.02,3.66,
+Interval by Interval,Pearson's R,.97,.02,3.69,
+N of Valid Cases,,10,,,
+])
+
+# Test 3.
+AT_DATA([pearson3.sps], [dnl
+* From http://learntech.uwe.ac.uk/da/Default.aspx?pageid=1442.
+DATA LIST FREE/x y.
+BEGIN DATA.
+56 87
+56 91
+65 85
+65 91
+50 75
+25 28
+87 122
+44 66
+35 58
+END DATA.
+CROSSTABS x BY y/STATISTICS=CORR.
+])
+AT_CHECK([pspp -O format=csv pearson3.sps], [0], [dnl
+Table: Summary.
+,Cases,,,,,
+,Valid,,Missing,,Total,
+,N,Percent,N,Percent,N,Percent
+x * y,9,100.0%,0,0.0%,9,100.0%
+
+Table: x * y [[count]].
+,y,,,,,,,,
+x,28.00,58.00,66.00,75.00,85.00,87.00,91.00,122.00,Total
+25.00,1.00,.00,.00,.00,.00,.00,.00,.00,1.00
+35.00,.00,1.00,.00,.00,.00,.00,.00,.00,1.00
+44.00,.00,.00,1.00,.00,.00,.00,.00,.00,1.00
+50.00,.00,.00,.00,1.00,.00,.00,.00,.00,1.00
+56.00,.00,.00,.00,.00,.00,1.00,1.00,.00,2.00
+65.00,.00,.00,.00,.00,1.00,.00,1.00,.00,2.00
+87.00,.00,.00,.00,.00,.00,.00,.00,1.00,1.00
+Total,1.00,1.00,1.00,1.00,1.00,1.00,2.00,1.00,9.00
+
+Table: Symmetric measures.
+Category,Statistic,Value,Asymp. Std. Error,Approx. T,Approx. Sig.
+Ordinal by Ordinal,Spearman Correlation,.91,.07,3.45,
+Interval by Interval,Pearson's R,.97,.02,5.00,
+N of Valid Cases,,9,,,
+])
+
+# Test 4.
+AT_DATA([pearson4.sps], [dnl
+* From http://psychology.ucdavis.edu/faculty_sites/sommerb/sommerdemo/correlation/hand/pearson_hand.htm.
+DATA LIST FREE/x y.
+BEGIN DATA.
+5 5
+10 20
+6 4
+8 15
+4 11
+4 9
+3 12
+10 18
+2 7
+6 2
+7 14
+9 17
+END DATA.
+CROSSTABS x BY y/STATISTICS=CORR.
+])
+AT_CHECK([pspp -O format=csv pearson4.sps], [0], [dnl
+Table: Summary.
+,Cases,,,,,
+,Valid,,Missing,,Total,
+,N,Percent,N,Percent,N,Percent
+x * y,12,100.0%,0,0.0%,12,100.0%
+
+Table: x * y [[count]].
+,y,,,,,,,,,,,,
+x,2.00,4.00,5.00,7.00,9.00,11.00,12.00,14.00,15.00,17.00,18.00,20.00,Total
+2.00,.00,.00,.00,1.00,.00,.00,.00,.00,.00,.00,.00,.00,1.00
+3.00,.00,.00,.00,.00,.00,.00,1.00,.00,.00,.00,.00,.00,1.00
+4.00,.00,.00,.00,.00,1.00,1.00,.00,.00,.00,.00,.00,.00,2.00
+5.00,.00,.00,1.00,.00,.00,.00,.00,.00,.00,.00,.00,.00,1.00
+6.00,1.00,1.00,.00,.00,.00,.00,.00,.00,.00,.00,.00,.00,2.00
+7.00,.00,.00,.00,.00,.00,.00,.00,1.00,.00,.00,.00,.00,1.00
+8.00,.00,.00,.00,.00,.00,.00,.00,.00,1.00,.00,.00,.00,1.00
+9.00,.00,.00,.00,.00,.00,.00,.00,.00,.00,1.00,.00,.00,1.00
+10.00,.00,.00,.00,.00,.00,.00,.00,.00,.00,.00,1.00,1.00,2.00
+Total,1.00,1.00,1.00,1.00,1.00,1.00,1.00,1.00,1.00,1.00,1.00,1.00,12.00
+
+Table: Symmetric measures.
+Category,Statistic,Value,Asymp. Std. Error,Approx. T,Approx. Sig.
+Ordinal by Ordinal,Spearman Correlation,.66,.14,2.59,
+Interval by Interval,Pearson's R,.67,.13,2.93,
+N of Valid Cases,,12,,,
+])
+
+# Test 5.
+AT_DATA([pearson5.sps], [dnl
+* From http://www.statisticslectures.com/topics/pearsonr/.
+DATA LIST FREE/x y.
+BEGIN DATA.
+18 15000
+25 29000
+57 68000
+45 52000
+26 32000
+64 80000
+37 41000
+40 45000
+24 26000
+33 33000
+END DATA.
+CROSSTABS x BY y/STATISTICS=CORR.
+])
+AT_CHECK([pspp -O format=csv pearson5.sps], [0], [dnl
+Table: Summary.
+,Cases,,,,,
+,Valid,,Missing,,Total,
+,N,Percent,N,Percent,N,Percent
+x * y,10,100.0%,0,0.0%,10,100.0%
+
+Table: x * y [[count]].
+,y,,,,,,,,,,
+x,15000.00,26000.00,29000.00,32000.00,33000.00,41000.00,45000.00,52000.00,68000.00,80000.00,Total
+18.00,1.00,.00,.00,.00,.00,.00,.00,.00,.00,.00,1.00
+24.00,.00,1.00,.00,.00,.00,.00,.00,.00,.00,.00,1.00
+25.00,.00,.00,1.00,.00,.00,.00,.00,.00,.00,.00,1.00
+26.00,.00,.00,.00,1.00,.00,.00,.00,.00,.00,.00,1.00
+33.00,.00,.00,.00,.00,1.00,.00,.00,.00,.00,.00,1.00
+37.00,.00,.00,.00,.00,.00,1.00,.00,.00,.00,.00,1.00
+40.00,.00,.00,.00,.00,.00,.00,1.00,.00,.00,.00,1.00
+45.00,.00,.00,.00,.00,.00,.00,.00,1.00,.00,.00,1.00
+57.00,.00,.00,.00,.00,.00,.00,.00,.00,1.00,.00,1.00
+64.00,.00,.00,.00,.00,.00,.00,.00,.00,.00,1.00,1.00
+Total,1.00,1.00,1.00,1.00,1.00,1.00,1.00,1.00,1.00,1.00,10.00
+
+Table: Symmetric measures.
+Category,Statistic,Value,Asymp. Std. Error,Approx. T,Approx. Sig.
+Ordinal by Ordinal,Spearman Correlation,1.00,.00,3.76,
+Interval by Interval,Pearson's R,.99,.00,3.86,
+N of Valid Cases,,10,,,
+])
+AT_CLEANUP
Unnamed repository; edit this file to name it for gitweb.