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