tests/language/stats/glm.at

   1 AT_BANNER([GLM procedure])
   2
   3 AT_SETUP([GLM latin square design])
   4
   5 dnl This example comes from :
   6 dnl  http://ssnds.uwo.ca/statsexamples/spssanova/latinsquareresults.html
   7 AT_DATA([latin.sps], [dnl
   8 set format = F20.3.
   9 data list notable  fixed /a 1 b 3 c 5 y 7-10(2).
  10 begin data.
  11 1 1 6  3.5
  12 1 2 2  8.9
  13 1 3 3  9.6
  14 1 4 4 10.5
  15 1 5 5  3.1
  16 1 6 1  5.9
  17 2 1 2  4.2
  18 2 2 6  1.9
  19 2 3 5  3.7
  20 2 4 3 10.2
  21 2 5 1  7.2
  22 2 6 4  7.6
  23 3 1 1  6.7
  24 3 2 4  5.8
  25 3 3 6 -2.7
  26 3 4 2  4.6
  27 3 5 3  4.0
  28 3 6 5 -0.7
  29 4 1 4  6.6
  30 4 2 1  4.5
  31 4 3 2  3.7
  32 4 4 5  3.7
  33 4 5 6 -3.3
  34 4 6 3  3.0
  35 5 1 3  4.1
  36 5 2 5  2.4
  37 5 3 4  6.0
  38 5 4 1  5.1
  39 5 5 2  3.5
  40 5 6 6  4.0
  41 6 1 5  3.8
  42 6 2 3  5.8
  43 6 3 1  7.0
  44 6 4 6  3.8
  45 6 5 4  5.0
  46 6 6 2  8.6
  47 end data.
  48
  49 variable labels a 'Factor A' b 'Factor B' c 'Factor C' y 'Criterion'.
  50
  51 glm y by   b a c
  52   /intercept=include
  53   /criteria=alpha(.05)
  54   /design = a b c
  55   .
  56 ])
  57
  58 AT_CHECK([pspp -O format=csv latin.sps | sed 's/329.62[[678]]/329.62/'], [0],
  59   [dnl
  60 warning: GLM is experimental.  Do not rely on these results.
  61
  62 Table: Tests of Between-Subjects Effects
  63 Source,Type III Sum of Squares,df,Mean Square,F,Sig.
  64 Corrected Model,263.064,15,17.538,5.269,.000
  65 Intercept,815.103,1,815.103,244.910,.000
  66 Factor A,78.869,5,15.774,4.739,.005
  67 Factor B,28.599,5,5.720,1.719,.176
  68 Factor C,155.596,5,31.119,9.350,.000
  69 Error,66.563,20,3.328,,
  70 Total,1144.730,36,,,
  71 Corrected Total,329.62,35,,,
  72 ])
  73
  74 AT_CLEANUP
  75
  76 AT_SETUP([GLM 2 by 2 factorial design])
  77
  78 AT_DATA([2by2.sps], [dnl
  79 set format = F20.3.
  80 data list notable  list /Factor0 * Factor1 * points (F10).
  81 begin data.
  82 1       4       332
  83 1       4       380
  84 1       4       371
  85 1       4       366
  86 1       4       354
  87 1       0       259.5
  88 1       0       302.5
  89 1       0       296
  90 1       0       349
  91 1       0       309
  92 2       4       354.67
  93 2       4       353.5
  94 2       4       304
  95 2       4       365
  96 2       4       339
  97 2       0       306
  98 2       0       339
  99 2       0       353
 100 2       0       351
 101 2       0       333
 102 end data.
 103
 104 glm points by Factor0 Factor1
 105   /intercept=include
 106   /criteria=alpha(.05)
 107   .
 108 ])
 109
 110
 111 AT_CHECK([pspp -O format=csv 2by2.sps ], [0],
 112   [dnl
 113 warning: GLM is experimental.  Do not rely on these results.
 114
 115 Table: Tests of Between-Subjects Effects
 116 Source,Type III Sum of Squares,df,Mean Square,F,Sig.
 117 Corrected Model,8667.053,3,2889.018,5.043,.012
 118 Intercept,2256018.640,1,2256018.640,3937.693,.000
 119 Factor0,313.394,1,313.394,.547,.470
 120 Factor1,5157.508,1,5157.508,9.002,.008
 121 Factor0 * Factor1,3196.150,1,3196.150,5.579,.031
 122 Error,9166.865,16,572.929,,
 123 Total,2273852.559,20,,,
 124 Corrected Total,17833.918,19,,,
 125 ])
 126
 127 AT_CLEANUP
 128
 129
 130 AT_SETUP([GLM Type 1 Sums of Squares])
 131
 132 dnl  The following example comes from
 133 dnl  http://www.uvm.edu/~dhowell/StatPages/More_Stuff/Type1-3.pdf
 134 AT_DATA([type1.sps], [dnl
 135 set decimal = dot.
 136 set format=F20.3.
 137 data list notable list /dv * Agrp * B0 * B1 * B2 * i0 * i1 * i2 * sss *.
 138 begin data.
 139 5 1 1 0 0 1 0 0 1.00
 140 7 1 1 0 0 1 0 0 1.00
 141 9 1 1 0 0 1 0 0 1.00
 142 8 1 1 0 0 1 0 0 1.00
 143 2 1 0 1 0 0 1 0 1.00
 144 5 1 0 1 0 0 1 0 1.00
 145 7 1 0 1 0 0 1 0 1.00
 146 3 1 0 1 0 0 1 0 1.00
 147 9 1 0 1 0 0 1 0 1.00
 148 8 1 0 0 1 0 0 1 1.00
 149 11 1 0 0 1 0 0 1 1.00
 150 12 1 0 0 1 0 0 1 1.00
 151 14 1 0 0 1 0 0 1 1.00
 152 11 1 -1 -1 -1 -1 -1 -1 1.00
 153 15 1 -1 -1 -1 -1 -1 -1 1.00
 154 16 1 -1 -1 -1 -1 -1 -1 1.00
 155 10 1 -1 -1 -1 -1 -1 -1 1.00
 156 9 1 -1 -1 -1 -1 -1 -1 1.00
 157 7 -1 1 0 0 -1 0 0 2.00
 158 9 -1 1 0 0 -1 0 0 2.00
 159 10 -1 1 0 0 -1 0 0 2.00
 160 9 -1 1 0 0 -1 0 0 2.00
 161 3 -1 0 1 0 0 -1 0 2.00
 162 8 -1 0 1 0 0 -1 0 2.00
 163 9 -1 0 1 0 0 -1 0 2.00
 164 11 -1 0 1 0 0 -1 0 2.00
 165 9 -1 0 0 1 0 0 -1 2.00
 166 12 -1 0 0 1 0 0 -1 2.00
 167 14 -1 0 0 1 0 0 -1 2.00
 168 8 -1 0 0 1 0 0 -1 2.00
 169 7 -1 0 0 1 0 0 -1 2.00
 170 11 -1 -1 -1 -1 1 1 1 2.00
 171 14 -1 -1 -1 -1 1 1 1 2.00
 172 10 -1 -1 -1 -1 1 1 1 2.00
 173 12 -1 -1 -1 -1 1 1 1 2.00
 174 13 -1 -1 -1 -1 1 1 1 2.00
 175 11 -1 -1 -1 -1 1 1 1 2.00
 176 12 -1 -1 -1 -1 1 1 1 2.00
 177 end data.
 178
 179 do if B0 = -1 AND B1 = -1 AND B2 = -1.
 180 compute Bgrp = -1.
 181 end if.
 182
 183 do if B0 = 0 AND B1 = 0 AND B2 = 1.
 184 compute Bgrp = 1.
 185 end if.
 186
 187 do if B0 = 0 AND B1 = 1 AND B2 = 0.
 188 compute Bgrp = 2.
 189 end if.
 190
 191 do if B0 = 1 AND B1 = 0 AND B2 = 0.
 192 compute Bgrp = 3.
 193 end if.
 194
 195
 196 do if B0 = 0 AND B1 = 1 AND B2 = 0.
 197 compute Bgrp = 4.
 198 end if.
 199
 200
 201 glm dv by Agrp Bgrp
 202         /method = sstype (1)
 203         .
 204
 205 glm dv by Agrp Bgrp
 206         /method = sstype (1)
 207         /design Bgrp Agrp Bgrp * Agrp
 208         .
 209 ])
 210
 211
 212 AT_CHECK([pspp -O format=csv type1.sps], [0],
 213   [dnl
 214 warning: GLM is experimental.  Do not rely on these results.
 215
 216 Table: Tests of Between-Subjects Effects
 217 Source,Type I Sum of Squares,df,Mean Square,F,Sig.
 218 Corrected Model,216.017,7,30.860,5.046,.001
 219 Intercept,3410.526,1,3410.526,557.709,.000
 220 Agrp,9.579,1,9.579,1.566,.220
 221 Bgrp,186.225,3,62.075,10.151,.000
 222 Agrp * Bgrp,20.212,3,6.737,1.102,.364
 223 Error,183.457,30,6.115,,
 224 Total,3810.000,38,,,
 225 Corrected Total,399.474,37,,,
 226
 227 warning: GLM is experimental.  Do not rely on these results.
 228
 229 Table: Tests of Between-Subjects Effects
 230 Source,Type I Sum of Squares,df,Mean Square,F,Sig.
 231 Corrected Model,216.017,7,30.860,5.046,.001
 232 Intercept,3410.526,1,3410.526,557.709,.000
 233 Bgrp,193.251,3,64.417,10.534,.000
 234 Agrp,2.553,1,2.553,.418,.523
 235 Bgrp * Agrp,20.212,3,6.737,1.102,.364
 236 Error,183.457,30,6.115,,
 237 Total,3810.000,38,,,
 238 Corrected Total,399.474,37,,,
 239 ])
 240
 241 AT_CLEANUP