tests/language/stats/npar.at

   1 AT_BANNER([NPAR TESTS])
   2
   3 AT_SETUP([NPAR TESTS BINOMIAL, P < 0.5; N1/N2 < 1])
   4 AT_DATA([npar.sps], [dnl
   5 SET FORMAT F8.3.
   6
   7 DATA LIST LIST NOTABLE /x * w *.
   8 BEGIN DATA.
   9 1   6
  10 2   15
  11 END DATA.
  12
  13 WEIGHT BY w.
  14
  15 NPAR TESTS
  16         /BINOMIAL(0.3) = x
  17         .
  18 ])
  19 AT_CHECK([pspp -O format=csv npar.sps], [0], [dnl
  20 Table: Binomial Test
  21 ,,Category,N,Observed Prop.,Test Prop.,Exact Sig. (1-tailed)
  22 x,Group1,1.000,6.000,.286,.300,.551
  23 ,Group2,2.000,15.000,.714,,
  24 ,Total,,21.000,1.000,,
  25 ])
  26 AT_CLEANUP
  27
  28 AT_SETUP([NPAR TESTS BINOMIAL, P < 0.5; N1/N2 > 1])
  29 AT_DATA([npar.sps], [dnl
  30 SET FORMAT F8.3.
  31
  32 DATA LIST LIST NOTABLE /x (F8.0) w (F8.0).
  33 BEGIN DATA.
  34 1   7
  35 2   6
  36 END DATA.
  37
  38 WEIGHT BY w.
  39
  40 NPAR TESTS
  41         /BINOMIAL(0.4) = x
  42         .
  43 ])
  44 AT_CHECK([pspp -O format=csv npar.sps], [0], [dnl
  45 Table: Binomial Test
  46 ,,Category,N,Observed Prop.,Test Prop.,Exact Sig. (1-tailed)
  47 x,Group1,1,7,.538,.400,.229
  48 ,Group2,2,6,.462,,
  49 ,Total,,13,1.000,,
  50 ])
  51 AT_CLEANUP
  52
  53 AT_SETUP([NPAR TESTS BINOMIAL, P < 0.5; N1/N2 = 1])
  54 AT_DATA([npar.sps], [dnl
  55 SET FORMAT F8.3.
  56
  57 DATA LIST LIST NOTABLE /x (F8.0) w (F8.0).
  58 BEGIN DATA.
  59 1   8
  60 2   8
  61 END DATA.
  62
  63 WEIGHT BY w.
  64
  65 NPAR TESTS
  66         /BINOMIAL(0.4) = x
  67         .
  68 ])
  69 AT_CHECK([pspp -O format=csv npar.sps], [0], [dnl
  70 Table: Binomial Test
  71 ,,Category,N,Observed Prop.,Test Prop.,Exact Sig. (1-tailed)
  72 x,Group1,1,8,.500,.400,.284
  73 ,Group2,2,8,.500,,
  74 ,Total,,16,1.000,,
  75 ])
  76 AT_CLEANUP
  77
  78 AT_SETUP([NPAR TESTS BINOMIAL, P > 0.5; N1/N2 < 1])
  79 AT_DATA([npar.sps], [dnl
  80 SET FORMAT F8.3.
  81
  82 DATA LIST LIST NOTABLE /x (F8.0) w (F8.0).
  83 BEGIN DATA.
  84 1   11
  85 2   12
  86 END DATA.
  87
  88 WEIGHT BY w.
  89
  90 NPAR TESTS
  91         /BINOMIAL(0.6) = x
  92         .
  93 ])
  94 AT_CHECK([pspp -O format=csv npar.sps], [0], [dnl
  95 Table: Binomial Test
  96 ,,Category,N,Observed Prop.,Test Prop.,Exact Sig. (1-tailed)
  97 x,Group1,1,11,.478,.600,.164
  98 ,Group2,2,12,.522,,
  99 ,Total,,23,1.000,,
 100 ])
 101 AT_CLEANUP
 102
 103 AT_SETUP([NPAR TESTS BINOMIAL, P > 0.5; N1/N2 > 1])
 104 AT_DATA([npar.sps], [dnl
 105 SET FORMAT F8.3.
 106
 107 DATA LIST LIST NOTABLE /x (F8.0) w (F8.0).
 108 BEGIN DATA.
 109 1   11
 110 2   9
 111 END DATA.
 112
 113 WEIGHT BY w.
 114
 115 NPAR TESTS
 116         /BINOMIAL(0.6) = x.
 117 ])
 118 AT_CHECK([pspp -O format=csv npar.sps], [0], [dnl
 119 Table: Binomial Test
 120 ,,Category,N,Observed Prop.,Test Prop.,Exact Sig. (1-tailed)
 121 x,Group1,1,11,.550,.600,.404
 122 ,Group2,2,9,.450,,
 123 ,Total,,20,1.000,,
 124 ])
 125 AT_CLEANUP
 126
 127 AT_SETUP([NPAR TESTS BINOMIAL, P > 0.5; N1/N2 = 1])
 128 AT_DATA([npar.sps], [dnl
 129 SET FORMAT F8.3.
 130
 131 DATA LIST LIST NOTABLE /x (F8.0) w (F8.0).
 132 BEGIN DATA.
 133 1   11
 134 2   11
 135 END DATA.
 136
 137 WEIGHT BY w.
 138
 139 NPAR TESTS
 140         /BINOMIAL(0.6) = x.
 141 ])
 142 AT_CHECK([pspp -O format=csv npar.sps], [0], [dnl
 143 Table: Binomial Test
 144 ,,Category,N,Observed Prop.,Test Prop.,Exact Sig. (1-tailed)
 145 x,Group1,1,11,.500,.600,.228
 146 ,Group2,2,11,.500,,
 147 ,Total,,22,1.000,,
 148 ])
 149 AT_CLEANUP
 150
 151 AT_SETUP([NPAR TESTS BINOMIAL, P = 0.5; N1/N2 < 1])
 152 AT_DATA([npar.sps], [dnl
 153 SET FORMAT F8.3.
 154
 155 DATA LIST LIST NOTABLE /x (F8.0) w (F8.0).
 156 BEGIN DATA.
 157 1   8
 158 2   15
 159 END DATA.
 160
 161 WEIGHT BY w.
 162
 163 NPAR TESTS
 164         /BINOMIAL = x
 165         .
 166 ])
 167 AT_CHECK([pspp -O format=csv npar.sps], [0], [dnl
 168 Table: Binomial Test
 169 ,,Category,N,Observed Prop.,Test Prop.,Exact Sig. (2-tailed)
 170 x,Group1,1,8,.348,.500,.210
 171 ,Group2,2,15,.652,,
 172 ,Total,,23,1.000,,
 173 ])
 174 AT_CLEANUP
 175
 176 AT_SETUP([NPAR TESTS BINOMIAL, P = 0.5; N1/N2 > 1])
 177 AT_DATA([npar.sps], [dnl
 178 SET FORMAT F8.3.
 179
 180 DATA LIST LIST NOTABLE /x (F8.0) w (F8.0).
 181 BEGIN DATA.
 182 1   12
 183 2   6
 184 END DATA.
 185
 186 WEIGHT BY w.
 187
 188 NPAR TESTS
 189         /BINOMIAL(0.5) = x.
 190 ])
 191 AT_CHECK([pspp -O format=csv npar.sps], [0], [dnl
 192 Table: Binomial Test
 193 ,,Category,N,Observed Prop.,Test Prop.,Exact Sig. (2-tailed)
 194 x,Group1,1,12,.667,.500,.238
 195 ,Group2,2,6,.333,,
 196 ,Total,,18,1.000,,
 197 ])
 198 AT_CLEANUP
 199
 200 AT_SETUP([NPAR TESTS BINOMIAL, P = 0.5; N1/N2 = 1])
 201 AT_DATA([npar.sps], [dnl
 202 SET FORMAT F8.3.
 203
 204 DATA LIST LIST NOTABLE /x (F8.0) w (F8.0).
 205 BEGIN DATA.
 206 1   10
 207 2   10
 208 END DATA.
 209
 210 WEIGHT BY w.
 211
 212 NPAR TESTS
 213         /BINOMIAL(0.5) = x
 214         .
 215 ])
 216 AT_CHECK([pspp -O format=csv npar.sps], [0], [dnl
 217 Table: Binomial Test
 218 ,,Category,N,Observed Prop.,Test Prop.,Exact Sig. (2-tailed)
 219 x,Group1,1,10,.500,.500,1.000
 220 ,Group2,2,10,.500,,
 221 ,Total,,20,1.000,,
 222 ])
 223 AT_CLEANUP
 224
 225 AT_SETUP([NPAR TESTS BINOMIAL, P = 0.5; N1/N2 = 1 Cutpoint])
 226 AT_DATA([npar.sps], [dnl
 227 SET FORMAT F8.3.
 228
 229 DATA LIST LIST NOTABLE /x * w *.
 230 BEGIN DATA.
 231 9    3
 232 10   7
 233 11   16
 234 END DATA.
 235
 236 WEIGHT BY w.
 237
 238 NPAR TESTS
 239         /BINOMIAL(0.5) = x (10)
 240         .
 241 ])
 242 AT_CHECK([pspp -O format=csv npar.sps], [0], [dnl
 243 Table: Binomial Test
 244 ,,Category,N,Observed Prop.,Test Prop.,Exact Sig. (2-tailed)
 245 x,Group1,<= 10,10.000,.385,.500,.327
 246 ,Group2,,16.000,.615,,
 247 ,Total,,26.000,1.000,,
 248 ])
 249 AT_CLEANUP
 250
 251 AT_SETUP([NPAR TESTS BINOMIAL, P = 0.5; N1/N2 = 1 Named values])
 252 AT_DATA([npar.sps], [dnl
 253 SET FORMAT F8.3.
 254
 255 DATA LIST LIST NOTABLE /x * w *.
 256 BEGIN DATA.
 257 10   10
 258 15   45
 259 20   13
 260 END DATA.
 261
 262 WEIGHT BY w.
 263
 264 NPAR TESTS
 265         /BINOMIAL(0.5) = x (10, 20)
 266         .
 267 ])
 268 AT_CHECK([pspp -O format=csv npar.sps], [0], [dnl
 269 Table: Binomial Test
 270 ,,Category,N,Observed Prop.,Test Prop.,Exact Sig. (2-tailed)
 271 x,Group1,10.000,10.000,.435,.500,.678
 272 ,Group2,20.000,13.000,.565,,
 273 ,Total,,23.000,1.000,,
 274 ])
 275 AT_CLEANUP
 276
 277 AT_SETUP([NPAR TESTS CHISQUARE])
 278 AT_DATA([npar.sps], [dnl
 279 DATA LIST NOTABLE LIST /x * y * w *.
 280 BEGIN DATA.
 281 1   2  1
 282 2   1  3
 283 3.1 1  4
 284 3.2 2  1
 285 4   2  2
 286 5   3  1
 287 1   4  2
 288 END DATA.
 289
 290 WEIGHT BY w.
 291
 292 NPAR TESTS
 293   CHISQUARE=x y
 294   .
 295
 296 NPAR TESTS
 297   CHISQUARE=y
 298   /EXPECTED=3 4 5 4
 299   .
 300
 301 NPAR TESTS
 302   CHISQUARE=x y(2, 4)
 303   /EXPECTED = 6 10 3
 304   .
 305 ])
 306 AT_CHECK([pspp -O format=csv npar.sps], [0], [dnl
 307 Table: x
 308 ,Observed N,Expected N,Residual
 309 1.00,3.00,2.33,.67
 310 2.00,3.00,2.33,.67
 311 3.10,4.00,2.33,1.67
 312 3.20,1.00,2.33,-1.33
 313 4.00,2.00,2.33,-.33
 314 5.00,1.00,2.33,-1.33
 315 Total,14.00,,
 316
 317 Table: y
 318 ,Observed N,Expected N,Residual
 319 1.00,7.00,3.50,3.50
 320 2.00,4.00,3.50,.50
 321 3.00,1.00,3.50,-2.50
 322 4.00,2.00,3.50,-1.50
 323 Total,14.00,,
 324
 325 Table: Test Statistics
 326 ,x,y
 327 Chi-Square,3.14,6.00
 328 df,5,3
 329 Asymp. Sig.,.68,.11
 330
 331 Table: y
 332 ,Observed N,Expected N,Residual
 333 1.00,7.00,2.63,4.38
 334 2.00,4.00,3.50,.50
 335 3.00,1.00,4.38,-3.38
 336 4.00,2.00,3.50,-1.50
 337 Total,14.00,,
 338
 339 Table: Test Statistics
 340 ,y
 341 Chi-Square,10.61
 342 df,3
 343 Asymp. Sig.,.01
 344
 345 Table: Frequencies
 346 ,x,,,,y,,,
 347 ,Category,Observed N,Expected N,Residual,Category,Observed N,Expected N,Residual
 348 1,2.00,3.00,3.16,-.16,2.00,4.00,2.21,1.79
 349 2,3.00,5.00,5.26,-.26,3.00,1.00,3.68,-2.68
 350 3,4.00,2.00,1.58,.42,4.00,2.00,1.11,.89
 351 Total,,10.00,,,,7.00,,
 352
 353 Table: Test Statistics
 354 ,x,y
 355 Chi-Square,.13,4.13
 356 df,2,2
 357 Asymp. Sig.,.94,.13
 358 ])
 359 AT_CLEANUP
 360
 361 AT_SETUP([NPAR TESTS CHISQUARE expected values missing])
 362 AT_DATA([npar.sps], [dnl
 363 DATA LIST NOTABLE LIST /x * y * w *.
 364 BEGIN DATA.
 365 1   2  1
 366 2   1  3
 367 3.1 1  4
 368 3.2 2  1
 369 4   2  2
 370 5   3  1
 371 1   4  2
 372 END DATA.
 373
 374 WEIGHT BY w.
 375
 376 NPAR TESTS
 377   CHISQUARE=y
 378   /EXPECTED = 3 4 5 4 3 1
 379   .
 380 ])
 381 AT_CHECK([pspp -O format=csv npar.sps], [1], [dnl
 382 "error: CHISQUARE test specified 6 expected values, but 4 distinct values were encountered in variable y."
 383
 384 Table: Test Statistics
 385 ,y
 386 Chi-Square,.00
 387 df,0
 388 Asymp. Sig.,1.00
 389 ])
 390 AT_CLEANUP
 391
 392 AT_SETUP([NPAR TESTS CHISQUARE with DESCRIPTIVES])
 393 AT_DATA([npar.sps], [dnl
 394 DATA LIST NOTABLE LIST /x * y * w * .
 395 BEGIN DATA.
 396 1   2  1
 397 2   1  3
 398 3.1 1  4
 399 3.2 2  1
 400 4   2  2
 401 5   3  1
 402 1   4  2
 403 .   5  1
 404 END DATA.
 405
 406 WEIGHT BY w.
 407
 408 MISSING VALUES x (4).
 409
 410 NPAR TESTS
 411   CHISQUARE=x y(-2,5)
 412   /MISSING=ANALYSIS
 413   /STATISTICS=DESCRIPTIVES
 414   .
 415 ])
 416 AT_CHECK([pspp -O format=csv npar.sps], [0], [dnl
 417 Table: Frequencies
 418 ,x,,,,y,,,
 419 ,Category,Observed N,Expected N,Residual,Category,Observed N,Expected N,Residual
 420 1,-2.00,.00,1.50,-1.50,-2.00,.00,1.88,-1.88
 421 2,-1.00,.00,1.50,-1.50,-1.00,.00,1.88,-1.88
 422 3,.00,.00,1.50,-1.50,.00,.00,1.88,-1.88
 423 4,1.00,3.00,1.50,1.50,1.00,7.00,1.88,5.13
 424 5,2.00,3.00,1.50,1.50,2.00,4.00,1.88,2.13
 425 6,3.00,5.00,1.50,3.50,3.00,1.00,1.88,-.88
 426 7,4.00,.00,1.50,-1.50,4.00,2.00,1.88,.13
 427 8,5.00,1.00,1.50,-.50,5.00,1.00,1.88,-.88
 428 Total,,12.00,,,,15.00,,
 429
 430 Table: Test Statistics
 431 ,x,y
 432 Chi-Square,17.33,22.87
 433 df,7,7
 434 Asymp. Sig.,.02,.00
 435
 436 Table: Descriptive Statistics
 437 ,N,Mean,Std. Deviation,Minimum,Maximum
 438 ,,,,,
 439 x,12.00,2.47,1.19,1.00,5.00
 440 y,15.00,2.07,1.33,1.00,5.00
 441 ])
 442 AT_CLEANUP
 443
 444 AT_SETUP([NPAR TESTS CHISQUARE, listwise missing])
 445 AT_DATA([npar.sps], [dnl
 446 DATA LIST NOTABLE LIST /x * y * w * .
 447 BEGIN DATA.
 448 1   2  1
 449 2   1  3
 450 3.1 1  4
 451 3.2 2  1
 452 4   2  2
 453 5   3  1
 454 1   4  2
 455 .   5  1
 456 END DATA.
 457
 458 WEIGHT BY w.
 459
 460 * MISSING VALUES x (4).
 461
 462 NPAR TESTS
 463   CHISQUARE=x y(-2,5)
 464   /MISSING=LISTWISE
 465   /STATISTICS=DESCRIPTIVES
 466   .
 467 ])
 468 AT_CHECK([pspp -O format=csv npar.sps], [0], [dnl
 469 Table: Frequencies
 470 ,x,,,,y,,,
 471 ,Category,Observed N,Expected N,Residual,Category,Observed N,Expected N,Residual
 472 1,-2.00,.00,1.75,-1.75,-2.00,.00,1.75,-1.75
 473 2,-1.00,.00,1.75,-1.75,-1.00,.00,1.75,-1.75
 474 3,.00,.00,1.75,-1.75,.00,.00,1.75,-1.75
 475 4,1.00,3.00,1.75,1.25,1.00,7.00,1.75,5.25
 476 5,2.00,3.00,1.75,1.25,2.00,4.00,1.75,2.25
 477 6,3.00,5.00,1.75,3.25,3.00,1.00,1.75,-.75
 478 7,4.00,2.00,1.75,.25,4.00,2.00,1.75,.25
 479 8,5.00,1.00,1.75,-.75,5.00,.00,1.75,-1.75
 480 Total,,14.00,,,,14.00,,
 481
 482 Table: Test Statistics
 483 ,x,y
 484 Chi-Square,13.43,26.00
 485 df,7,7
 486 Asymp. Sig.,.06,.00
 487
 488 Table: Descriptive Statistics
 489 ,N,Mean,Std. Deviation,Minimum,Maximum
 490 ,,,,,
 491 x,14.00,2.69,1.23,1.00,5.00
 492 y,14.00,1.86,1.10,1.00,4.00
 493 ])
 494 AT_CLEANUP
 495
 496 AT_SETUP([NPAR TESTS WILCOXON])
 497 AT_DATA([npar.sps], [dnl
 498 data list notable list /foo * bar * w (f8.0).
 499 begin data.
 500 1.00     1.00   1
 501 1.00     2.00   1
 502 2.00     1.00   1
 503 1.00     4.00   1
 504 2.00     5.00   1
 505 1.00    19.00   1
 506 2.00     7.00   1
 507 4.00     5.00   1
 508 1.00    12.00   1
 509 2.00    13.00   1
 510 2.00     2.00   1
 511 12.00      .00  2
 512 12.00     1.00  1
 513 13.00     1.00  1
 514 end data
 515
 516 variable labels foo "first" bar "second".
 517
 518 weight by w.
 519
 520 npar test
 521  /wilcoxon=foo with bar (paired)
 522  /missing analysis
 523  /method=exact.
 524 ])
 525 AT_CHECK([pspp -o pspp.csv npar.sps])
 526 AT_CHECK([cat pspp.csv], [0], [dnl
 527 Table: Ranks
 528 ,,N,Mean Rank,Sum of Ranks
 529 second - first,Negative Ranks,5,8.60,43.00
 530 ,Positive Ranks,8,6.00,48.00
 531 ,Ties,2,,
 532 ,Total,15,,
 533
 534 Table: Test Statistics
 535 ,second - first
 536 Z,-.18
 537 Asymp. Sig. (2-tailed),.86
 538 Exact Sig. (2-tailed),.89
 539 Exact Sig. (1-tailed),.45
 540 ])
 541 AT_CLEANUP
 542
 543 AT_SETUP([NPAR TESTS WILCOXON with missing values])
 544 AT_DATA([npar.sps], [dnl
 545 data list notable list /foo * bar * dummy *.
 546 begin data.
 547 1.00     1.00    1
 548 1.00     2.00    1
 549 2.00     1.00    1
 550 1.00     4.00    .
 551 2.00     5.00    .
 552 1.00    19.00    .
 553 2.00     7.00    1
 554 4.00     5.00    1
 555 1.00    12.00    1
 556 2.00    13.00    1
 557 2.00     2.00    1
 558 12.00      .00   1
 559 12.00      .00   1
 560 34.2       .     1
 561 12.00     1.00   1
 562 13.00     1.00   1
 563 end data
 564
 565 variable labels foo "first" bar "second".
 566
 567 npar test
 568  /wilcoxon=foo with bar (paired)
 569  /missing analysis
 570  /method=exact.
 571
 572 ])
 573 AT_CHECK([pspp -o pspp.csv npar.sps])
 574 dnl This is the same output as the previous test.
 575 AT_CHECK([cat pspp.csv], [0], [dnl
 576 Table: Ranks
 577 ,,N,Mean Rank,Sum of Ranks
 578 second - first,Negative Ranks,5,8.60,43.00
 579 ,Positive Ranks,8,6.00,48.00
 580 ,Ties,2,,
 581 ,Total,15,,
 582
 583 Table: Test Statistics
 584 ,second - first
 585 Z,-.18
 586 Asymp. Sig. (2-tailed),.86
 587 Exact Sig. (2-tailed),.89
 588 Exact Sig. (1-tailed),.45
 589 ])
 590 AT_CLEANUP
 591
 592 AT_SETUP([NPAR TESTS SIGN])
 593 AT_DATA([npar.sps], [dnl
 594 set format = F9.3.
 595
 596 data list notable list /age * height rank *.
 597 begin data.
 598 10 12 11
 599 12 13 13
 600 13 14 12
 601 12 12 10
 602 9   9 10
 603 10.3 10.2 12
 604 end data.
 605
 606 npar tests
 607         /sign=age height WITH height rank (PAIRED)
 608         /MISSING ANALYSIS
 609         /METHOD=EXACT
 610         .
 611 ])
 612 AT_CHECK([pspp -o pspp.csv npar.sps])
 613 dnl Some machines return .313 instead of .312 for the Point Probability
 614 dnl (see bug #31611).
 615 AT_CHECK([sed 's/\.313$/.312/' pspp.csv], [0], [dnl
 616 Table: Frequencies
 617 ,,N
 618 height - age,Negative Differences,1
 619 ,Positive Differences,3
 620 ,Ties,2
 621 ,Total,6
 622 rank - height,Negative Differences,3
 623 ,Positive Differences,2
 624 ,Ties,1
 625 ,Total,6
 626
 627 Table: Test Statistics
 628 ,height - age,rank - height
 629 Exact Sig. (2-tailed),.625,1.000
 630 Exact Sig. (1-tailed),.312,.500
 631 Point Probability,.250,.312
 632 ])
 633 AT_CLEANUP
 634
 635
 636 AT_SETUP([NPAR Kruskal-Wallis test])
 637
 638 dnl Simple case
 639 AT_DATA([kw-simple.sps], [dnl
 640 set format = F9.3.
 641
 642 data list notable list /gv * xscore *.
 643 begin data
 644 1 96
 645 1 128
 646 1 83
 647 2 132
 648 2 135
 649 2 109
 650 3 115
 651 1 61
 652 1 101
 653 2 82
 654 2 124
 655 3 149
 656 3 166
 657 3 147
 658 end data.
 659
 660 value label /gv
 661        1 "timed out"
 662        2 "hit wicket"
 663        3 "handled the ball".
 664
 665 npar tests
 666         /kruskal-wallis xscore by gv (1, 3)
 667         .
 668 ])
 669
 670 AT_CHECK([pspp -o pspp.csv kw-simple.sps])
 671 AT_CHECK([cat pspp.csv], [0], [dnl
 672 Table: Ranks
 673 ,gv,N,Mean Rank
 674 xscore,timed out,5,4.400
 675 ,handled the ball,4,11.500
 676 ,hit wicket,5,7.400
 677 ,Total,14,
 678
 679 Table: Test Statistics
 680 ,,xscore
 681 Chi-Square,,6.406
 682 df,,2
 683 Asymp. Sig.,,.041
 684 ])
 685
 686
 687 dnl Now try a missing value in the group variable
 688 AT_DATA([kw-missing-group.sps], [dnl
 689 set format = F9.3.
 690
 691 data list notable list /gv * xscore *.
 692 begin data
 693 1 96
 694 1 128
 695 1 83
 696 1 61
 697 1 101
 698 2 82
 699 2 124
 700 2 132
 701 2 135
 702 2 109
 703 3 115
 704 3 149
 705 3 166
 706 3 147
 707 2.5 344
 708 end data.
 709
 710 missing values gv (2.5).
 711
 712 value label /gv
 713        1 "timed out"
 714        2 "hit wicket"
 715        3 "handled the ball".
 716
 717 npar tests
 718         /kruskal-wallis xscore by gv (1, 3)
 719         /missing=exclude
 720         .
 721 ])
 722
 723 AT_CHECK([pspp -o pspp2.csv kw-missing-group.sps])
 724
 725 dnl The result should be the same as before
 726 AT_CHECK([diff pspp.csv pspp2.csv], [0])
 727
 728 AT_CLEANUP
 729
 730
 731 AT_SETUP([NPAR Kruskal-Wallis multiple-variables])
 732
 733 AT_DATA([kw-multi.sps], [dnl
 734 set format = F9.3.
 735
 736 data list notable list /gv * xscore * yscore.
 737 begin data
 738 1 96   .
 739 1 128  .
 740 1 83   .
 741 2 132  132
 742 2 135  135
 743 2 109  109
 744 3 115  115
 745 1 61   .
 746 1 101  .
 747 2 82   82
 748 2 124  124
 749 3 149  149
 750 3 166  166
 751 3 147  147
 752 4 .    96
 753 4 .    128
 754 4 .    83
 755 4 .    61
 756 4 .    101
 757 end data.
 758
 759 value label /gv
 760        1 "timed out"
 761        2 "hit wicket"
 762        3 "handled the ball"
 763        4 "bowled"
 764        5 "lbw"
 765        .
 766
 767 npar tests
 768         /k-w xscore yscore by gv (1, 5)
 769         .
 770
 771 ])
 772
 773
 774 AT_CHECK([pspp -o pspp.csv kw-multi.sps])
 775 AT_CHECK([cat pspp.csv], [0], [dnl
 776 Table: Ranks
 777 ,gv,N,Mean Rank
 778 xscore,timed out,5,4.400
 779 ,handled the ball,4,11.500
 780 ,hit wicket,5,7.400
 781 ,Total,14,
 782 yscore,handled the ball,4,11.500
 783 ,bowled,5,4.400
 784 ,hit wicket,5,7.400
 785 ,Total,14,
 786
 787 Table: Test Statistics
 788 ,,xscore,yscore,
 789 Chi-Square,,6.406,6.406,
 790 df,,2,2,
 791 Asymp. Sig.,,.041,.041,
 792 ])
 793
 794 AT_CLEANUP
 795
 796
 797
 798 AT_SETUP([NPAR TESTS Runs])
 799 AT_DATA([npar-runs.sps], [dnl
 800 set format F11.4.
 801 data list notable list /score * w *.
 802 begin data
 803 4     6
 804 .     4
 805 4     3
 806 3    20
 807 2    29
 808 1    42
 809 6    18
 810 5     7
 811 6    78
 812 5    10
 813 6    46
 814 5     5
 815 6    17
 816 5     1
 817 6    11
 818 4     2
 819 3     7
 820 2     6
 821 1    10
 822 4    13
 823 3    22
 824 3    11
 825 2    24
 826 1    18
 827 4     4
 828 3    12
 829 2    10
 830 1    25
 831 4     4
 832 3     7
 833 2     3
 834 1     4
 835 4     2
 836 3     3
 837 2     2
 838 1     4
 839 end data.
 840
 841 weight by w.
 842
 843 npar tests
 844         /runs (MEDIAN) = score
 845         /runs (MEAN) = score
 846         /runs (MODE) = score
 847         .
 848 ])
 849
 850 AT_CHECK([pspp -o pspp.csv npar-runs.sps])
 851
 852 AT_CHECK([cat pspp.csv], [0], [dnl
 853 Table: Runs Test
 854 ,score
 855 Test Value (median),3.0000
 856 Cases < Test Value,177.0000
 857 Cases >= Test Value,309.0000
 858 Total Cases,486.0000
 859 Number of Runs,12
 860 Z,-20.9931
 861 Asymp. Sig. (2-tailed),.0000
 862
 863 Table: Runs Test
 864 ,score
 865 Test Value (mean),3.6379
 866 Cases < Test Value,259.0000
 867 Cases >= Test Value,227.0000
 868 Total Cases,486.0000
 869 Number of Runs,12
 870 Z,-21.0650
 871 Asymp. Sig. (2-tailed),.0000
 872
 873 Table: Runs Test
 874 ,score
 875 Test Value (mode),6.0000
 876 Cases < Test Value,316.0000
 877 Cases >= Test Value,170.0000
 878 Total Cases,486.0000
 879 Number of Runs,11
 880 Z,-21.0742
 881 Asymp. Sig. (2-tailed),.0000
 882 ])
 883
 884 AT_CLEANUP
 885
 886
 887 AT_SETUP([NPAR TESTS Friedman])
 888 AT_DATA([npar-friedman.sps], [dnl
 889 set format F15.4.
 890 data list notable list /x * y * z.
 891 begin data
 892 9.5 6.5 8.1
 893 8.0 6.0 6.0
 894 7.0 6.5 4.2
 895 9.5 5.0 7.3
 896 9.0 7.0 6.2
 897 8.5 6.9 6.5
 898 7.5 8.0 6.5
 899 6.0 8.0 3.1
 900 5.0 6.0 4.9
 901 7.5 7.5 6.2
 902 end data.
 903
 904 npar tests
 905      /friedman = x y z.
 906 ])
 907
 908 AT_CHECK([pspp -o pspp.csv npar-friedman.sps])
 909
 910 AT_CHECK([cat pspp.csv], [0], [dnl
 911 Table: Ranks
 912 ,Mean Rank
 913 x,2.6500
 914 y,2.1000
 915 z,1.2500
 916
 917 Table: Test Statistics
 918 N,10
 919 Chi-Square,10.4737
 920 df,2
 921 Asymp. Sig.,.0053
 922 ])
 923
 924 AT_CLEANUP
 925
 926
 927
 928 AT_SETUP([NPAR TESTS Mann-Whitney])
 929 AT_DATA([npar-mann-whitney.sps], [dnl
 930 SET FORMAT     = F11.4
 931
 932 data list notable list /height * sex (f1.0).
 933 begin data.
 934 201 1
 935 84 1
 936 83 1
 937 94 1
 938 88 0
 939 99 0
 940 55 0
 941 69 0
 942 86 1
 943 79 1
 944 91 0
 945 201 0
 946 88 1
 947 85 1
 948 82 1
 949 88 0
 950 75 0
 951 99 0
 952 81 0
 953 72 1
 954 89 1
 955 92 1
 956 80 0
 957 82 0
 958 76 0
 959 65 0
 960 85 0
 961 76 1
 962 145 1
 963 24 1
 964 end data.
 965
 966 NPAR TESTS
 967      /M-W = height BY sex (0,1).
 968 ])
 969
 970 AT_CHECK([pspp -o pspp.csv npar-mann-whitney.sps])
 971
 972 AT_CHECK([cat pspp.csv], [0], [dnl
 973 Table: Ranks
 974 ,N,,,Mean Rank,,Sum of Ranks,
 975 ,0,1,Total,0,1,0,1
 976 height,15.0000,15.0000,30.0000,14.5333,16.4667,218.0000,247.0000
 977
 978 Table: Test Statistics
 979 ,Mann-Whitney U,Wilcoxon W,Z,Asymp. Sig. (2-tailed)
 980 height,98.0000,218.0000,-.6020,.5472
 981 ])
 982
 983
 984 AT_CLEANUP
 985
 986
 987 AT_SETUP([NPAR TESTS Cochran])
 988 AT_DATA([npar-cochran.sps], [dnl
 989 set format f11.3.
 990
 991 data list notable list /v1 * v2 * v3 * v4 * v5 * v6 * v7 *.
 992 begin data.
 993 2 1 1 2 1 1 2
 994 2 2 2 2 1 1 1
 995 1 1 2 2 1 1 2
 996 2 2 2 2 1 1 2
 997 2 1 2 1 1 2 1
 998 1 2 2 1 1 1 1
 999 1 2 2 2 2 2 2
1000 2 2 1 2 1 1 1
1001 1 2 1 2 1 1 2
1002 end data.
1003
1004 npar tests
1005         /cochran = v1 to v7 .
1006
1007 ])
1008
1009 AT_CHECK([pspp -o pspp.csv npar-cochran.sps])
1010
1011 AT_CHECK([cat pspp.csv], [0], [dnl
1012 Table: Frequencies
1013 ,Value,
1014 ,Success (2),Failure (1)
1015 v1,5,4
1016 v2,6,3
1017 v3,6,3
1018 v4,7,2
1019 v5,1,8
1020 v6,2,7
1021 v7,5,4
1022
1023 Table: Test Statistics
1024 N,9
1025 Cochran's Q,12.735
1026 df,6
1027 Asymp. Sig.,.047
1028 ])
1029
1030 AT_CLEANUP
1031
1032
1033
1034 AT_SETUP([NPAR TESTS Kendall])
1035 AT_DATA([npar-kendall.sps], [dnl
1036 SET FORMAT F14.3.
1037
1038 data list notable list /v1 * v2 * v3
1039 begin data.
1040  7  7  2
1041  5  6  5
1042  8  6  4
1043  5  7  4
1044  5  4  4
1045  8  6  5
1046  6  3  5
1047  7  6  5
1048  8  5  5
1049  .  2  2
1050  5  4  5
1051  3  4  4
1052  5  1  2
1053  5  2  1
1054  7  6  5
1055  6  3  4
1056  6  6  6
1057  5  4  5
1058  4  3  4
1059  9  1  1
1060  6  2  1
1061  3  7  8
1062  6  3  4
1063  4  4  4
1064  5  4  3
1065  6  5  2
1066  4  4  8
1067  4  6  4
1068  6  5  5
1069  7  8  6
1070  5  3  5
1071 end data.
1072
1073 npar tests
1074         /kendall = all
1075         .
1076 ])
1077
1078 AT_CHECK([pspp -o pspp.csv npar-kendall.sps])
1079
1080 AT_CHECK([cat pspp.csv], [0], [dnl
1081 Table: Ranks
1082 ,Mean Rank
1083 v1,2.500
1084 v2,1.817
1085 v3,1.683
1086
1087 Table: Test Statistics
1088 N,30
1089 Kendall's W,.233
1090 Chi-Square,13.960
1091 df,2
1092 Asymp. Sig.,.001
1093 ])
1094
1095 AT_CLEANUP