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], [0], [dnl
  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.627,35,,,
  70 ])
  71
  72 AT_CLEANUP