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   /method=sstype(3)
  53   /intercept=include
  54   /criteria=alpha(.05)
  55   /design = a b c
  56   .
  57 ])
  58
  59 AT_CHECK([pspp -O format=csv latin.sps | sed 's/329.62[[678]]/329.62/'], [0],
  60   [Table: Tests of Between-Subjects Effects
  61 Source,Type III Sum of Squares,df,Mean Square,F,Sig.
  62 Corrected Model,263.064,15,17.538,5.269,.000
  63 Intercept,815.103,1,815.103,244.910,.000
  64 Factor A,78.869,5,15.774,4.739,.005
  65 Factor B,28.599,5,5.720,1.719,.176
  66 Factor C,155.596,5,31.119,9.350,.000
  67 Error,66.563,20,3.328,,
  68 Total,1144.730,36,,,
  69 Corrected Total,329.62,35,,,
  70 ])
  71
  72 AT_CLEANUP
  73
  74 AT_SETUP([GLM 2 by 2 factorial design])
  75
  76 AT_DATA([2by2.sps], [dnl
  77 set format = F20.3.
  78 data list notable  list /Factor0 * Factor1 * points (F10).
  79 begin data.
  80 1       4       332
  81 1       4       380
  82 1       4       371
  83 1       4       366
  84 1       4       354
  85 1       0       259.5
  86 1       0       302.5
  87 1       0       296
  88 1       0       349
  89 1       0       309
  90 2       4       354.67
  91 2       4       353.5
  92 2       4       304
  93 2       4       365
  94 2       4       339
  95 2       0       306
  96 2       0       339
  97 2       0       353
  98 2       0       351
  99 2       0       333
 100 end data.
 101
 102 glm points by Factor0 Factor1
 103   /method=sstype(3)
 104   /intercept=include
 105   /criteria=alpha(.05)
 106   .
 107 ])
 108
 109
 110 AT_CHECK([pspp -O format=csv 2by2.sps ], [0],
 111   [Table: Tests of Between-Subjects Effects
 112 Source,Type III Sum of Squares,df,Mean Square,F,Sig.
 113 Corrected Model,8667.053,3,2889.018,5.043,.012
 114 Intercept,2256018.640,1,2256018.640,3937.693,.000
 115 Factor0,313.394,1,313.394,.547,.470
 116 Factor1,5157.508,1,5157.508,9.002,.008
 117 Factor0 * Factor1,3196.150,1,3196.150,5.579,.031
 118 Error,9166.865,16,572.929,,
 119 Total,2273852.559,20,,,
 120 Corrected Total,17833.918,19,,,
 121 ])
 122
 123 AT_CLEANUP