1 dnl PSPP - a program for statistical analysis.
2 dnl Copyright (C) 2017 Free Software Foundation, Inc.
4 dnl This program is free software: you can redistribute it and/or modify
5 dnl it under the terms of the GNU General Public License as published by
6 dnl the Free Software Foundation, either version 3 of the License, or
7 dnl (at your option) any later version.
9 dnl This program is distributed in the hope that it will be useful,
10 dnl but WITHOUT ANY WARRANTY; without even the implied warranty of
11 dnl MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12 dnl GNU General Public License for more details.
14 dnl You should have received a copy of the GNU General Public License
15 dnl along with this program. If not, see <http://www.gnu.org/licenses/>.
17 AT_BANNER([GLM procedure])
19 AT_SETUP([GLM latin square design])
21 dnl This example comes from :
22 dnl http://ssnds.uwo.ca/statsexamples/spssanova/latinsquareresults.html
23 AT_DATA([latin.sps], [dnl
25 data list notable fixed /a 1 b 3 c 5 y 7-10(2).
65 variable labels a 'Factor A' b 'Factor B' c 'Factor C' y 'Criterion'.
74 AT_CHECK([pspp -O format=csv latin.sps | sed 's/329.62[[678]]/329.62/'], [0],
76 Table: Tests of Between-Subjects Effects
77 Source,Type III Sum of Squares,df,Mean Square,F,Sig.
78 Corrected Model,263.064,15,17.538,5.269,.000
79 Intercept,815.103,1,815.103,244.910,.000
80 Factor A,78.869,5,15.774,4.739,.005
81 Factor B,28.599,5,5.720,1.719,.176
82 Factor C,155.596,5,31.119,9.350,.000
83 Error,66.563,20,3.328,,
85 Corrected Total,329.62,35,,,
90 AT_SETUP([GLM 2 by 2 factorial design])
92 AT_DATA([2by2.sps], [dnl
94 data list notable list /Factor0 * Factor1 * points (F10).
118 glm points by Factor0 Factor1
125 AT_CHECK([pspp -O format=csv 2by2.sps ], [0],
127 Table: Tests of Between-Subjects Effects
128 Source,Type III Sum of Squares,df,Mean Square,F,Sig.
129 Corrected Model,8667.053,3,2889.018,5.043,.012
130 Intercept,2256018.640,1,2256018.640,3937.693,.000
131 Factor0,313.394,1,313.394,.547,.470
132 Factor1,5157.508,1,5157.508,9.002,.008
133 Factor0 * Factor1,3196.150,1,3196.150,5.579,.031
134 Error,9166.865,16,572.929,,
135 Total,2273852.559,20,,,
136 Corrected Total,17833.918,19,,,
142 AT_SETUP([GLM Type I and II Sums of Squares])
144 dnl The following example comes from
145 dnl http://www.uvm.edu/~dhowell/StatPages/More_Stuff/Type1-3.pdf
146 AT_DATA([data-inc.sps], [dnl
149 data list notable list /dv * Agrp * B0 * B1 * B2 * i0 * i1 * i2 * sss *.
161 11 1 0 0 1 0 0 1 1.00
162 12 1 0 0 1 0 0 1 1.00
163 14 1 0 0 1 0 0 1 1.00
164 11 1 -1 -1 -1 -1 -1 -1 1.00
165 15 1 -1 -1 -1 -1 -1 -1 1.00
166 16 1 -1 -1 -1 -1 -1 -1 1.00
167 10 1 -1 -1 -1 -1 -1 -1 1.00
168 9 1 -1 -1 -1 -1 -1 -1 1.00
169 7 -1 1 0 0 -1 0 0 2.00
170 9 -1 1 0 0 -1 0 0 2.00
171 10 -1 1 0 0 -1 0 0 2.00
172 9 -1 1 0 0 -1 0 0 2.00
173 3 -1 0 1 0 0 -1 0 2.00
174 8 -1 0 1 0 0 -1 0 2.00
175 9 -1 0 1 0 0 -1 0 2.00
176 11 -1 0 1 0 0 -1 0 2.00
177 9 -1 0 0 1 0 0 -1 2.00
178 12 -1 0 0 1 0 0 -1 2.00
179 14 -1 0 0 1 0 0 -1 2.00
180 8 -1 0 0 1 0 0 -1 2.00
181 7 -1 0 0 1 0 0 -1 2.00
182 11 -1 -1 -1 -1 1 1 1 2.00
183 14 -1 -1 -1 -1 1 1 1 2.00
184 10 -1 -1 -1 -1 1 1 1 2.00
185 12 -1 -1 -1 -1 1 1 1 2.00
186 13 -1 -1 -1 -1 1 1 1 2.00
187 11 -1 -1 -1 -1 1 1 1 2.00
188 12 -1 -1 -1 -1 1 1 1 2.00
191 do if B0 = -1 AND B1 = -1 AND B2 = -1.
195 do if B0 = 0 AND B1 = 0 AND B2 = 1.
199 do if B0 = 0 AND B1 = 1 AND B2 = 0.
203 do if B0 = 1 AND B1 = 0 AND B2 = 0.
208 AT_DATA([type1.sps], [dnl
209 include 'data-inc.sps'.
217 /design Bgrp Agrp Bgrp * Agrp
222 AT_CHECK([pspp -O format=csv type1.sps], [0],
224 Table: Tests of Between-Subjects Effects
225 Source,Type I Sum of Squares,df,Mean Square,F,Sig.
226 Corrected Model,216.017,7,30.860,5.046,.001
228 Agrp,9.579,1,9.579,1.566,.220
229 Bgrp,186.225,3,62.075,10.151,.000
230 Agrp * Bgrp,20.212,3,6.737,1.102,.364
231 Error,183.457,30,6.115,,
233 Corrected Total,399.474,37,,,
235 Table: Tests of Between-Subjects Effects
236 Source,Type I Sum of Squares,df,Mean Square,F,Sig.
237 Corrected Model,216.017,7,30.860,5.046,.001
239 Bgrp,193.251,3,64.417,10.534,.000
240 Agrp,2.553,1,2.553,.418,.523
241 Bgrp * Agrp,20.212,3,6.737,1.102,.364
242 Error,183.457,30,6.115,,
244 Corrected Total,399.474,37,,,
248 AT_DATA([type2.sps], [dnl
249 include 'data-inc.sps'.
257 AT_CHECK([pspp -O format=csv type2.sps], [0],
259 Table: Tests of Between-Subjects Effects
260 Source,Type II Sum of Squares,df,Mean Square,F,Sig.
261 Corrected Model,216.017,7,30.860,5.046,.001
263 Agrp,2.553,1,2.553,.418,.523
264 Bgrp,186.225,3,62.075,10.151,.000
265 Agrp * Bgrp,20.212,3,6.737,1.102,.364
266 Error,183.457,30,6.115,,
268 Corrected Total,399.474,37,,,
275 AT_SETUP([GLM excluded intercept])
277 dnl The following example comes from
279 dnl Rudolf N. Cardinal
280 dnl Graduate-level statistics for psychology and neuroscience
281 dnl ANOVA in practice, and complex ANOVA designs
282 dnl Version of 2 May 2004
284 dnl Downloaded from: http://egret.psychol.cam.ac.uk/psychology/graduate/Guide_to_ANOVA.pdf
286 AT_DATA([intercept-exclude.sps], [dnl
289 data list notable list /depvar * A *.
324 AT_CHECK([pspp -O format=csv intercept-exclude.sps], [0],
326 Table: Tests of Between-Subjects Effects
327 Source,Type III Sum of Squares,df,Mean Square,F,Sig.
328 Model,1636.826,2,818.413,43.556,.000
329 A,1636.826,2,818.413,43.556,.000
330 Error,338.216,18,18.790,,
333 Table: Tests of Between-Subjects Effects
334 Source,Type III Sum of Squares,df,Mean Square,F,Sig.
335 Corrected Model,98.568,1,98.568,5.246,.034
336 Intercept,1538.258,1,1538.258,81.867,.000
337 A,98.568,1,98.568,5.246,.034
338 Error,338.216,18,18.790,,
340 Corrected Total,436.784,19,,,
346 AT_SETUP([GLM missing values])
348 AT_DATA([glm.data], [dnl
387 AT_DATA([glm-miss.sps], [dnl
389 data list file='glm.data' notable fixed /a 1 b 3 c 5 y 7-10(2).
392 recode y (else=SYSMIS).
401 AT_CHECK([pspp -O format=csv glm-miss.sps], [0], [dnl
402 Table: Tests of Between-Subjects Effects
403 Source,Type III Sum of Squares,df,Mean Square,F,Sig.
404 Corrected Model,251.621,14,17.973,4.969,.002
405 Intercept,628.376,1,628.376,173.737,.000
406 a,72.929,4,18.232,5.041,.009
407 b,20.703,5,4.141,1.145,.380
408 c,135.179,5,27.036,7.475,.001
409 Error,54.253,15,3.617,,
411 Corrected Total,305.874,29,,,
416 AT_DATA([glm-miss2.sps], [dnl
418 data list file='glm.data' notable fixed /a 1 b 3 c 5 y 7-10(2).
428 AT_CHECK([pspp -O format=csv glm-miss2.sps], [0], [dnl
429 Table: Tests of Between-Subjects Effects
430 Source,Type III Sum of Squares,df,Mean Square,F,Sig.
431 Corrected Model,251.621,14,17.973,4.969,.002
432 Intercept,628.376,1,628.376,173.737,.000
433 a,72.929,4,18.232,5.041,.009
434 b,20.703,5,4.141,1.145,.380
435 c,135.179,5,27.036,7.475,.001
436 Error,54.253,15,3.617,,
438 Corrected Total,305.874,29,,,
442 dnl Now for some missing values in the factor variables.
444 AT_DATA([glm-miss3.sps], [dnl
446 data list file=glm.data notable fixed /a 1 b 3 c 5 y 7-10(2).
449 recode a (else=SYSMIS).
458 AT_CHECK([pspp -O format=csv glm-miss3.sps], [0], [dnl
459 Table: Tests of Between-Subjects Effects
460 Source,Type III Sum of Squares,df,Mean Square,F,Sig.
461 Corrected Model,251.621,14,17.973,4.969,.002
462 Intercept,628.376,1,628.376,173.737,.000
463 a,72.929,4,18.232,5.041,.009
464 b,20.703,5,4.141,1.145,.380
465 c,135.179,5,27.036,7.475,.001
466 Error,54.253,15,3.617,,
468 Corrected Total,305.874,29,,,