pintos-os.org Git - pspp/blob - 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 7  3
   8 4  3
   9 3  1
  10 2  1
  11 1  1
  12 4  2
  13 2  2
  14 3  2
  15 5  3
  16 1  1
  17 4  1
  18 5  2
  19 2  2
  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,Sig.
  48 Breaking Strain,.09,2,12,.91
  49
  50 Table: ANOVA
  51 ,,Sum of Squares,df,Mean Square,F,Sig.
  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,Sig.
 121 Breaking Strain,1.09,1,5,.35
 122
 123 Table: ANOVA
 124 ,,Sum of Squares,df,Mean Square,F,Sig.
 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,Sig.
 154 Breaking Strain,.92,1,6,.37
 155
 156 Table: ANOVA
 157 ,,Sum of Squares,df,Mean Square,F,Sig.
 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,Sig.
 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,Sig.
 365 QUALITY,.09,2,12,.91
 366
 367 Table: ANOVA
 368 ,,Sum of Squares,df,Mean Square,F,Sig.
 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 AT_CHECK([pspp -o pspp.csv multivar.sps])
 418 dnl Some machines return 3.88 instead of 3.87 below (see bug #31611).
 419 AT_CHECK([sed 's/^,Within Groups,3.88/,Within Groups,3.87/' pspp.csv], [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,Sig.
 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,Sig.
 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,.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,.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,.37
 476 ,,2,11.73,14.53,.81,9.88,.44
 477 ])
 478
 479 AT_CLEANUP
 480
 481
 482
 483 dnl Tests that everything treats weights properly
 484 AT_SETUP([ONEWAY vs. weights])
 485
 486 AT_DATA([oneway-unweighted.sps],
 487 [DATA LIST NOTABLE LIST /QUALITY * BRAND * W *.
 488 BEGIN DATA
 489 3  1   1
 490 2  1   1
 491 1  1   1
 492 1  1   1
 493 4  1   1
 494 5  2   1
 495 2  2   1
 496 4  2   1
 497 4  2   1
 498 4  2   1
 499 2  2   1
 500 2  2   1
 501 3  2   1
 502 7  3   1
 503 4  3   1
 504 5  3   1
 505 5  3   1
 506 3  3   1
 507 6  3   1
 508 END DATA.
 509
 510 WEIGHT BY W.
 511
 512 ONEWAY
 513         quality BY brand
 514         /STATISTICS descriptives homogeneity
 515         .
 516 ])
 517
 518 AT_CHECK([pspp -o pspp-unweighted.csv oneway-unweighted.sps], [0], [ignore], [ignore])
 519
 520 AT_DATA([oneway-weighted.sps],
 521 [DATA LIST NOTABLE LIST /QUALITY * BRAND * W *.
 522 BEGIN DATA
 523 3  1   1
 524 2  1   1
 525 1  1   2
 526 4  1   1
 527 5  2   1
 528 2  2   1
 529 4  2   3
 530 2  2   2
 531 3  2   1
 532 7  3   1
 533 4  3   1
 534 5  3   2
 535 3  3   1
 536 6  3   1
 537 END DATA.
 538
 539 WEIGHT BY W.
 540
 541 ONEWAY
 542         quality BY brand
 543         /STATISTICS descriptives homogeneity
 544         .
 545 ])
 546
 547 AT_CHECK([pspp -o pspp-weighted.csv oneway-weighted.sps], [0], [ignore], [ignore])
 548
 549 AT_CHECK([diff pspp-weighted.csv pspp-unweighted.csv], [0])
 550
 551 AT_CLEANUP
 552
 553
 554
 555 AT_SETUP([ONEWAY posthoc LSD and BONFERRONI])
 556 AT_DATA([oneway-pig.sps],[dnl
 557 SET FORMAT F12.3.
 558 data list notable list /pigmentation * family *.
 559 begin data.
 560 36 1
 561 39 1
 562 43 1
 563 38 1
 564 37 1
 565 46 2
 566 47 2
 567 47 2
 568 47 2
 569 43 2
 570 40 3
 571 50 3
 572 44 3
 573 48 3
 574 50 3
 575 45 4
 576 53 4
 577 56 4
 578 52 4
 579 56 4
 580 end data.
 581
 582
 583 oneway pigmentation by family
 584         /statistics = descriptives
 585         /posthoc = lsd bonferroni alpha (0.05)
 586          .
 587 ])
 588
 589 AT_CHECK([pspp -O format=csv oneway-pig.sps], [0],
 590 [Table: Descriptives
 591 ,,,,,,95% Confidence Interval for Mean,,,
 592 ,,N,Mean,Std. Deviation,Std. Error,Lower Bound,Upper Bound,Minimum,Maximum
 593 pigmentation,1.000,5,38.600,2.702,1.208,35.245,41.955,36.000,43.000
 594 ,2.000,5,46.000,1.732,.775,43.849,48.151,43.000,47.000
 595 ,3.000,5,46.400,4.336,1.939,41.016,51.784,40.000,50.000
 596 ,4.000,5,52.400,4.506,2.015,46.806,57.994,45.000,56.000
 597 ,Total,20,45.850,5.967,1.334,43.057,48.643,36.000,56.000
 598
 599 Table: ANOVA
 600 ,,Sum of Squares,df,Mean Square,F,Sig.
 601 pigmentation,Between Groups,478.950,3,159.650,12.927,.000
 602 ,Within Groups,197.600,16,12.350,,
 603 ,Total,676.550,19,,,
 604
 605 Table: Multiple Comparisons (pigmentation)
 606 ,,,Mean Difference,,,95% Confidence Interval,
 607 ,(I) family,(J) family,(I - J),Std. Error,Sig.,Lower Bound,Upper Bound
 608 LSD,1.000,2.000,-7.400,2.223,.004,-12.112,-2.688
 609 ,,3.000,-7.800,2.223,.003,-12.512,-3.088
 610 ,,4.000,-13.800,2.223,.000,-18.512,-9.088
 611 ,2.000,1.000,7.400,2.223,.004,2.688,12.112
 612 ,,3.000,-.400,2.223,.859,-5.112,4.312
 613 ,,4.000,-6.400,2.223,.011,-11.112,-1.688
 614 ,3.000,1.000,7.800,2.223,.003,3.088,12.512
 615 ,,2.000,.400,2.223,.859,-4.312,5.112
 616 ,,4.000,-6.000,2.223,.016,-10.712,-1.288
 617 ,4.000,1.000,13.800,2.223,.000,9.088,18.512
 618 ,,2.000,6.400,2.223,.011,1.688,11.112
 619 ,,3.000,6.000,2.223,.016,1.288,10.712
 620 Bonferroni,1.000,2.000,-7.400,2.223,.025,-14.086,-.714
 621 ,,3.000,-7.800,2.223,.017,-14.486,-1.114
 622 ,,4.000,-13.800,2.223,.000,-20.486,-7.114
 623 ,2.000,1.000,7.400,2.223,.025,.714,14.086
 624 ,,3.000,-.400,2.223,1.000,-7.086,6.286
 625 ,,4.000,-6.400,2.223,.065,-13.086,.286
 626 ,3.000,1.000,7.800,2.223,.017,1.114,14.486
 627 ,,2.000,.400,2.223,1.000,-6.286,7.086
 628 ,,4.000,-6.000,2.223,.095,-12.686,.686
 629 ,4.000,1.000,13.800,2.223,.000,7.114,20.486
 630 ,,2.000,6.400,2.223,.065,-.286,13.086
 631 ,,3.000,6.000,2.223,.095,-.686,12.686
 632 ])
 633
 634 AT_CLEANUP
 635
 636
 637 AT_SETUP([ONEWAY posthoc Tukey HSD and Games-Howell])
 638 AT_DATA([oneway-tukey.sps],[dnl
 639 set format = f11.3.
 640 data list notable list /libido * dose *.
 641 begin data.
 642 3 0
 643 2 0
 644 1 0
 645 1 0
 646 4 0
 647 5 1
 648 2 1
 649 4 1
 650 2 1
 651 3 1
 652 7 2
 653 4 2
 654 5 2
 655 3 2
 656 6 2
 657 end data.
 658
 659 variable label dose 'Dose of Viagra'.
 660
 661 add value labels dose 0 'Placebo' 1 '1 Dose' 2 '2 Doses'.
 662
 663 oneway libido by dose
 664         /posthoc tukey gh.
 665 ])
 666
 667 AT_CHECK([pspp -O format=csv oneway-tukey.sps], [0],
 668 [Table: ANOVA
 669 ,,Sum of Squares,df,Mean Square,F,Sig.
 670 libido,Between Groups,20.133,2,10.067,5.119,.025
 671 ,Within Groups,23.600,12,1.967,,
 672 ,Total,43.733,14,,,
 673
 674 Table: Multiple Comparisons (libido)
 675 ,,,Mean Difference,,,95% Confidence Interval,
 676 ,(I) Dose of Viagra,(J) Dose of Viagra,(I - J),Std. Error,Sig.,Lower Bound,Upper Bound
 677 Tukey HSD,Placebo,1 Dose,-1.000,.887,.516,-3.366,1.366
 678 ,,2 Doses,-2.800,.887,.021,-5.166,-.434
 679 ,1 Dose,Placebo,1.000,.887,.516,-1.366,3.366
 680 ,,2 Doses,-1.800,.887,.147,-4.166,.566
 681 ,2 Doses,Placebo,2.800,.887,.021,.434,5.166
 682 ,,1 Dose,1.800,.887,.147,-.566,4.166
 683 Games-Howell,Placebo,1 Dose,-1.000,.887,.479,-3.356,1.356
 684 ,,2 Doses,-2.800,.887,.039,-5.439,-.161
 685 ,1 Dose,Placebo,1.000,.887,.479,-1.356,3.356
 686 ,,2 Doses,-1.800,.887,.185,-4.439,.839
 687 ,2 Doses,Placebo,2.800,.887,.039,.161,5.439
 688 ,,1 Dose,1.800,.887,.185,-.839,4.439
 689 ])
 690
 691 AT_CLEANUP
 692
 693 AT_SETUP([ONEWAY posthoc Sidak])
 694 AT_DATA([oneway-sidak.sps],[dnl
 695 SET FORMAT F20.4.
 696
 697 DATA LIST notable LIST /program score.
 698 BEGIN DATA.
 699 1   9
 700 1  12
 701 1  14
 702 1  11
 703 1  13
 704 2  10
 705 2   6
 706 2   9
 707 2   9
 708 2  10
 709 3  12
 710 3  14
 711 3  11
 712 3  13
 713 3  11
 714 4   9
 715 4   8
 716 4  11
 717 4   7
 718 4   8
 719 END DATA.
 720
 721 ONEWAY
 722   score BY program
 723   /MISSING ANALYSIS
 724   /POSTHOC = SIDAK.
 725 ])
 726
 727 AT_CHECK([pspp -O format=csv oneway-sidak.sps], [0],
 728 [Table: ANOVA
 729 ,,Sum of Squares,df,Mean Square,F,Sig.
 730 score,Between Groups,54.9500,3,18.3167,7.0449,.0031
 731 ,Within Groups,41.6000,16,2.6000,,
 732 ,Total,96.5500,19,,,
 733
 734 Table: Multiple Comparisons (score)
 735 ,,,Mean Difference,,,95% Confidence Interval,
 736 ,(I) program,(J) program,(I - J),Std. Error,Sig.,Lower Bound,Upper Bound
 737 Šidák,1.0000,2.0000,3.0000,1.0198,.0561,-.0575,6.0575
 738 ,,3.0000,-.4000,1.0198,.9993,-3.4575,2.6575
 739 ,,4.0000,3.2000,1.0198,.0375,.1425,6.2575
 740 ,2.0000,1.0000,-3.0000,1.0198,.0561,-6.0575,.0575
 741 ,,3.0000,-3.4000,1.0198,.0250,-6.4575,-.3425
 742 ,,4.0000,.2000,1.0198,1.0000,-2.8575,3.2575
 743 ,3.0000,1.0000,.4000,1.0198,.9993,-2.6575,3.4575
 744 ,,2.0000,3.4000,1.0198,.0250,.3425,6.4575
 745 ,,4.0000,3.6000,1.0198,.0166,.5425,6.6575
 746 ,4.0000,1.0000,-3.2000,1.0198,.0375,-6.2575,-.1425
 747 ,,2.0000,-.2000,1.0198,1.0000,-3.2575,2.8575
 748 ,,3.0000,-3.6000,1.0198,.0166,-6.6575,-.5425
 749 ])
 750
 751 AT_CLEANUP
 752
 753 AT_SETUP([ONEWAY posthoc Scheffe])
 754 AT_DATA([oneway-scheffe.sps],[dnl
 755 set format = f11.3.
 756 data list notable list /usage * group *.
 757 begin data.
 758 21.00     1
 759 19.00     1
 760 18.00     1
 761 25.00     1
 762 14.00     1
 763 13.00     1
 764 24.00     1
 765 19.00     1
 766 20.00     1
 767 21.00     1
 768 15.00     2
 769 10.00     2
 770 13.00     2
 771 16.00     2
 772 14.00     2
 773 24.00     2
 774 16.00     2
 775 14.00     2
 776 18.00     2
 777 16.00     2
 778 10.00     3
 779  7.00     3
 780 13.00     3
 781 20.00     3
 782   .00     3
 783  8.00     3
 784  6.00     3
 785  1.00     3
 786 12.00     3
 787 14.00     3
 788 18.00     4
 789 15.00     4
 790  3.00     4
 791 27.00     4
 792  6.00     4
 793 14.00     4
 794 13.00     4
 795 11.00     4
 796  9.00     4
 797 18.00     4
 798 end data.
 799
 800 variable label usage 'Days of Use'.
 801
 802 add value labels group 0 'none' 1 'one' 2 'two' 3 'three' 4 'four'.
 803
 804 oneway usage by group
 805         /posthoc scheffe.
 806 ])
 807
 808 AT_CHECK([pspp -O format=csv oneway-scheffe.sps], [0],
 809 [Table: ANOVA
 810 ,,Sum of Squares,df,Mean Square,F,Sig.
 811 Days of Use,Between Groups,555.275,3,185.092,6.663,.001
 812 ,Within Groups,1000.100,36,27.781,,
 813 ,Total,1555.375,39,,,
 814
 815 Table: Multiple Comparisons (Days of Use)
 816 ,,,Mean Difference,,,95% Confidence Interval,
 817 ,(I) group,(J) group,(I - J),Std. Error,Sig.,Lower Bound,Upper Bound
 818 Scheffé,one,two,3.800,2.357,.467,-3.112,10.712
 819 ,,three,10.300,2.357,.001,3.388,17.212
 820 ,,four,6.000,2.357,.110,-.912,12.912
 821 ,two,one,-3.800,2.357,.467,-10.712,3.112
 822 ,,three,6.500,2.357,.072,-.412,13.412
 823 ,,four,2.200,2.357,.832,-4.712,9.112
 824 ,three,one,-10.300,2.357,.001,-17.212,-3.388
 825 ,,two,-6.500,2.357,.072,-13.412,.412
 826 ,,four,-4.300,2.357,.358,-11.212,2.612
 827 ,four,one,-6.000,2.357,.110,-12.912,.912
 828 ,,two,-2.200,2.357,.832,-9.112,4.712
 829 ,,three,4.300,2.357,.358,-2.612,11.212
 830 ])
 831
 832 AT_CLEANUP
 833
 834
 835 AT_SETUP([ONEWAY bad contrast count])
 836
 837 AT_DATA([oneway-bad-contrast.sps],[dnl
 838 DATA LIST NOTABLE LIST /height * weight * temperature * sex *.
 839 BEGIN DATA.
 840 1884     88.6       39.97     0
 841 1801     90.9       39.03     0
 842 1801     91.7       38.98     0
 843 1607     56.3       36.26     1
 844 1608     46.3       46.26     1
 845 1607     55.9       37.84     1
 846 1604     56.6       36.81     1
 847 1606     56.1       34.56     1
 848 END DATA.
 849
 850 ONEWAY /VARIABLES= height weight temperature BY sex
 851  /CONTRAST = -1  1
 852  /CONTRAST = -3  3
 853  /CONTRAST =  2 -2  1
 854  /CONTRAST = -9  9
 855  .
 856 ])
 857
 858
 859 AT_CHECK([pspp -O format=csv oneway-bad-contrast.sps], [0], [dnl
 860 "oneway-bad-contrast.sps:18: warning: ONEWAY: In contrast list 3, the number of coefficients (3) does not equal the number of groups (2). This contrast list will be ignored."
 861
 862 Table: ANOVA
 863 ,,Sum of Squares,df,Mean Square,F,Sig.
 864 height,Between Groups,92629.63,1,92629.63,120.77,.00
 865 ,Within Groups,4601.87,6,766.98,,
 866 ,Total,97231.50,7,,,
 867 weight,Between Groups,2451.65,1,2451.65,174.59,.00
 868 ,Within Groups,84.25,6,14.04,,
 869 ,Total,2535.90,7,,,
 870 temperature,Between Groups,1.80,1,1.80,.13,.73
 871 ,Within Groups,84.55,6,14.09,,
 872 ,Total,86.36,7,,,
 873
 874 Table: Contrast Coefficients
 875 ,,sex,
 876 ,,.00,1.00
 877 Contrast,1,-1,1
 878 ,2,-3,3
 879 ,3,-9,9
 880
 881 Table: Contrast Tests
 882 ,,Contrast,Value of Contrast,Std. Error,t,df,Sig. (2-tailed)
 883 height,Assume equal variances,1,-222.27,20.23,10.99,6,.00
 884 ,,2,-666.80,60.68,10.99,6,.00
 885 ,,3,-2000.40,182.03,10.99,6,.00
 886 ,Does not assume equal,1,-222.27,27.67,-8.03,2.00,.02
 887 ,,2,-666.80,83.02,-8.03,2.00,.02
 888 ,,3,-2000.40,249.07,-8.03,2.00,.02
 889 weight,Assume equal variances,1,-36.16,2.74,13.21,6,.00
 890 ,,2,-108.48,8.21,13.21,6,.00
 891 ,,3,-325.44,24.63,13.21,6,.00
 892 ,Does not assume equal,1,-36.16,2.19,-16.48,5.42,.00
 893 ,,2,-108.48,6.58,-16.48,5.42,.00
 894 ,,3,-325.44,19.75,-16.48,5.42,.00
 895 temperature,Assume equal variances,1,-.98,2.74,.36,6,.73
 896 ,,2,-2.94,8.22,.36,6,.73
 897 ,,3,-8.83,24.67,.36,6,.73
 898 ,Does not assume equal,1,-.98,2.07,-.47,4.19,.66
 899 ,,2,-2.94,6.22,-.47,4.19,.66
 900 ,,3,-8.83,18.66,-.47,4.19,.66
 901 ])
 902
 903 AT_CLEANUP
 904
 905
 906 AT_SETUP([ONEWAY crash on single category independent variable])
 907 AT_DATA([crash.sps],[
 908 input program.
 909 loop #i = 1 to 10.
 910 compute test = #i.
 911 end case.
 912 end loop.
 913 end file.
 914 end input program.
 915
 916 compute x = 1.
 917
 918 oneway test by x.
 919 ])
 920
 921 AT_CHECK([pspp -O format=csv crash.sps], [0], [ignore])
 922
 923 AT_CLEANUP
 924
 925
 926
 927 AT_SETUP([ONEWAY crash on missing dependent variable])
 928 AT_DATA([crash2.sps],[dnl
 929 data list notable list /dv1 * dv2  *  y * .
 930 begin data.
 931 2  .  2
 932 1  .  2
 933 1  .  1
 934 2  .  4
 935 3  .  4
 936 4  .  4
 937 5  .  4
 938 end data.
 939
 940 ONEWAY
 941        /VARIABLES= dv1 dv2  BY y
 942        /STATISTICS = DESCRIPTIVES
 943        /POSTHOC = BONFERRONI LSD SCHEFFE SIDAK TUKEY
 944        /MISSING = ANALYSIS
 945        .
 946 ])
 947
 948 AT_CHECK([pspp -O format=csv crash2.sps], [0], [ignore])
 949
 950 AT_CLEANUP
 951
 952
 953
 954
 955 AT_SETUP([ONEWAY Games-Howell test with few cases])
 956 AT_DATA([crash3.sps],[dnl
 957 data list notable list /dv * y * .
 958 begin data.
 959 2 2
 960 1 2
 961 1 1
 962 2 4
 963 3 4
 964 end data.
 965
 966 ONEWAY
 967  /VARIABLES= dv BY y
 968  /POSTHOC = GH
 969  .
 970 ])
 971
 972 AT_CHECK([pspp -O format=csv crash3.sps], [0], [ignore])
 973
 974 AT_CLEANUP