Oneway: Use covariance matrix and sweep operator
[pspp] / tests / language / stats / oneway.at
1 AT_BANNER([ONEWAY procedure])
2
3 AT_SETUP([ONEWAY basic operation])
4 AT_DATA([oneway.sps],
5   [DATA LIST NOTABLE LIST /QUALITY * BRAND * .
6 BEGIN DATA
7 3 1
8 2 1
9 1 1
10 1 1
11 4 1
12 5 2
13 2 2
14 4 2
15 2 2
16 3 2
17 7  3
18 4  3
19 5  3
20 3  3
21 6  3
22 END DATA
23
24 VARIABLE LABELS brand 'Manufacturer'.
25 VARIABLE LABELS quality 'Breaking Strain'.
26
27 VALUE LABELS /brand 1 'Aspeger' 2 'Bloggs' 3 'Charlies'.
28
29 ONEWAY
30         quality BY brand
31         /STATISTICS descriptives homogeneity
32         /CONTRAST =  -2 1 1 
33         /CONTRAST = 0 -1 1
34         .
35 ])
36
37 AT_CHECK([pspp -O format=csv oneway.sps], [0], 
38 [Table: Descriptives
39 ,,,,,,95% Confidence Interval for Mean,,,
40 ,,N,Mean,Std. Deviation,Std. Error,Lower Bound,Upper Bound,Minimum,Maximum
41 Breaking Strain,Aspeger,5,2.20,1.30,.58,.58,3.82,1.00,4.00
42 ,Bloggs,5,3.20,1.30,.58,1.58,4.82,2.00,5.00
43 ,Charlies,5,5.00,1.58,.71,3.04,6.96,3.00,7.00
44 ,Total,15,3.47,1.77,.46,2.49,4.45,1.00,7.00
45
46 Table: Test of Homogeneity of Variances
47 ,Levene Statistic,df1,df2,Significance
48 Breaking Strain,.09,2,12,.91
49
50 Table: ANOVA
51 ,,Sum of Squares,df,Mean Square,F,Significance
52 Breaking Strain,Between Groups,20.13,2,10.07,5.12,.02
53 ,Within Groups,23.60,12,1.97,,
54 ,Total,43.73,14,,,
55
56 Table: Contrast Coefficients
57 ,,Manufacturer,,
58 ,,Aspeger,Bloggs,Charlies
59 Contrast,1,-2,1,1
60 ,2,0,-1,1
61
62 Table: Contrast Tests
63 ,,Contrast,Value of Contrast,Std. Error,t,df,Sig. (2-tailed)
64 Breaking Strain,Assume equal variances,1,3.80,1.54,2.47,12,.03
65 ,,2,1.80,.89,2.03,12,.07
66 ,Does not assume equal,1,3.80,1.48,2.56,8.74,.03
67 ,,2,1.80,.92,1.96,7.72,.09
68 ])
69 AT_CLEANUP
70
71
72 AT_SETUP([ONEWAY with splits])
73 AT_DATA([oneway-splits.sps],
74 [DATA LIST NOTABLE LIST /QUALITY * BRAND * S *.
75 BEGIN DATA
76 3 1 1
77 2 1 1
78 1 1 1
79 1 1 1
80 4 1 1
81 5 2 1
82 2 2 1
83 4 2 2
84 2 2 2
85 3 2 2
86 7  3 2
87 4  3 2
88 5  3 2
89 3  3 2
90 6  3 2
91 END DATA
92
93 VARIABLE LABELS brand 'Manufacturer'.
94 VARIABLE LABELS quality 'Breaking Strain'.
95
96 VALUE LABELS /brand 1 'Aspeger' 2 'Bloggs' 3 'Charlies'.
97
98 SPLIT FILE by s.
99
100 ONEWAY
101         quality BY brand
102         /STATISTICS descriptives homogeneity
103         /CONTRAST =  -2 2
104         /CONTRAST = -1 1
105         .
106 ])
107
108 AT_CHECK([pspp -O format=csv oneway-splits.sps], [0], 
109 [Variable,Value,Label
110 S,1.00,
111
112 Table: Descriptives
113 ,,,,,,95% Confidence Interval for Mean,,,
114 ,,N,Mean,Std. Deviation,Std. Error,Lower Bound,Upper Bound,Minimum,Maximum
115 Breaking Strain,Aspeger,5,2.20,1.30,.58,.58,3.82,1.00,4.00
116 ,Bloggs,2,3.50,2.12,1.50,-15.56,22.56,2.00,5.00
117 ,Total,7,2.57,1.51,.57,1.17,3.97,1.00,5.00
118
119 Table: Test of Homogeneity of Variances
120 ,Levene Statistic,df1,df2,Significance
121 Breaking Strain,1.09,1,5,.35
122
123 Table: ANOVA
124 ,,Sum of Squares,df,Mean Square,F,Significance
125 Breaking Strain,Between Groups,2.41,1,2.41,1.07,.35
126 ,Within Groups,11.30,5,2.26,,
127 ,Total,13.71,6,,,
128
129 Table: Contrast Coefficients
130 ,,Manufacturer,
131 ,,Aspeger,Bloggs
132 Contrast,1,-2,2
133 ,2,-1,1
134
135 Table: Contrast Tests
136 ,,Contrast,Value of Contrast,Std. Error,t,df,Sig. (2-tailed)
137 Breaking Strain,Assume equal variances,1,2.60,2.52,1.03,5,.35
138 ,,2,1.30,1.26,1.03,5,.35
139 ,Does not assume equal,1,2.60,3.22,.81,1.32,.54
140 ,,2,1.30,1.61,.81,1.32,.54
141
142 Variable,Value,Label
143 S,2.00,
144
145 Table: Descriptives
146 ,,,,,,95% Confidence Interval for Mean,,,
147 ,,N,Mean,Std. Deviation,Std. Error,Lower Bound,Upper Bound,Minimum,Maximum
148 Breaking Strain,Bloggs,3,3.00,1.00,.58,.52,5.48,2.00,4.00
149 ,Charlies,5,5.00,1.58,.71,3.04,6.96,3.00,7.00
150 ,Total,8,4.25,1.67,.59,2.85,5.65,2.00,7.00
151
152 Table: Test of Homogeneity of Variances
153 ,Levene Statistic,df1,df2,Significance
154 Breaking Strain,.92,1,6,.37
155
156 Table: ANOVA
157 ,,Sum of Squares,df,Mean Square,F,Significance
158 Breaking Strain,Between Groups,7.50,1,7.50,3.75,.10
159 ,Within Groups,12.00,6,2.00,,
160 ,Total,19.50,7,,,
161
162 Table: Contrast Coefficients
163 ,,Manufacturer,
164 ,,Bloggs,Charlies
165 Contrast,1,-2,2
166 ,2,-1,1
167
168 Table: Contrast Tests
169 ,,Contrast,Value of Contrast,Std. Error,t,df,Sig. (2-tailed)
170 Breaking Strain,Assume equal variances,1,4.00,2.07,1.94,6,.10
171 ,,2,2.00,1.03,1.94,6,.10
172 ,Does not assume equal,1,4.00,1.83,2.19,5.88,.07
173 ,,2,2.00,.91,2.19,5.88,.07
174 ])
175
176 AT_CLEANUP
177
178
179 AT_SETUP([ONEWAY with missing values])
180 dnl Check that missing are treated properly
181 AT_DATA([oneway-missing1.sps],
182 [DATA LIST NOTABLE LIST /v1 * v2 * dep * vn *.
183 BEGIN DATA
184 . .  1  4
185 3 3  1  2
186 2 2  1  2
187 1 1  1  2
188 1 1  1  2
189 4 4  1  2
190 5 5  2  2
191 2 2  2  2
192 4 4  2  2
193 2 2  2  2
194 3 3  2  2
195 7 7  3  2
196 4 4  3  2
197 5 5  3  2
198 3 3  3  2
199 6 6  3  2
200 END DATA
201
202 ONEWAY
203         v1 v2 BY dep
204         /STATISTICS descriptives homogeneity
205         /MISSING ANALYSIS 
206         .
207 ])
208
209 AT_DATA([oneway-missing2.sps],
210 [DATA LIST NOTABLE LIST /v1 * v2 * dep * vn * .
211 BEGIN DATA
212 4 .  1  2 
213 3 3  1  2
214 2 2  1  2
215 1 1  1  2
216 1 1  1  2
217 4 4  1  2
218 5 5  2  2
219 2 2  2  2
220 4 4  2  2
221 2 2  2  2
222 3 3  2  2
223 7 7  3  2
224 4 4  3  2
225 5 5  3  2
226 3 3  3  2
227 6 6  3  2
228 END DATA
229
230 ONEWAY
231         v1 v2 BY dep
232         /STATISTICS descriptives homogeneity
233         /MISSING LISTWISE
234         .
235 ])
236
237
238
239 AT_CHECK([pspp -O format=csv oneway-missing1.sps > first.out], [0])
240
241 AT_CHECK([pspp -O format=csv oneway-missing2.sps > second.out], [0])
242
243 AT_CHECK([diff first.out second.out], [0], [])
244
245 dnl Now a test with missing values in the independent variable
246 AT_DATA([oneway-missing3.sps],
247 [DATA LIST NOTABLE LIST /v1 * v2 * dep * vn * .
248 BEGIN DATA
249 4 2  .  2 
250 3 3  1  2
251 2 2  1  2
252 1 1  1  2
253 1 1  1  2
254 4 4  1  2
255 5 5  2  2
256 2 2  2  2
257 4 4  2  2
258 2 2  2  2
259 3 3  2  2
260 7 7  3  2
261 4 4  3  2
262 5 5  3  4
263 3 3  3  2
264 6 6  3  2
265 END DATA
266
267 ONEWAY
268         v1 v2 BY dep
269         /STATISTICS descriptives homogeneity
270         /MISSING ANALYSIS
271         .
272 ])
273
274 AT_CHECK([pspp -O format=csv oneway-missing3.sps > third.out], [0])
275
276 AT_CHECK([diff first.out third.out], [0], [])
277
278 AT_CLEANUP
279
280
281
282
283
284 AT_SETUP([ONEWAY descriptives subcommand])
285
286 AT_DATA([oneway-descriptives.sps],
287   [DATA LIST NOTABLE LIST /QUALITY * BRAND * .
288 BEGIN DATA
289 13 11
290 12 11
291 11 11
292 11 11
293 14 11
294 15 25
295 12 25
296 14 25
297 12 25
298 13 25
299 17  301
300 14  301
301 15  301
302 13  301
303 16  301
304 END DATA
305
306
307 ONEWAY
308         quality BY brand
309         /STATISTICS descriptives
310         .
311 ])
312
313 AT_CHECK([pspp -O format=csv oneway-descriptives.sps], [0], 
314 [Table: Descriptives
315 ,,,,,,95% Confidence Interval for Mean,,,
316 ,,N,Mean,Std. Deviation,Std. Error,Lower Bound,Upper Bound,Minimum,Maximum
317 QUALITY,11.00,5,12.20,1.30,.58,10.58,13.82,11.00,14.00
318 ,25.00,5,13.20,1.30,.58,11.58,14.82,12.00,15.00
319 ,301.00,5,15.00,1.58,.71,13.04,16.96,13.00,17.00
320 ,Total,15,13.47,1.77,.46,12.49,14.45,11.00,17.00
321
322 Table: ANOVA
323 ,,Sum of Squares,df,Mean Square,F,Significance
324 QUALITY,Between Groups,20.13,2,10.07,5.12,.02
325 ,Within Groups,23.60,12,1.97,,
326 ,Total,43.73,14,,,
327 ])
328
329 AT_CLEANUP
330
331
332
333 AT_SETUP([ONEWAY homogeneity subcommand])
334
335 AT_DATA([oneway-homogeneity.sps],
336   [DATA LIST NOTABLE LIST /QUALITY * BRAND * .
337 BEGIN DATA
338 13 11
339 12 11
340 11 11
341 11 11
342 14 11
343 15 25
344 12 25
345 14 25
346 12 25
347 13 25
348 17  301
349 14  301
350 15  301
351 13  301
352 16  301
353 END DATA
354
355
356 ONEWAY
357         quality BY brand
358         /STATISTICS homogeneity
359         .
360 ])
361
362 AT_CHECK([pspp -O format=csv oneway-homogeneity.sps], [0], 
363 [Table: Test of Homogeneity of Variances
364 ,Levene Statistic,df1,df2,Significance
365 QUALITY,.09,2,12,.91
366
367 Table: ANOVA
368 ,,Sum of Squares,df,Mean Square,F,Significance
369 QUALITY,Between Groups,20.13,2,10.07,5.12,.02
370 ,Within Groups,23.60,12,1.97,,
371 ,Total,43.73,14,,,
372 ])
373
374 AT_CLEANUP
375
376
377
378 AT_SETUP([ONEWAY multiple variables])
379 dnl check that everything works ok when several different dependent variables are specified.
380 dnl This of course does not mean that we're doing a multivariate analysis.  It's just like
381 dnl running several tests at once.
382 AT_DATA([multivar.sps],
383 [DATA LIST notable LIST /x * y * z * g *.
384 begin data.
385 1 1 0 10
386 1 1 9 10
387 9 1 2 10
388 1 1 3 20
389 1 1 8 20
390 1 1 1 20
391 1 1 2 20
392 0 1 3 20
393 1 1 4 30
394 0 1 5 30
395 1 1 6 30
396 0 1 7 30
397 1 2 8 30
398 2 2 9 30
399 1 2 1 30
400 1 2 0 30
401 1 2 2 40
402 8 2 3 40
403 1 2 4 40
404 1 2 9 40
405 9 2 8 40
406 7 3 7 40
407 2 3 6 40
408 3 3 5 40
409 end data.
410
411 ONEWAY x y z by g
412         /STATISTICS = DESCRIPTIVES HOMOGENEITY
413         /CONTRAST 3  2 0 -5
414         /CONTRAST 2 -9 7  0
415         .
416 ])
417
418
419 AT_CHECK([pspp -O format=csv multivar.sps], [0], 
420 [Table: Descriptives
421 ,,,,,,95% Confidence Interval for Mean,,,
422 ,,N,Mean,Std. Deviation,Std. Error,Lower Bound,Upper Bound,Minimum,Maximum
423 x,10.00,3,3.67,4.62,2.67,-7.81,15.14,1.00,9.00
424 ,20.00,5,.80,.45,.20,.24,1.36,.00,1.00
425 ,30.00,8,.88,.64,.23,.34,1.41,.00,2.00
426 ,40.00,8,4.00,3.42,1.21,1.14,6.86,1.00,9.00
427 ,Total,24,2.25,2.83,.58,1.05,3.45,.00,9.00
428 y,10.00,3,1.00,.00,.00,1.00,1.00,1.00,1.00
429 ,20.00,5,1.00,.00,.00,1.00,1.00,1.00,1.00
430 ,30.00,8,1.50,.53,.19,1.05,1.95,1.00,2.00
431 ,40.00,8,2.38,.52,.18,1.94,2.81,2.00,3.00
432 ,Total,24,1.63,.71,.15,1.32,1.93,1.00,3.00
433 z,10.00,3,3.67,4.73,2.73,-8.07,15.41,.00,9.00
434 ,20.00,5,3.40,2.70,1.21,.05,6.75,1.00,8.00
435 ,30.00,8,5.00,3.21,1.13,2.32,7.68,.00,9.00
436 ,40.00,8,5.50,2.45,.87,3.45,7.55,2.00,9.00
437 ,Total,24,4.67,2.99,.61,3.40,5.93,.00,9.00
438
439 Table: Test of Homogeneity of Variances
440 ,Levene Statistic,df1,df2,Significance
441 x,18.76,3,20,.00
442 y,71.41,3,20,.00
443 z,.89,3,20,.46
444
445 Table: ANOVA
446 ,,Sum of Squares,df,Mean Square,F,Significance
447 x,Between Groups,56.16,3,18.72,2.92,.06
448 ,Within Groups,128.34,20,6.42,,
449 ,Total,184.50,23,,,
450 y,Between Groups,7.75,3,2.58,13.33,.00
451 ,Within Groups,3.87,20,.19,,
452 ,Total,11.63,23,,,
453 z,Between Groups,17.47,3,5.82,.62,.61
454 ,Within Groups,187.87,20,9.39,,
455 ,Total,205.33,23,,,
456
457 Table: Contrast Coefficients
458 ,,g,,,
459 ,,10.00,20.00,30.00,40.00
460 Contrast,1,3,2,0,-5
461 ,2,2,-9,7,0
462
463 Table: Contrast Tests
464 ,,Contrast,Value of Contrast,Std. Error,t,df,Sig. (2-tailed)
465 x,Assume equal variances,1,-7.40,6.67,1.11,20,.28
466 ,,2,6.26,12.32,.51,20,.62
467 ,Does not assume equal,1,-7.40,10.04,-.74,4.53,1.50
468 ,,2,6.26,5.85,1.07,2.87,.37
469 y,Assume equal variances,1,-6.88,1.16,5.94,20,.00
470 ,,2,3.50,2.14,1.63,20,.12
471 ,Does not assume equal,1,-6.88,.91,-7.51,7.00,2.00
472 ,,2,3.50,1.32,2.65,7.00,.03
473 z,Assume equal variances,1,-9.70,8.07,1.20,20,.24
474 ,,2,11.73,14.91,.79,20,.44
475 ,Does not assume equal,1,-9.70,9.57,-1.01,3.64,1.63
476 ,,2,11.73,14.53,.81,9.88,.44
477 ])
478
479 AT_CLEANUP