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