pintos-os.org Git - pspp/blob - tests/language/stats/regression.at

   1 AT_BANNER([LINEAR REGRESSION])
   2
   3 AT_SETUP([LINEAR REGRESSION - basic])
   4 AT_DATA([regression.sps], [dnl
   5 set format = F22.3.
   6 data list notable list / v0 to v2.
   7 begin data
   8  0.65377128  7.735648 -23.97588
   9 -0.13087553  6.142625 -19.63854
  10  0.34880368  7.651430 -25.26557
  11  0.69249021  6.125125 -16.57090
  12 -0.07368178  8.245789 -25.80001
  13 -0.34404919  6.031540 -17.56743
  14  0.75981559  9.832291 -28.35977
  15 -0.46958313  5.343832 -16.79548
  16 -0.06108490  8.838262 -29.25689
  17  0.56154863  6.200189 -18.58219
  18 end data
  19 regression /variables=v0 v1 v2 /statistics defaults /dependent=v2 /method=enter /save=pred resid.
  20 list.
  21 ])
  22
  23 AT_CHECK([pspp -O format=csv regression.sps], [0], [dnl
  24 Table: Model Summary (v2)
  25 ,R,R Square,Adjusted R Square,Std. Error of the Estimate
  26 ,.971,.942,.925,1.337
  27
  28 Table: ANOVA (v2)
  29 ,,Sum of Squares,df,Mean Square,F,Sig.
  30 ,Regression,202.753,2,101.376,56.754,.000
  31 ,Residual,12.504,7,1.786,,
  32 ,Total,215.256,9,,,
  33
  34 Table: Coefficients (v2)
  35 ,,Unstandardized Coefficients,,Standardized Coefficients,,
  36 ,,B,Std. Error,Beta,t,Sig.
  37 ,(Constant),2.191,2.357,.000,.930,.380
  38 ,v0,1.813,1.053,.171,1.722,.129
  39 ,v1,-3.427,.332,-1.026,-10.334,.000
  40
  41 Table: Data List
  42 v0,v1,v2,RES1,PRED1
  43 .654,7.736,-23.976,-.84,-23.13
  44 -.131,6.143,-19.639,-.54,-19.10
  45 .349,7.651,-25.266,-1.87,-23.40
  46 .692,6.125,-16.571,.97,-17.54
  47 -.074,8.246,-25.800,.40,-26.20
  48 -.344,6.032,-17.567,1.53,-19.10
  49 .760,9.832,-28.360,1.77,-30.13
  50 -.470,5.344,-16.795,.18,-16.97
  51 -.061,8.838,-29.257,-1.05,-28.21
  52 .562,6.200,-18.582,-.54,-18.04
  53 ])
  54 AT_CLEANUP
  55
  56
  57 # Test to ensure that the /SAVE subcommand works properly when SPLIT is active
  58 AT_SETUP([LINEAR REGRESSION - SAVE vs SPLITS])
  59
  60 # Generate some test data based on a linear model
  61 AT_DATA([gen-data.sps], [dnl
  62 set seed = 1.
  63 input program.
  64 loop #c = 1 to 20.
  65      compute x0 = rv.normal (0,1).
  66      compute x1 = rv.normal (0,2).
  67      compute err = rv.normal (0,0.1).
  68      compute y = 4 - 2 * x0 + 3 * x1 + err.
  69      compute g = (#c > 10).
  70      end case.
  71 end loop.
  72 end file.
  73 end input program.
  74
  75 print outfile='regdata.txt' /g x0 x1 y err *.
  76 execute.
  77 ])
  78
  79 AT_CHECK([pspp -O format=csv gen-data.sps], [0], [ignore])
  80
  81 # Use our test data to create a predictor and a residual variable
  82 # for G == 0
  83 AT_DATA([regression0.sps], [dnl
  84 data list notable file='regdata.txt' list /g x0 x1 y err *.
  85
  86 select if (g = 0).
  87
  88 regression
  89            /variables = x0 x1
  90            /dependent = y
  91            /statistics = all
  92            /save = pred resid.
  93            .
  94
  95 print outfile='outdata-g0.txt' /g x0 x1 y err res1 pred1 *.
  96 execute.
  97 ])
  98
  99
 100 AT_CHECK([pspp -O format=csv regression0.sps], [0], [ignore])
 101
 102 # Use our test data to create a predictor and a residual variable
 103 # for G == 1
 104 AT_DATA([regression1.sps], [dnl
 105 data list notable file='regdata.txt' list /g x0 x1 y err *.
 106
 107 select if (g = 1).
 108
 109 regression
 110            /variables = x0 x1
 111            /dependent = y
 112            /statistics = all
 113            /save = pred resid.
 114            .
 115
 116 print outfile='outdata-g1.txt' /g x0 x1 y err res1 pred1 *.
 117 execute.
 118 ])
 119
 120
 121 AT_CHECK([pspp -O format=csv regression1.sps], [0], [ignore])
 122
 123 # Use our test data to create a predictor and a residual variable
 124 # The data is split on G
 125 AT_DATA([regression-split.sps], [dnl
 126 data list notable file='regdata.txt' list /g x0 x1 y err *.
 127
 128 split file by g.
 129
 130 regression
 131            /variables = x0 x1
 132            /dependent = y
 133            /statistics = all
 134            /save = pred resid.
 135            .
 136
 137 print outfile='outdata-split.txt' /g x0 x1 y err res1 pred1 *.
 138 execute.
 139 ])
 140
 141 AT_CHECK([pspp -O format=csv regression-split.sps], [0], [ignore])
 142
 143 # The concatenation of G==0 and G==1 should be identical to the SPLIT data
 144 AT_CHECK([cat outdata-g0.txt outdata-g1.txt | diff outdata-split.txt - ], [0], [])
 145
 146 AT_CLEANUP
 147
 148
 149 # Test that the procedure behaves sensibly when presented with
 150 # multiple dependent variables
 151 AT_SETUP([LINEAR REGRESSION multiple dependent variables])
 152 AT_DATA([regression.sps], [dnl
 153 set seed = 2.
 154 input program.
 155 loop #c = 1 to 200.
 156      compute x0 = rv.normal (0, 1).
 157      compute x1 = rv.normal (0, 2).
 158      compute err = rv.normal (0, 0.8).
 159      compute y = 2 - 1.5 * x0 + 8.4 * x1 + err.
 160      compute ycopy = y.
 161      end case.
 162 end loop.
 163 end file.
 164 end input program.
 165
 166 regression
 167            /variables = x0 x1
 168            /dependent = y ycopy
 169            /statistics = default.
 170 ])
 171
 172 AT_CHECK([pspp -O format=csv regression.sps > output], [0], [ignore])
 173
 174 AT_CHECK([head -16 output > first], [0], [])
 175 AT_CHECK([tail -16 output > second], [0], [])
 176
 177 AT_CHECK([sed -e 's/ycopy/y/g' second | diff first -], [0], [])
 178
 179
 180 AT_CLEANUP
 181
 182 # Tests the QR decomposition used by the REGRESSION command.
 183 AT_SETUP([LINEAR REGRESSION test of QR decomposition])
 184 AT_DATA([regression.sps], [dnl
 185 data list list / v0 to v1.
 186 begin data
 187 -12.84099361 0.873270778
 188  16.64932538 0.371315664
 189  -1.88061907 0.505503722
 190  -6.20952354 0.734698282
 191   0.33272576 0.891224610
 192  -5.54912717 0.052318165
 193   6.11832417 0.448853404
 194  11.78124974 0.470447593
 195   0.75960353 0.565082303
 196   6.06432768 0.149316743
 197  -2.64919436 0.752532411
 198 -10.32250712 0.798263603
 199   2.06355038 0.469129797
 200  -9.71851742 0.927162270
 201   4.65582553 0.250629262
 202   9.54574474 0.847032310
 203   7.35544368 0.197028541
 204  -2.09609740 0.400584261
 205  10.30101161 0.671546480
 206  -5.24501039 0.929962876
 207   1.73412473 0.758161354
 208  -3.12732732 0.569785505
 209  12.66261501 0.630640223
 210  -2.90956805 0.576067804
 211   4.89649177 0.624483995
 212  13.64613114 0.591089881
 213  14.03198397 0.544587572
 214   2.23566810 0.967898139
 215   5.37367760 0.916246929
 216   9.01346888 0.451702743
 217   0.75378683 0.235544137
 218  -3.47470624 0.742668194
 219  -1.02063266 0.860311687
 220  -2.67132813 0.082460702
 221  23.67661680 0.932553932
 222   7.95061359 0.430161125
 223   2.05300558 0.066331375
 224  -2.01332644 0.163705417
 225  20.00663784 0.587292630
 226   3.06099417 0.161411889
 227  -3.46115358 0.216684625
 228  -6.85287183 0.548714855
 229  -4.27923809 0.630997663
 230  -0.94863395 0.880612945
 231   4.47481747 0.359885215
 232 -12.80962955 0.886070341
 233   9.35753086 0.187176558
 234   2.81002235 0.063035095
 235   0.01532424 0.964327101
 236   0.29867732 0.866408063
 237  -2.89035649 0.812135868
 238   4.17352811 0.608884061
 239  18.15502183 0.920568258
 240  -2.92662792 0.550792959
 241  -6.08090449 0.965036595
 242  -1.09135397 0.862548019
 243   7.02816784 0.042277017
 244 -21.20245068 0.430673493
 245  -8.83397584 0.724976162
 246  -0.89055843 0.017934904
 247   7.03871587 0.308829557
 248   3.84286316 0.685105924
 249   4.50280692 0.447635420
 250  11.39207346 0.875177896
 251  10.86673874 0.518530912
 252   7.09853081 0.588367569
 253 -12.82864915 0.184667098
 254  13.74888760 0.610891139
 255   0.37379146 0.557720134
 256  -9.79020267 0.942839981
 257   0.71574466 0.564570338
 258 -17.56040637 0.182061777
 259   2.52620466 0.306875011
 260   5.37718673 0.366807049
 261  -1.83964300 0.465772898
 262   6.04848363 0.644501799
 263   4.57402403 0.121419591
 264   8.55606848 0.373011464
 265  -8.46827907 0.491176571
 266  -1.77989798 0.734722847
 267  -0.68661121 0.540984182
 268   1.55798880 0.822587656
 269   5.22810831 0.333747878
 270   9.50280477 0.068100934
 271  -3.74521465 0.248537644
 272   1.36045068 0.851827791
 273   4.41604088 0.197207162
 274  -3.72568327 0.726916693
 275  -5.36123334 0.906513529
 276   3.61594583 0.414340595
 277 -10.01952852 0.140372658
 278  25.48681482 0.354309660
 279  -3.34529093 0.090075388
 280 -18.00437582 0.461438059
 281  -5.29782460 0.004362856
 282   2.79608522 0.861294398
 283  -1.64076209 0.345775481
 284   6.82802334 0.137933862
 285  -0.45416818 0.404379208
 286  -1.66868582 0.797685201
 287 -10.02820292 0.075876582
 288   5.68232031 0.404815042
 289   8.25113850 0.769173748
 290  -2.83544237 0.076583474
 291   0.87659945 0.092751009
 292   6.60270870 0.530444351
 293 -12.63924989 0.362099960
 294  -6.24451253 0.641993458
 295   3.53339015 0.461991892
 296  -0.74012232 0.437409755
 297  15.37311996 0.974913038
 298  -8.09464797 0.543308711
 299  -9.61320222 0.221564578
 300   0.21843662 0.856512540
 301  -1.56958954 0.610709221
 302   6.44977372 0.200382138
 303 -13.29136274 0.093222309
 304   6.46257214 0.024135196
 305  -3.82727990 0.601335801
 306   0.43081953 0.268230667
 307  19.06654416 0.219972815
 308  17.02906651 0.996849502
 309 -10.18073139 0.012543080
 310  12.72088788 0.910600764
 311  10.45328185 0.331285901
 312   7.14370922 0.896312020
 313  -2.81754334 0.048741266
 314   6.40217095 0.075796756
 315  -3.18030478 0.666325307
 316   8.64585957 0.120549153
 317   1.37952764 0.899991932
 318 -11.81143886 0.601949630
 319   0.03899706 0.363808260
 320 -10.63828243 0.031092967
 321  -6.66940972 0.246204205
 322  -5.07374962 0.951272057
 323   4.82281566 0.063928187
 324 -21.93693564 0.050972680
 325  -4.54569883 0.225839693
 326  -0.92422779 0.437796785
 327  -1.11683029 0.740215139
 328  16.77765554 0.851072372
 329   9.73614597 0.388180586
 330  14.05345168 0.063760129
 331   1.20512012 0.665964184
 332   8.00307080 0.102447114
 333   8.01252623 0.580929209
 334 -13.54924183 0.438420739
 335   9.87164361 0.970859344
 336  17.63437095 0.250501797
 337  -3.42503574 0.873290220
 338  -2.45873197 0.847756049
 339  17.29212092 0.411683187
 340   1.15496098 0.530658504
 341  -2.14438907 0.592255367
 342  -1.79942021 0.517773009
 343  -1.30677990 0.830860762
 344   1.70233874 0.291826660
 345  -3.05532536 0.801767829
 346  -4.06732625 0.092294501
 347   6.34665476 0.270426235
 348   9.46946411 0.196915311
 349  14.50919907 0.480357167
 350   8.93767237 0.778228613
 351   1.90298854 0.903146151
 352  18.50500507 0.598561307
 353   4.45123027 0.555898218
 354  11.37344114 0.616557707
 355 -12.14693218 0.409187285
 356  18.27198688 0.141619222
 357  -5.75939569 0.056989619
 358  -4.05515382 0.369281201
 359  16.69882098 0.946885257
 360   6.39050536 0.679704228
 361   4.04213339 0.662792380
 362   6.89608366 0.419877433
 363   1.56496633 0.358227958
 364   5.16679947 0.095144366
 365  -3.06280456 0.883265975
 366   2.76279175 0.866571973
 367   1.84969249 0.264869828
 368  21.79840498 0.702650979
 369   1.42450528 0.719308635
 370   0.96797046 0.111937435
 371  18.26840323 0.075621738
 372  13.38288377 0.573399086
 373   2.41101500 0.766238677
 374   3.83866337 0.499888953
 375  -1.56577367 0.695244089
 376  -0.90342790 0.671654151
 377  10.83775583 0.026041124
 378  -9.89767935 0.745297991
 379  11.74840150 0.309144074
 380   1.73069359 0.814063985
 381  -5.27966183 0.591005828
 382   3.33030043 0.559401806
 383   1.31427975 0.520950237
 384 -10.04588558 0.507008362
 385  10.41228345 0.425867272
 386   1.71961097 0.595783108
 387 -17.54904427 0.328788939
 388  -2.23545419 0.223377350
 389  -8.68774333 0.980964240
 390  -3.48048220 0.008877675
 391  -3.69635326 0.090236718
 392   9.76114237 0.769375983
 393 -10.25662038 0.508137553
 394   0.11155446 0.468504431
 395  -8.06824580 0.414098962
 396   3.10031660 0.327130207
 397  -3.33393146 0.756896774
 398  -3.96276749 0.530956360
 399  14.53610268 0.846474699
 400   1.70505918 0.754662464
 401  -1.93495001 0.656650411
 402   5.01974522 0.745337633
 403  13.41249973 0.489362476
 404  11.49288744 0.335924476
 405  12.59019763 0.155560469
 406 -10.17947298 0.677318449
 407   0.05556115 0.655090105
 408   3.82092860 0.051838719
 409   8.23041456 0.918272190
 410  -0.50314649 0.772015826
 411  20.05162157 0.880265258
 412   8.98816884 0.666646668
 413  -6.28312120 0.138534416
 414   3.68589909 0.274559458
 415   0.59699510 0.253180863
 416  -2.74783135 0.983525221
 417   0.32515065 0.839969577
 418  -3.60606166 0.330646732
 419  -0.82037740 0.129591173
 420   6.12444860 0.098536516
 421  10.95671074 0.033546728
 422  -2.84911174 0.720288722
 423   6.04597572 0.577061422
 424  -0.60147150 0.674096868
 425  -5.30458364 0.291468008
 426   2.68044943 0.379853840
 427   0.85986585 0.984214339
 428 -12.77906359 0.882390290
 429   7.21420144 0.550884826
 430   2.31817022 0.231021556
 431  11.60161950 0.888496654
 432  -0.19346228 0.242609713
 433   5.07478120 0.759161318
 434  14.54155003 0.040387654
 435   3.81039636 0.874572741
 436   2.23233049 0.448317248
 437   0.19481869 0.201906051
 438   2.81530451 0.132131690
 439  12.39893259 0.674693704
 440   0.47054642 0.632959494
 441   2.16152913 0.734480632
 442   0.33398836 0.315024718
 443   7.35509037 0.304570986
 444  -2.92336559 0.539062343
 445   5.79622573 0.392393310
 446  -2.37607425 0.403380474
 447   0.04498550 0.756875541
 448  -1.63674414 0.613789514
 449  11.80310547 0.832651469
 450   6.30630243 0.850689403
 451   1.48394652 0.096243229
 452   4.03361865 0.799660045
 453   3.54707273 0.408520520
 454   2.00327040 0.702944912
 455  17.30761707 0.380542812
 456   5.72738968 0.105447516
 457 -13.64604891 0.328506659
 458   8.35976334 0.702173924
 459  -7.41197443 0.134396488
 460 -15.95683040 0.618526462
 461   8.76889573 0.950243069
 462  -1.13482624 0.113477080
 463  -0.60311407 0.090444247
 464   4.95508365 0.612511543
 465   5.36934491 0.979213258
 466  -0.03554882 0.807185690
 467 -11.58131144 0.183341373
 468   4.46809041 0.796330582
 469  12.49741067 0.346860912
 470   8.63824488 0.073684997
 471   0.49990913 0.732519306
 472  12.82688360 0.109400213
 473  13.20375065 0.850369092
 474  -8.41110869 0.177717087
 475  16.31959963 0.727704840
 476  17.59203613 0.235311681
 477   0.32148420 0.842195936
 478   5.43148331 0.670904647
 479   7.14649727 0.028190029
 480   0.25410683 0.421535783
 481 -12.41047826 0.086404379
 482 -10.64180909 0.229659236
 483  -6.40185653 0.876365242
 484  15.63063324 0.667672536
 485   1.94280423 0.799266628
 486  -5.76507450 0.367344192
 487   8.60895533 0.154109357
 488   9.38306751 0.788742770
 489   3.43573528 0.284535277
 490   4.81848966 0.872283177
 491  11.65839314 0.234109111
 492  -5.57884822 0.030363060
 493  -3.94238060 0.325320686
 494   9.38133340 0.201141788
 495  -7.65003459 0.647734396
 496  11.23091019 0.084927159
 497  -6.07705432 0.037273791
 498   7.46380750 0.506897136
 499   7.42034855 0.869351148
 500  -4.43031973 0.231191152
 501  -1.07351537 0.480234836
 502  -1.40653281 0.690620421
 503  -3.82710168 0.990191328
 504   5.04583490 0.543427375
 505 -11.54265099 0.270542185
 506   0.49059479 0.991447248
 507  -1.40871469 0.555998766
 508   3.64241437 0.743840673
 509 -18.30031589 0.357478210
 510   4.27487959 0.770619738
 511   1.28805821 0.654787106
 512  -3.19542768 0.218110139
 513  12.53375654 0.011857644
 514  11.78889419 0.054127726
 515  -5.38392310 0.839309080
 516  16.38024181 0.228801038
 517  -0.59622631 0.134381782
 518  -0.74107258 0.258146632
 519 -12.31429450 0.020524447
 520  -0.79785028 0.968028764
 521   6.39899711 0.038162566
 522   7.42024044 0.716163692
 523  -3.62470664 0.018201813
 524  -2.55049724 0.162446610
 525 -10.79888854 0.683070478
 526  10.18490144 0.546461234
 527  -2.76979044 0.198830067
 528   4.85164813 0.094100357
 529   0.96477200 0.381801756
 530   8.13344336 0.639730450
 531   9.04684412 0.786084368
 532  10.41746272 0.828304181
 533   0.94334368 0.798419831
 534  10.13116556 0.191715972
 535  -4.12728628 0.575178239
 536  -9.59222379 0.876405375
 537   1.64680258 0.391003085
 538  -4.58897613 0.039176486
 539   0.38394379 0.511577564
 540  -4.80428215 0.222785463
 541   0.35363661 0.681658725
 542  -9.63685708 0.183035382
 543   3.54363414 0.766127414
 544   6.89610808 0.967514568
 545  -2.03781105 0.464416752
 546   8.67956196 0.421424078
 547  -1.09959038 0.061231448
 548   7.12587456 0.028601318
 549  -6.93064672 0.402561175
 550   8.57989199 0.925089270
 551  -9.55071810 0.454993099
 552  -8.11914736 0.509644286
 553  -5.41909698 0.077813151
 554 -17.03336572 0.875713545
 555  -1.27438609 0.602163625
 556   3.09834374 0.105599007
 557  -1.59865741 0.439939102
 558  11.82272089 0.754984309
 559   4.30969696 0.483834579
 560 -10.76886192 0.222486992
 561   7.05419803 0.903020271
 562   7.36096847 0.440357053
 563  -2.05864869 0.581170147
 564  -9.08366913 0.318677911
 565   8.57119930 0.605668919
 566   7.87702340 0.570206991
 567   5.22035786 0.542344385
 568   2.37238850 0.595969470
 569  -4.29809941 0.634313781
 570   4.51647479 0.796663089
 571  -0.62478780 0.562099444
 572   8.50866078 0.490014249
 573   3.46694991 0.122890089
 574  -7.31956453 0.885170890
 575   2.20259268 0.167180856
 576  -1.81003626 0.702563515
 577   8.44526939 0.973495019
 578   8.19767069 0.881261264
 579  -5.92422578 0.686557351
 580  -0.11826129 0.712798344
 581   5.66132869 0.922826429
 582  -5.40845018 0.642183516
 583   6.67839036 0.680978989
 584  11.88962825 0.487904896
 585   3.32266332 0.931709581
 586   0.24234019 0.405641313
 587 -12.79023339 0.361005489
 588 -13.57875491 0.266289733
 589   1.81304596 0.775093821
 590   0.36755600 0.400225605
 591  -9.15574205 0.518040748
 592  -3.90436548 0.396869908
 593   9.24764042 0.669374848
 594   0.74869385 0.609881390
 595  -3.62958907 0.928867495
 596  -0.02527232 0.557679930
 597   0.04433418 0.152565816
 598  11.76152632 0.865663501
 599  -4.62181124 0.007000650
 600   5.82271403 0.389678502
 601   0.33289002 0.532940826
 602  -7.65647076 0.681574524
 603  11.81023732 0.107165912
 604  11.42121999 0.989580324
 605  -5.47120641 0.762285550
 606   3.82311561 0.388755074
 607  16.91059711 0.461236022
 608   4.14012105 0.802420151
 609  -1.35278659 0.036646959
 610 -12.81733350 0.179096148
 611  -0.94395134 0.450959878
 612  -5.39002376 0.264783062
 613   4.16227017 0.780743762
 614   7.26179625 0.821382454
 615  15.10062276 0.469253936
 616   1.45877225 0.685434405
 617  -9.87966760 0.767201511
 618   7.03156071 0.195142483
 619 -11.71327419 0.774014869
 620  -4.55518706 0.973103604
 621  -1.75221406 0.175172193
 622  10.35631400 0.080670414
 623   4.97650495 0.597478189
 624   2.25703939 0.585949751
 625  10.72500409 0.339720931
 626  -5.02901029 0.997874377
 627   6.24398637 0.655067479
 628  -5.83880059 0.184948259
 629   2.17256077 0.746741866
 630  -5.59809380 0.277523381
 631   8.19384177 0.334565607
 632   3.35250431 0.952057263
 633  16.20874892 0.901400446
 634   1.63205839 0.235388475
 635  -1.07921163 0.608376332
 636   0.24315118 0.862639830
 637  15.61923078 0.050955422
 638   1.99639207 0.358905687
 639   8.14825538 0.190069662
 640   4.55210835 0.784025901
 641  13.51582298 0.973572910
 642  15.42415796 0.969080992
 643   2.23978124 0.551857514
 644   1.00858991 0.919566804
 645  -2.77293574 0.906998180
 646   7.10750420 0.934792213
 647  -8.01377290 0.682306063
 648   9.67873875 0.239576806
 649   7.54867950 0.065860266
 650  13.85701962 0.733823443
 651   8.48212853 0.285731085
 652   3.55278843 0.998255904
 653  21.94592206 0.205463912
 654  -2.07957143 0.948665109
 655   1.54169997 0.200109744
 656 -11.36934275 0.447122472
 657   3.07094572 0.815147945
 658   6.45818709 0.007849948
 659   1.85594578 0.818796540
 660  -2.43799564 0.962013689
 661 -17.96539549 0.654190963
 662  -0.93433746 0.454930236
 663 -11.06904368 0.898560975
 664  14.89733742 0.479152492
 665  -5.72390675 0.136197255
 666   9.46781102 0.669006610
 667   5.35954546 0.259381138
 668   3.78388994 0.933778797
 669   1.95373423 0.517555994
 670  10.96772341 0.666138826
 671   9.40585102 0.779906833
 672   0.75347502 0.142656741
 673   7.64803672 0.734297119
 674  -0.40051164 0.362230232
 675  10.00747057 0.660820381
 676 -12.86024975 0.072988046
 677   1.43515528 0.229672223
 678  -6.75981709 0.658534537
 679  -5.61355474 0.795897133
 680  -4.40596595 0.038787666
 681  -1.37033650 0.371835229
 682   6.66666573 0.560963737
 683   8.18430044 0.284787698
 684  -0.55742330 0.622783662
 685  -0.39757686 0.673551753
 686 -12.68628005 0.373038561
 687   4.06416215 0.760546238
 688   4.65859855 0.516761886
 689   3.55304076 0.266856843
 690  -7.35294817 0.615783196
 691   1.01222898 0.158266779
 692   9.91052610 0.285619547
 693  -6.42966726 0.573689954
 694  10.97425098 0.985095061
 695   5.79394599 0.404333309
 696  10.09106608 0.441037857
 697  -1.47295537 0.577661077
 698  -2.07959719 0.547176133
 699  -8.76910940 0.498979558
 700 -11.15658312 0.135862745
 701  -0.88456783 0.326480064
 702   9.71607440 0.998076370
 703  -8.76072622 0.386244511
 704 -19.26823092 0.461833959
 705  -0.11280313 0.064155908
 706   0.64625887 0.172078148
 707  -5.35323428 0.331153163
 708  -1.71034509 0.330955888
 709   4.27104744 0.590544244
 710   7.33843789 0.263171531
 711  -5.38121637 0.539675802
 712  -6.87566548 0.127313096
 713  -2.50161298 0.269417630
 714  10.99076986 0.097362729
 715   6.34017269 0.318528587
 716  -4.63672382 0.451038055
 717 -11.55122495 0.987073278
 718   4.78618612 0.297342215
 719   2.97547390 0.197312152
 720  -5.54495280 0.499701114
 721  17.67606173 0.810316588
 722  16.01578815 0.643667608
 723  -0.16258467 0.228284761
 724   7.92123340 0.784289369
 725  -2.26303900 0.270764770
 726   5.84136933 0.437763291
 727   4.96955217 0.389720490
 728  -8.09516710 0.829068548
 729  14.59207207 0.513593803
 730   2.80954688 0.650799867
 731   4.53653552 0.672326278
 732   4.49116737 0.807447691
 733  18.87549709 0.647303378
 734   9.80118464 0.932576117
 735 -13.02124969 0.038651904
 736  -6.15189291 0.697593318
 737  15.81920283 0.249825051
 738  10.81503188 0.152372300
 739 -23.58738366 0.593560367
 740   8.15716338 0.411680007
 741   3.45349379 0.351061414
 742  -6.39345334 0.374926213
 743   8.72924585 0.165759028
 744  22.17948804 0.003736780
 745  -4.73053410 0.582425257
 746  16.88289626 0.484899167
 747  -1.78826142 0.663273340
 748  -0.78106025 0.337039969
 749  -2.92461669 0.810174719
 750 -13.89224399 0.177428986
 751   4.56809819 0.025010350
 752  -1.07452825 0.649632933
 753   0.58148751 0.829606422
 754  12.13329525 0.354819526
 755  17.35605568 0.284862590
 756 -12.43678107 0.827661083
 757  -1.89492796 0.574929572
 758   2.18520382 0.846299917
 759  18.11449649 0.603173531
 760   4.34508582 0.484049042
 761  17.49394569 0.094811656
 762  10.67752350 0.166176400
 763  17.13374502 0.547208197
 764   4.42138123 0.768691494
 765   5.38445574 0.788597361
 766   0.79946671 0.851883720
 767  -4.67547904 0.995621191
 768  -5.61496422 0.523793593
 769 -20.52093126 0.881207308
 770  -8.95996814 0.851078124
 771  -7.63483710 0.739657373
 772  11.02131097 0.678060014
 773 -10.56228517 0.202393048
 774   6.48841788 0.143946271
 775   3.44853632 0.913249620
 776  -0.02080024 0.070765134
 777   2.08654297 0.032468089
 778   8.13415912 0.439470874
 779  11.19028936 0.944954026
 780   0.26670866 0.492724593
 781  -9.33692734 0.982611921
 782  17.23967092 0.313428994
 783   0.36906670 0.660669528
 784   7.89735684 0.977628886
 785  -4.00171487 0.379327632
 786   5.01615432 0.735627296
 787   0.42214214 0.092461754
 788  13.60634772 0.218359635
 789   6.57431413 0.067653525
 790  -1.77668341 0.717799276
 791   5.16227422 0.325502093
 792 -15.29091550 0.332815338
 793   3.33602480 0.594168551
 794  13.80131443 0.817724470
 795   5.92111679 0.947854666
 796   3.59747624 0.330860216
 797  -6.79722403 0.093518715
 798  -1.86606213 0.824179728
 799  17.05226458 0.466573672
 800  10.39712467 0.409067778
 801  -4.78536386 0.891470739
 802  11.92963128 0.719633060
 803  -1.44230992 0.232628002
 804 -12.31860616 0.834134222
 805   2.93439660 0.957842480
 806  14.27963295 0.546264646
 807   2.17488820 0.701170328
 808  10.78772417 0.612332448
 809  -0.47049341 0.378564293
 810  -0.35140634 0.034337429
 811   5.04887868 0.211697132
 812  -3.51562580 0.663243607
 813  -0.82013387 0.602497174
 814   2.78954743 0.325294790
 815   8.67905777 0.820296625
 816 -12.70343389 0.315467361
 817  -2.59373236 0.015571904
 818  -4.60369241 0.293737716
 819   1.58669084 0.671091860
 820 -10.44245103 0.501340276
 821   4.85215578 0.141572007
 822  10.46303284 0.801814632
 823   1.27898298 0.236929983
 824  -1.72225479 0.608500539
 825  20.18685735 0.827124630
 826   3.27308817 0.542065179
 827   1.01596956 0.254672115
 828  -8.88872881 0.460876757
 829 -11.31397349 0.636168639
 830   0.85294367 0.816417328
 831   3.54262337 0.944147626
 832 -10.53603202 0.675775741
 833   4.34832198 0.121988381
 834  11.56451662 0.283063133
 835  -7.36454369 0.500596540
 836  -8.23701113 0.379483261
 837  -8.36081323 0.219730782
 838  -6.39158860 0.739171315
 839  -1.40518544 0.478709398
 840  -4.01314821 0.460476388
 841  -7.34814047 0.406242873
 842  -7.80836711 0.730648091
 843  -0.57729135 0.152336258
 844   4.98352832 0.026424939
 845  -3.78181635 0.453598432
 846  20.16821827 0.845273124
 847   5.20758271 0.573569671
 848  -3.05534245 0.286828574
 849  -5.31306254 0.961990401
 850   1.09307567 0.006478724
 851  -3.75412572 0.598277695
 852  -2.38444245 0.777900122
 853   2.46837742 0.280363751
 854   9.72195519 0.041094463
 855   3.96271247 0.604775284
 856   2.14105354 0.400315328
 857   7.88645912 0.404573389
 858  -4.03565076 0.798377309
 859  10.80180959 0.932152434
 860 -10.89359212 0.446813857
 861   1.43144578 0.310194540
 862   4.79825196 0.504826858
 863  10.73201365 0.384306369
 864  -4.07526187 0.893893643
 865  -2.84330198 0.390202663
 866   5.81825057 0.830581384
 867  -2.77842745 0.382966910
 868  -7.70333673 0.157692966
 869  -3.25753058 0.726303603
 870   8.50032387 0.556524444
 871   2.35027236 0.857076526
 872  -1.70740565 0.194760923
 873  -3.40693880 0.696420946
 874  -8.03983352 0.514393263
 875  -1.85105344 0.609459979
 876  -9.01148029 0.526019631
 877  18.37344635 0.690793045
 878  16.46079416 0.811535334
 879   4.10224315 0.043403618
 880   7.06657672 0.831274577
 881  15.31421824 0.434558881
 882 -12.36760970 0.004215634
 883   1.95473415 0.277788662
 884  -0.93207006 0.368433415
 885  15.39919341 0.843189783
 886   5.23452387 0.626226925
 887  11.40805770 0.002417288
 888  -1.30282837 0.072493756
 889   3.92130690 0.675355182
 890   2.53148399 0.027222295
 891   4.92705318 0.934429364
 892   5.54978818 0.042268708
 893  -2.19608977 0.246743834
 894  -0.62565550 0.858214200
 895  -8.98329365 0.646827226
 896  12.78468146 0.533966352
 897   2.01061290 0.418710227
 898   1.03689579 0.019241741
 899   8.01166696 0.992268130
 900  -4.49786437 0.694127903
 901  -8.15387184 0.066275002
 902   2.22256207 0.083301613
 903 -12.27145086 0.535369809
 904   9.95709112 0.227692557
 905  14.58198717 0.667298058
 906   5.98046083 0.922503625
 907   1.25640725 0.632933575
 908   9.77623752 0.136171032
 909   5.57068426 0.374916651
 910 -10.07048336 0.470411379
 911   3.69267954 0.897278365
 912   2.22185354 0.212539549
 913   7.96155623 0.720525208
 914  -6.21272358 0.771491819
 915   2.63054735 0.474989115
 916  -2.81488890 0.675381020
 917   4.52747191 0.118615879
 918  -3.22975936 0.783991133
 919  11.42834761 0.423344824
 920   0.26512464 0.617515445
 921  -5.84322807 0.210915613
 922   9.61073028 0.988117333
 923  -6.11878012 0.492318959
 924   5.30581443 0.339379499
 925  -6.40132703 0.903540026
 926   1.22921808 0.122161655
 927   8.08547837 0.197296588
 928  -0.77943801 0.935963718
 929  11.43194858 0.828270943
 930  -5.41689395 0.556863468
 931  15.14667847 0.565186375
 932  -5.15327419 0.542802437
 933  -3.95903082 0.643379366
 934   5.78847793 0.391369361
 935  11.54430873 0.158789330
 936   1.90340148 0.841316129
 937  14.69680285 0.532022770
 938   0.68552840 0.367073827
 939  -8.72287967 0.250127491
 940   9.35401445 0.836083158
 941   5.32139524 0.996712598
 942 -14.53387897 0.825434481
 943  -2.93925146 0.513153861
 944  12.54386493 0.713306793
 945   2.04842442 0.993893406
 946   2.87461954 0.049843312
 947   4.89765230 0.376710062
 948  -6.23945314 0.321108142
 949  -3.45840168 0.854710947
 950   9.05807160 0.199992188
 951   3.33815006 0.787302467
 952   4.22244242 0.351841910
 953  15.75879160 0.268699469
 954   2.78549859 0.920299974
 955  -4.46643118 0.727283862
 956   0.48021298 0.428672083
 957   2.55814938 0.130915212
 958   5.00692968 0.062266047
 959   4.78801127 0.325124688
 960   6.39524485 0.693406744
 961 -10.46792584 0.458128441
 962  10.14111908 0.353412759
 963 -10.56424183 0.821588957
 964   7.60967746 0.267669137
 965  -2.34956688 0.434855697
 966  23.82269027 0.802311880
 967   8.37170447 0.445185000
 968  10.05024769 0.778687843
 969  -9.15753018 0.957292819
 970  12.17438228 0.774769426
 971   1.57960028 0.783591989
 972   0.06719501 0.849073924
 973  16.21114558 0.243444943
 974  -3.79808298 0.842994720
 975   8.98927715 0.020537113
 976   7.72362992 0.984168340
 977  11.25158442 0.152385348
 978 -21.23936903 0.909204114
 979   7.34995949 0.987249305
 980  -7.99435203 0.335456401
 981  -2.78218185 0.768517548
 982  11.59547596 0.466617637
 983  15.90870706 0.071892573
 984   5.58160897 0.554485703
 985  16.05253351 0.815206562
 986  -3.23103465 0.280495460
 987  -4.61108636 0.035757819
 988   5.41596511 0.746146856
 989   2.92445613 0.136743821
 990  11.23628254 0.681316365
 991 -12.93714705 0.838791576
 992   9.94668264 0.084457395
 993  -4.56061529 0.983605894
 994  -4.24795688 0.601732731
 995  -2.83740044 0.375102341
 996  -0.43078317 0.403870303
 997  15.19689584 0.114826374
 998 -10.29920266 0.731582141
 999   6.01686515 0.641655876
1000   6.69431335 0.496723697
1001   4.62223602 0.328118236
1002  -1.74309026 0.072604771
1003  14.31971261 0.827101483
1004  -1.86629155 0.613346722
1005   8.30971428 0.274948560
1006   8.50080711 0.059822908
1007  -7.94061422 0.121069240
1008  -2.72096492 0.710791774
1009   3.33259421 0.398621625
1010   1.73248470 0.488581205
1011   9.56008489 0.011104565
1012  12.71499762 0.038568985
1013   4.11512127 0.219846314
1014  -0.96707584 0.822646857
1015   4.98621667 0.633779997
1016   4.69384821 0.295708955
1017  10.16008645 0.778287787
1018  -7.72973800 0.097096969
1019   2.87264210 0.796538177
1020   4.56095440 0.862952770
1021   5.02621658 0.934628629
1022   3.18138681 0.805600816
1023  -1.02245780 0.317640678
1024  18.16001652 0.992503640
1025   4.13729026 0.941910149
1026   1.61211303 0.377271914
1027   1.71520009 0.735196094
1028   3.26325421 0.514432564
1029  12.94663819 0.591190711
1030  10.53239931 0.005877708
1031   8.06705056 0.340779884
1032  -5.09007267 0.332516161
1033  12.31973355 0.323119296
1034  -2.69957650 0.633232996
1035  12.51207803 0.377641090
1036   8.02081444 0.859293157
1037  -0.13098726 0.099370804
1038  -0.97757546 0.852873609
1039  16.73605399 0.595854575
1040   3.63219184 0.329310613
1041  -4.79105630 0.247760146
1042  -4.77209495 0.708235587
1043   0.92107647 0.924567254
1044  12.12724271 0.433550712
1045  -5.07731478 0.200109463
1046   9.16019579 0.897456586
1047  18.33260560 0.649877409
1048   1.93596773 0.584401505
1049   8.51254631 0.283154523
1050  11.41092928 0.698703314
1051  10.85035748 0.351078210
1052  12.62749979 0.570101319
1053  -2.32028296 0.313842122
1054  -2.45778301 0.007943144
1055   6.93102526 0.108737491
1056  -0.67304654 0.245399613
1057   9.27294774 0.204010286
1058  14.29292826 0.396294626
1059  14.05843185 0.864613328
1060  -3.73515954 0.305862948
1061   0.36606339 0.116802407
1062  -5.79235478 0.457308058
1063   8.70346900 0.858244380
1064  -8.91321043 0.077001581
1065   0.58499566 0.503209780
1066   0.39160153 0.324883353
1067   7.46715326 0.343451039
1068  12.36256009 0.679483638
1069   8.84283689 0.687359177
1070  -6.39396909 0.113065562
1071  -3.67844896 0.667335667
1072   9.36904962 0.009815419
1073  -3.25244888 0.213105120
1074  19.09389976 0.593130536
1075   7.28826611 0.829483570
1076  -5.44565944 0.956490203
1077   7.96993416 0.770961635
1078   0.20683778 0.006497153
1079  -3.73273760 0.037042812
1080 -10.64745846 0.813594448
1081   5.70578906 0.157678242
1082   4.05282218 0.224663656
1083  14.77711159 0.577586777
1084   0.89685942 0.297213941
1085   3.92600687 0.672347849
1086 -12.29347477 0.367072171
1087  -9.33603480 0.456544225
1088  -0.86683190 0.088696811
1089   4.65685745 0.779783359
1090   1.24438030 0.712958633
1091  11.43533814 0.920345548
1092 -10.18380242 0.044456697
1093  -1.20684029 0.992051648
1094  -9.78059038 0.611477837
1095   3.05588762 0.581933667
1096   3.47419279 0.769325101
1097   0.87528245 0.455214184
1098  -3.13185655 0.805887381
1099  -0.82283965 0.707668384
1100  -1.86717272 0.984060013
1101  16.56357048 0.217369677
1102  -2.11052646 0.474156371
1103  -1.39795364 0.958554209
1104   4.87468692 0.328779186
1105   2.69163553 0.401633221
1106   6.08640626 0.599963560
1107   7.41420081 0.240202007
1108   5.73729928 0.696034193
1109   7.31747120 0.569520861
1110  -6.20465547 0.214005920
1111  17.52477873 0.667125450
1112  12.97855692 0.796977778
1113  -3.35883428 0.379721403
1114  -2.90306972 0.552454626
1115   5.31617371 0.401625473
1116  -3.86414389 0.830986352
1117 -14.94107832 0.702705123
1118  -5.74060402 0.833328045
1119  -8.10116203 0.078855027
1120  23.48247017 0.568666620
1121  20.22005082 0.357069809
1122  -2.53387193 0.637455425
1123  15.72048831 0.845354124
1124  -4.41494567 0.934471473
1125  -8.02254420 0.378467959
1126  -0.13398716 0.489382793
1127   0.95967155 0.813667919
1128   0.14835664 0.215786848
1129  14.31875579 0.675145039
1130  -6.36589196 0.822037848
1131   8.25942906 0.156787526
1132  -7.33597529 0.051076292
1133  12.58936771 0.666507807
1134   2.34653798 0.626196518
1135  -0.69351398 0.050664564
1136   7.08738260 0.808776877
1137   5.19653521 0.779008623
1138   3.20900427 0.197212774
1139   7.81171331 0.744975548
1140   6.49008186 0.991318119
1141   7.27471854 0.839642650
1142  -7.68367290 0.880500743
1143  12.04846713 0.797754890
1144  14.93435279 0.190527791
1145  -3.83641079 0.075995951
1146   2.15090497 0.426560973
1147  -3.61166623 0.777188818
1148   8.49333248 0.891445999
1149  -7.46936100 0.148607446
1150  13.85406193 0.983656455
1151  12.20477754 0.499345090
1152  10.09415710 0.638127733
1153   5.37134772 0.110929011
1154   8.17660840 0.879411588
1155  -4.38804367 0.608933700
1156 -11.78145902 0.265134740
1157   6.18940186 0.970982743
1158 -16.24831477 0.844983635
1159   9.52790402 0.578152651
1160  16.44372225 0.264144422
1161  -2.48286428 0.893865621
1162   5.33297280 0.512990215
1163  -2.68912507 0.851636020
1164   9.94607707 0.644483197
1165  -1.93526852 0.550759844
1166   2.34310539 0.787853650
1167  11.79131608 0.983668283
1168   3.16689104 0.605394987
1169  10.47759320 0.919442774
1170   2.86973133 0.557835916
1171  10.30674302 0.442504870
1172  10.92820575 0.976183635
1173   7.98050212 0.139334994
1174  -0.64719705 0.981199028
1175  -2.63625596 0.341524563
1176  11.38799583 0.858987904
1177   1.37321916 0.202373294
1178  12.66698520 0.142127091
1179   5.83599540 0.864497670
1180   4.88659560 0.472598564
1181  13.00108599 0.961629827
1182   5.79514791 0.408377170
1183  -1.47807631 0.536772872
1184  -3.38142805 0.288956265
1185   6.25154986 0.828695103
1186   2.40919373 0.478123848
1187   3.72990486 0.056539500
1188  -9.90915815 0.603356617
1189   0.21737084 0.737251896
1190   5.36929388 0.026920178
1191  -1.05027354 0.034992509
1192  -4.97887624 0.506301429
1193  -6.40058435 0.014061876
1194  -0.14610837 0.619699963
1195   3.78483619 0.653952701
1196   3.84143365 0.162122572
1197   2.66030676 0.196542503
1198 -10.56809462 0.386200215
1199  -5.01140125 0.711703654
1200  -3.09809005 0.118120179
1201  -2.76110171 0.118809515
1202   2.85825107 0.129646974
1203   2.75993661 0.171779333
1204  11.55931169 0.372165133
1205   9.21211486 0.969079819
1206   6.02207148 0.498965865
1207  -3.52883224 0.954619249
1208  -2.60190803 0.069405278
1209   1.34183694 0.569402487
1210 -11.35155228 0.766344735
1211   1.04661568 0.023673810
1212  -1.90461932 0.179728300
1213  13.72465582 0.467775796
1214  19.14882438 0.476924297
1215  -1.07480326 0.944407858
1216  -8.44289331 0.059804028
1217   1.89732882 0.743225795
1218  -7.87832463 0.672539050
1219 -12.24163608 0.916803014
1220 -12.77212790 0.648129714
1221   6.39197262 0.622954436
1222   5.26261666 0.494421400
1223 -10.65239640 0.695527931
1224   4.63841458 0.499519163
1225  -2.94276544 0.429201572
1226   4.68788953 0.639613685
1227  -1.03031400 0.349342009
1228  -2.69946354 0.221796918
1229 -15.32237714 0.631289988
1230  -8.31962698 0.925363812
1231  -5.80897714 0.833536878
1232   7.16070989 0.832098478
1233 -10.99679727 0.794048223
1234   0.84514458 0.748014415
1235   2.23308495 0.111176288
1236   3.56351018 0.599805508
1237   0.88336430 0.746908710
1238 -14.63461670 0.314391808
1239   4.39039715 0.079604833
1240  -7.07001439 0.633705345
1241   2.11252583 0.461468123
1242   7.60219364 0.497389476
1243  -4.87713428 0.039952736
1244   2.17515292 0.421830084
1245   0.64302362 0.267982804
1246 -22.29371533 0.646257366
1247   0.31652779 0.060548371
1248   7.93445046 0.343570449
1249   0.28292029 0.651909785
1250   2.77775640 0.637679287
1251   6.22941586 0.291132945
1252  23.68567532 0.708513840
1253   9.49503014 0.645200206
1254   0.87405420 0.063154289
1255  -4.04931224 0.110797498
1256   8.91607239 0.732917195
1257  -5.77728018 0.635435595
1258 -16.37296319 0.343727613
1259   9.87409940 0.774177478
1260  -8.11360210 0.377765616
1261  14.54242540 0.204343527
1262   0.36239636 0.115528352
1263  19.51009176 0.181365423
1264   1.23592729 0.011676577
1265 -15.81877035 0.767155028
1266  -0.05911251 0.737944231
1267  -6.55395965 0.214062137
1268  -7.85591487 0.539865054
1269  -9.73010882 0.730924287
1270  11.79433862 0.267116856
1271  -8.84308360 0.088069165
1272  -5.56689174 0.405987947
1273   7.59010135 0.655631611
1274  10.07629305 0.031106157
1275  -6.19331485 0.052350502
1276  -4.58626710 0.326901540
1277  -5.19431549 0.125740555
1278  -2.08129025 0.034657174
1279 -10.48798034 0.153632237
1280  13.04657686 0.317295703
1281   1.94142067 0.731437668
1282  -1.62470735 0.701070475
1283 -12.27046912 0.505781742
1284  -2.96095228 0.122808075
1285  -2.91847765 0.372668438
1286 -14.83230131 0.100749725
1287  16.57350659 0.707854947
1288  10.05473238 0.244046174
1289   5.50858969 0.070691273
1290   7.65309196 0.245393047
1291   7.16359996 0.056261015
1292   4.33026356 0.318855549
1293  -4.65721575 0.271249938
1294   2.85909691 0.309377566
1295   3.02736080 0.553944209
1296   6.22796768 0.763945813
1297  -4.47036396 0.197721195
1298   2.78901176 0.441166128
1299  -9.94574794 0.964660659
1300   1.86451969 0.704635530
1301 -10.38926659 0.772304221
1302  -5.36565800 0.029527218
1303   1.99230152 0.578448308
1304  13.65547415 0.936050102
1305  -2.05229879 0.851521142
1306   0.99504588 0.974891334
1307  -1.46027404 0.320227281
1308  -8.45614275 0.727910071
1309  -0.95201934 0.199101032
1310  -2.46642929 0.462252060
1311  -6.44060430 0.703637604
1312  -2.58115910 0.084948525
1313   0.76248197 0.125769097
1314  12.00603845 0.675927328
1315   1.97538215 0.782502470
1316   2.23331320 0.870228155
1317  -3.10226060 0.485056198
1318  12.59337170 0.584729095
1319  -2.42247402 0.387588168
1320   9.41981063 0.374604221
1321   6.26806243 0.727453335
1322  -5.30630356 0.427294265
1323  13.81542647 0.394246994
1324   1.05647858 0.646684666
1325 -12.25005208 0.010531726
1326  -4.58162076 0.077133994
1327   0.58094190 0.400275636
1328   5.79443858 0.641731247
1329   8.87635216 0.913593476
1330   9.71048520 0.955285711
1331  13.10563373 0.908471848
1332   4.99194220 0.967014095
1333  15.88178853 0.041518216
1334   9.35962068 0.864770023
1335  -7.53095731 0.300106124
1336  12.18427585 0.248876997
1337   9.22034502 0.450149366
1338  -1.02861237 0.684246939
1339  -2.98140404 0.326901490
1340  -4.64316598 0.425381055
1341  15.35233259 0.630774937
1342  -1.85655250 0.226889991
1343  15.43748330 0.584219351
1344  10.39060893 0.387854461
1345   2.80705696 0.564024865
1346   3.48201221 0.787103673
1347   7.03787977 0.112019552
1348   8.41853061 0.798376796
1349  15.63925527 0.873984550
1350  -4.05742183 0.699131238
1351   6.56954685 0.720018710
1352   2.44007265 0.232697343
1353   3.75597926 0.975133449
1354   2.92362149 0.290975435
1355  -4.74372257 0.003738451
1356   1.28365940 0.987536495
1357  15.65288265 0.179629701
1358 -11.76385004 0.850614822
1359   1.56331228 0.592017435
1360  -9.64774741 0.024951969
1361  -9.44879860 0.993960270
1362  29.33340056 0.913358233
1363   7.97233120 0.021820585
1364 -12.10837953 0.401535846
1365  -1.20729618 0.984977268
1366   3.63219301 0.491142613
1367   2.79853507 0.663823888
1368   3.19584583 0.612511282
1369  -0.81790885 0.908769330
1370  -1.67795944 0.611690031
1371   2.55137163 0.109447998
1372   4.36572889 0.382049700
1373  -6.35667866 0.162787163
1374  -0.76239101 0.892383562
1375  -3.99558996 0.466572017
1376  -0.47513018 0.457760464
1377 -10.69568261 0.544872910
1378   4.30943512 0.982456072
1379   2.91825703 0.823403368
1380   0.10753188 0.945676881
1381  -8.38623073 0.923085521
1382   4.95690232 0.188128654
1383   5.39956649 0.331692462
1384  -1.47421789 0.327711090
1385  -1.81689665 0.713285385
1386   5.15137860 0.414906436
1387   5.68897151 0.110799415
1388   0.78825159 0.396824099
1389  -1.78376652 0.929264595
1390   0.76991060 0.950124414
1391  15.81469073 0.951245195
1392  -4.33820920 0.009896093
1393   1.67174323 0.821983745
1394   0.38997945 0.928857784
1395   1.97848484 0.175680230
1396  -5.81067801 0.772580245
1397 -10.45208478 0.418845035
1398  13.34024524 0.905645046
1399  -8.79585122 0.906516178
1400   2.89093397 0.113010960
1401   2.22324289 0.799940482
1402   8.95497981 0.984663669
1403   0.93288527 0.277914575
1404 -17.35306978 0.455587022
1405  -3.26914604 0.406757639
1406   8.75871227 0.067059659
1407   1.79914932 0.784879863
1408  -0.67305388 0.006393497
1409   1.66805704 0.039614073
1410   9.03868439 0.601066847
1411   4.29458670 0.015772820
1412  -8.15564320 0.939633197
1413  12.50538902 0.766347628
1414  -0.45547258 0.464314122
1415  -9.47180656 0.640114882
1416 -13.25567198 0.125841930
1417   2.87660101 0.381931128
1418   7.37834152 0.648958712
1419  -0.45874073 0.303139498
1420   4.87941450 0.500090729
1421   4.50344891 0.311329309
1422  -6.14257896 0.061368838
1423   5.98243271 0.873804882
1424  -2.64694079 0.080493398
1425 -17.79727796 0.188420116
1426 -13.52552336 0.798403568
1427   2.29086373 0.518700767
1428   5.21493652 0.788828533
1429  -8.09641615 0.775041349
1430   5.87005782 0.079447757
1431  10.74795720 0.955691540
1432  -8.01115709 0.004508053
1433  -7.53735064 0.054195934
1434  -6.79130165 0.877193354
1435  -1.26419539 0.837772170
1436   8.31082852 0.967509866
1437  21.83090247 0.261529880
1438  11.20453234 0.913858875
1439   7.19128396 0.541942489
1440  -2.93623595 0.860095891
1441  -3.61446403 0.022418065
1442  -6.59997709 0.532998307
1443   3.71486934 0.522669434
1444  18.03420874 0.295064126
1445  -8.75452291 0.390175021
1446  -7.83680812 0.263760724
1447  -1.10263921 0.501819826
1448   2.05633484 0.338684642
1449   5.25636848 0.558667384
1450  -7.33260497 0.457327559
1451   3.86721296 0.612182242
1452 -15.94331373 0.478329365
1453   3.71501899 0.264241588
1454  18.26175822 0.023212661
1455  -5.21093378 0.184378036
1456  -2.44074986 0.297114134
1457 -11.88339919 0.875956945
1458   7.52127093 0.927322099
1459   5.31597834 0.416344968
1460  11.42012314 0.952491078
1461  13.64950955 0.794183413
1462  12.50861255 0.390723282
1463  -3.48142207 0.538708662
1464  14.32910902 0.085221990
1465   5.76196699 0.313860477
1466   2.63751452 0.917424732
1467   2.99975231 0.208662214
1468   7.09852825 0.798246964
1469  -1.95742636 0.352166210
1470   7.80534904 0.386523123
1471  -4.47229047 0.290188493
1472  -4.35535158 0.761527294
1473   1.47083860 0.447897289
1474   3.09504296 0.048513534
1475   3.50446804 0.925072429
1476  12.00487617 0.574499971
1477  10.35171466 0.934193962
1478  -5.63256003 0.968833982
1479   7.15625220 0.467160468
1480  -7.81378393 0.790220187
1481   4.52101003 0.014459969
1482  12.90773453 0.990835752
1483   7.70737851 0.785329264
1484  -3.37196794 0.066025357
1485  -5.12793918 0.347459322
1486  -7.96083724 0.216608294
1487  12.81247279 0.287880308
1488   4.63872463 0.426881173
1489  -3.85439309 0.336532356
1490   4.55633320 0.108001536
1491  -2.40824634 0.135247519
1492   1.65932541 0.005108006
1493   3.26129578 0.093163961
1494  -3.52114597 0.544041275
1495  -9.08479260 0.111212700
1496   7.75150456 0.942726234
1497   7.44829768 0.396996218
1498  14.44430576 0.525470762
1499  -2.13457508 0.207577358
1500   9.74871681 0.537177845
1501  -4.53338693 0.625854028
1502  16.15962824 0.947933141
1503  -0.17711664 0.480940902
1504  -7.21470818 0.006952612
1505  -6.27644212 0.737909602
1506  -0.81648165 0.230003567
1507  -2.10429152 0.209671398
1508   7.69291241 0.987903443
1509  -0.32284504 0.183904658
1510  -0.90833921 0.782169770
1511  10.35542238 0.201758865
1512  10.40788689 0.540802365
1513  10.80011849 0.298263948
1514   7.39943598 0.785716539
1515 -12.71674257 0.154135834
1516  16.67139866 0.116794235
1517  12.47832985 0.998179468
1518   3.24041348 0.653080096
1519  -5.50381593 0.995396942
1520 -10.41952811 0.576472768
1521   4.33514092 0.146434686
1522   4.41294276 0.507165968
1523  -0.14746982 0.698144836
1524  -1.33323051 0.466481571
1525  -7.01201350 0.797150114
1526   6.58669848 0.942287809
1527   7.19444974 0.053569397
1528  -5.66046997 0.435728340
1529  -5.07828702 0.497727572
1530  13.72045272 0.324222944
1531   9.99111984 0.355713969
1532  -4.42363728 0.071790181
1533  -2.34300923 0.618434528
1534   4.98594041 0.605667438
1535  -2.45307721 0.894546647
1536  -3.52276424 0.760086779
1537  -3.69489441 0.960758209
1538  13.04792817 0.511273320
1539  21.61433486 0.236270637
1540  -9.57303968 0.964235539
1541  11.70400744 0.391045695
1542   4.25170422 0.411577090
1543   7.87516537 0.952858161
1544   9.89202673 0.971106687
1545   7.51554467 0.505791978
1546   2.17944879 0.893835908
1547   0.82420351 0.213912155
1548   2.47121932 0.688019842
1549 -14.88503628 0.640950883
1550   8.16032283 0.277742858
1551  -4.65776244 0.129415853
1552   7.48274838 0.074213153
1553   7.70537066 0.476778957
1554   5.88202944 0.351838898
1555   1.95618325 0.106331699
1556   7.22064623 0.511434587
1557 -18.64632081 0.009314188
1558  -6.16794611 0.526204245
1559   2.14042033 0.675800465
1560  -1.89535048 0.916845690
1561  -7.77156605 0.069742819
1562  -6.84078801 0.865082345
1563   8.17539904 0.095895629
1564  10.75170309 0.821383078
1565  14.31498805 0.117893208
1566  -2.82264467 0.809086411
1567   0.11117380 0.400587471
1568   6.43898314 0.333163663
1569   9.48110784 0.465173316
1570   5.39395511 0.695273081
1571  -2.05636570 0.508060862
1572   0.68666117 0.647109494
1573  -6.41880322 0.267530762
1574  -0.12096589 0.986901165
1575   6.46062643 0.588580914
1576  -4.20926136 0.550783675
1577   4.07354300 0.907963701
1578  10.84045143 0.900920521
1579   2.64504664 0.767700269
1580  10.34578229 0.197810342
1581  -0.19222273 0.281932395
1582   3.47400952 0.977555902
1583  11.04549389 0.694010579
1584   6.79729856 0.056652433
1585   9.28300628 0.598930531
1586  -3.53453282 0.183412212
1587   8.04028248 0.250746943
1588   5.75964045 0.424692336
1589  -4.98252741 0.867446071
1590  12.00352175 0.289615080
1591   7.53497791 0.350915526
1592   2.54579776 0.655113837
1593  -9.29572208 0.900136667
1594   6.41659249 0.100570650
1595   7.37095646 0.907179211
1596 -12.78417775 0.262214556
1597   0.87962107 0.624657444
1598  -5.96939907 0.296725805
1599  -2.56857339 0.604065931
1600   4.27131811 0.962952479
1601  21.72603838 0.442485270
1602   6.10056565 0.418383130
1603  10.48099521 0.333593221
1604  19.28363092 0.382408442
1605   2.12080726 0.601206970
1606  -6.82450704 0.740158518
1607  11.32395692 0.627015570
1608   5.00040701 0.476274658
1609 -11.64750733 0.105099095
1610   5.77442654 0.576560214
1611   0.31340364 0.516479036
1612  -2.09881449 0.146089191
1613   5.12411327 0.368130477
1614   1.70530391 0.621828438
1615 -12.95649749 0.355726301
1616   8.43735652 0.275383759
1617 -15.56161079 0.413160084
1618   5.28942694 0.069125495
1619   5.96040877 0.438716686
1620  -2.59318107 0.571116303
1621   6.95988992 0.650760909
1622  14.00074797 0.623645969
1623   1.66101456 0.558763985
1624  -2.57968349 0.648185379
1625  -5.47584253 0.716901151
1626   6.37222581 0.060563130
1627   2.83664864 0.842419730
1628   1.48926558 0.620280308
1629   0.33471689 0.170312461
1630   5.21648412 0.317639631
1631   0.51733642 0.843867329
1632   9.86005834 0.306036746
1633  -5.81145791 0.975655452
1634  -5.43219061 0.303385368
1635   5.87157118 0.677369776
1636   2.08889926 0.310200439
1637  -2.53433085 0.194730908
1638   7.01359575 0.674259533
1639  -2.00936260 0.682056466
1640  -2.98240739 0.787899917
1641  -7.43289210 0.357483044
1642 -12.58905988 0.981387385
1643   5.78095517 0.533526274
1644  -1.23065889 0.687266774
1645  -6.82309960 0.293249774
1646   8.47000829 0.842056399
1647  -5.81624772 0.303700280
1648 -14.83571031 0.311387926
1649   4.66808472 0.091222946
1650  -2.90144463 0.438301785
1651  10.62458662 0.828335698
1652   7.88002491 0.990156110
1653  10.27680283 0.251087079
1654  -9.42498970 0.292462244
1655   6.73027640 0.213065205
1656   1.28169895 0.353152789
1657 -14.29203733 0.264563048
1658  20.35772711 0.265208837
1659   3.55095071 0.242905653
1660 -17.97067670 0.373951756
1661  10.53141139 0.247520698
1662   0.05293205 0.579940423
1663  12.79674707 0.288031751
1664  -5.44235185 0.075899079
1665  14.29464811 0.960707538
1666  -1.36753291 0.124265178
1667  -4.25946974 0.521720352
1668 -12.46519252 0.385503339
1669  -6.65343143 0.540942219
1670   5.55949184 0.143194404
1671  -1.20480594 0.515905644
1672  -4.13839908 0.164461445
1673  -2.21345425 0.812969725
1674   3.94223380 0.229238952
1675 -10.78661097 0.395049514
1676   3.06997341 0.791234255
1677  24.82205477 0.110859039
1678   6.28791249 0.867125744
1679  -2.80296119 0.703583849
1680  13.24274039 0.425951975
1681  -0.19577471 0.361568727
1682  -2.34894781 0.954814545
1683  19.76339577 0.635462177
1684  -1.87591480 0.149121567
1685  -7.70962391 0.711708342
1686  -2.46291902 0.390902746
1687 end data
1688 regression /variables=v0 v1 /statistics defaults /dependent=v0 /method=enter.
1689 ])
1690
1691 AT_CHECK([pspp -O format=csv regression.sps], [0], [dnl
1692 Table: Reading free-form data from INLINE.
1693 Variable,Format
1694 v0,F8.0
1695 v1,F8.0
1696
1697 Table: Model Summary (v0)
1698 ,R,R Square,Adjusted R Square,Std. Error of the Estimate
1699 ,.05,.00,.00,8.11
1700
1701 Table: ANOVA (v0)
1702 ,,Sum of Squares,df,Mean Square,F,Sig.
1703 ,Regression,235.23,1,235.23,3.58,.059
1704 ,Residual,98438.40,1498,65.71,,
1705 ,Total,98673.63,1499,,,
1706
1707 Table: Coefficients (v0)
1708 ,,Unstandardized Coefficients,,Standardized Coefficients,,
1709 ,,B,Std. Error,Beta,t,Sig.
1710 ,(Constant),1.24,.42,.00,2.95,.003
1711 ,v1,1.37,.72,.05,1.89,.059
1712 ])
1713
1714 AT_CLEANUP
1715
1716 AT_SETUP([LINEAR REGRESSION no crash on all missing])
1717 AT_DATA([regcrash.sps], [dnl
1718 data list list /x * y.
1719 begin data.
1720  . .
1721  . .
1722  . .
1723  . .
1724  . .
1725  . .
1726  . .
1727  . .
1728  . .
1729  . .
1730 end data.
1731
1732
1733 regression /variables=x y /dependent=y.
1734 ])
1735
1736 AT_CHECK([pspp -o pspp.csv regcrash.sps], [1], [ignore], [ignore])
1737
1738 AT_CLEANUP
1739
1740
1741
1742 AT_SETUP([LINEAR REGRESSION missing dependent variable])
1743
1744 dnl Test for a bug where missing values in the dependent variable were not being
1745 dnl ignored like they should have been.
1746 AT_DATA([reg-mdv-ref.sps], [dnl
1747 data list notable list / v0 to v2.
1748 begin data
1749  0.65377128  7.735648 -23.97588
1750 -0.13087553  6.142625 -19.63854
1751  0.34880368  7.651430 -25.26557
1752  0.69249021  6.125125 -16.57090
1753 -0.07368178  8.245789 -25.80001
1754 -0.34404919  6.031540 -17.56743
1755  0.75981559  9.832291 -28.35977
1756 -0.46958313  5.343832 -16.79548
1757 -0.06108490  8.838262 -29.25689
1758  0.56154863  6.200189 -18.58219
1759 end data
1760 regression /variables=v0 v1
1761              /statistics defaults
1762              /dependent=v2
1763              /method=enter.
1764 ])
1765
1766 AT_CHECK([pspp -o pspp-ref.csv reg-mdv-ref.sps])
1767
1768 AT_DATA([reg-mdv.sps], [dnl
1769 data list notable list / v0 to v2.
1770 begin data
1771  0.65377128  7.735648 -23.97588
1772 -0.13087553  6.142625 -19.63854
1773  0.34880368  7.651430 -25.26557
1774  0.69249021  6.125125 -16.57090
1775 -0.07368178  8.245789 -25.80001
1776 -0.34404919  6.031540 -17.56743
1777  0.75981559  9.832291 -28.35977
1778 -0.46958313  5.343832 -16.79548
1779 -0.06108490  8.838262 -29.25689
1780  0.56154863  6.200189 -18.58219
1781  0.5         8         9
1782 end data
1783
1784 missing values v2 (9).
1785
1786 regression /variables=v0 v1
1787              /statistics defaults
1788              /dependent=v2
1789              /method=enter.
1790 ])
1791
1792 AT_CHECK([pspp -o pspp.csv reg-mdv.sps])
1793
1794 AT_CHECK([diff pspp.csv pspp-ref.csv])
1795
1796
1797 AT_CLEANUP
1798
1799 AT_SETUP([LINEAR REGRESSION with invalid syntax (and empty dataset)])
1800
1801 AT_DATA([ss.sps], [dnl
1802 data list notable list / v0 to v2.
1803 begin data
1804 end data.
1805
1806 regression /variables=v0 v1
1807              /statistics r coeff anova
1808              /dependent=v2
1809              /method=enter v2.
1810 ])
1811
1812 AT_CHECK([pspp ss.sps], [1], [ignore])
1813
1814 AT_CLEANUP
1815
1816
1817 dnl The following example comes from
1818 dnl http://www.ats.ucla.edu/stat/spss/output/reg_spss%28long%29.htm
1819 AT_SETUP([LINEAR REGRESSION coefficient confidence interval])
1820
1821 AT_DATA([conf.sps], [dnl
1822 set format = F22.3.
1823
1824 data list notable list /math female socst read science *
1825 begin data.
1826     41.00       .00     57.00     57.00     47.00
1827     53.00      1.00     61.00     68.00     63.00
1828     54.00       .00     31.00     44.00     58.00
1829     47.00       .00     56.00     63.00     53.00
1830     57.00       .00     61.00     47.00     53.00
1831     51.00       .00     61.00     44.00     63.00
1832     42.00       .00     61.00     50.00     53.00
1833     45.00       .00     36.00     34.00     39.00
1834     54.00       .00     51.00     63.00     58.00
1835     52.00       .00     51.00     57.00     50.00
1836     51.00       .00     61.00     60.00     53.00
1837     51.00       .00     61.00     57.00     63.00
1838     71.00       .00     71.00     73.00     61.00
1839     57.00       .00     46.00     54.00     55.00
1840     50.00       .00     56.00     45.00     31.00
1841     43.00       .00     56.00     42.00     50.00
1842     51.00       .00     56.00     47.00     50.00
1843     60.00       .00     56.00     57.00     58.00
1844     62.00       .00     61.00     68.00     55.00
1845     57.00       .00     46.00     55.00     53.00
1846     35.00       .00     41.00     63.00     66.00
1847     75.00       .00     66.00     63.00     72.00
1848     45.00       .00     56.00     50.00     55.00
1849     57.00       .00     61.00     60.00     61.00
1850     45.00       .00     46.00     37.00     39.00
1851     46.00       .00     31.00     34.00     39.00
1852     66.00       .00     66.00     65.00     61.00
1853     57.00       .00     46.00     47.00     58.00
1854     49.00       .00     46.00     44.00     39.00
1855     49.00       .00     41.00     52.00     55.00
1856     57.00       .00     51.00     42.00     47.00
1857     64.00       .00     61.00     76.00     64.00
1858     63.00       .00     71.00     65.00     66.00
1859     57.00       .00     31.00     42.00     72.00
1860     50.00       .00     61.00     52.00     61.00
1861     58.00       .00     66.00     60.00     61.00
1862     75.00       .00     66.00     68.00     66.00
1863     68.00       .00     66.00     65.00     66.00
1864     44.00       .00     36.00     47.00     36.00
1865     40.00       .00     51.00     39.00     39.00
1866     41.00       .00     51.00     47.00     42.00
1867     62.00       .00     51.00     55.00     58.00
1868     57.00       .00     51.00     52.00     55.00
1869     43.00       .00     41.00     42.00     50.00
1870     48.00       .00     66.00     65.00     63.00
1871     63.00       .00     46.00     55.00     69.00
1872     39.00       .00     47.00     50.00     49.00
1873     70.00       .00     51.00     65.00     63.00
1874     63.00       .00     46.00     47.00     53.00
1875     59.00       .00     51.00     57.00     47.00
1876     61.00       .00     56.00     53.00     57.00
1877     38.00       .00     41.00     39.00     47.00
1878     61.00       .00     46.00     44.00     50.00
1879     49.00       .00     71.00     63.00     55.00
1880     73.00       .00     66.00     73.00     69.00
1881     44.00       .00     42.00     39.00     26.00
1882     42.00       .00     32.00     37.00     33.00
1883     39.00       .00     46.00     42.00     56.00
1884     55.00       .00     41.00     63.00     58.00
1885     52.00       .00     51.00     48.00     44.00
1886     45.00       .00     61.00     50.00     58.00
1887     61.00       .00     66.00     47.00     69.00
1888     39.00       .00     46.00     44.00     34.00
1889     41.00       .00     36.00     34.00     36.00
1890     50.00       .00     61.00     50.00     36.00
1891     40.00       .00     26.00     44.00     50.00
1892     60.00       .00     66.00     60.00     55.00
1893     47.00       .00     26.00     47.00     42.00
1894     59.00       .00     44.00     63.00     65.00
1895     49.00       .00     36.00     50.00     44.00
1896     46.00       .00     51.00     44.00     39.00
1897     58.00       .00     61.00     60.00     58.00
1898     71.00       .00     66.00     73.00     63.00
1899     58.00       .00     66.00     68.00     74.00
1900     46.00       .00     51.00     55.00     58.00
1901     43.00       .00     31.00     47.00     45.00
1902     54.00       .00     61.00     55.00     49.00
1903     56.00       .00     66.00     68.00     63.00
1904     46.00       .00     46.00     31.00     39.00
1905     54.00       .00     56.00     47.00     42.00
1906     57.00       .00     56.00     63.00     55.00
1907     54.00       .00     36.00     36.00     61.00
1908     71.00       .00     56.00     68.00     66.00
1909     48.00       .00     56.00     63.00     63.00
1910     40.00       .00     41.00     55.00     44.00
1911     64.00       .00     66.00     55.00     63.00
1912     51.00       .00     56.00     52.00     53.00
1913     39.00       .00     56.00     34.00     42.00
1914     40.00       .00     31.00     50.00     34.00
1915     61.00       .00     56.00     55.00     61.00
1916     66.00       .00     46.00     52.00     47.00
1917     49.00       .00     46.00     63.00     66.00
1918     65.00      1.00     61.00     68.00     69.00
1919     52.00      1.00     48.00     39.00     44.00
1920     46.00      1.00     51.00     44.00     47.00
1921     61.00      1.00     51.00     50.00     63.00
1922     72.00      1.00     56.00     71.00     66.00
1923     71.00      1.00     71.00     63.00     69.00
1924     40.00      1.00     41.00     34.00     39.00
1925     69.00      1.00     61.00     63.00     61.00
1926     64.00      1.00     66.00     68.00     69.00
1927     56.00      1.00     61.00     47.00     66.00
1928     49.00      1.00     41.00     47.00     33.00
1929     54.00      1.00     51.00     63.00     50.00
1930     53.00      1.00     51.00     52.00     61.00
1931     66.00      1.00     56.00     55.00     42.00
1932     67.00      1.00     56.00     60.00     50.00
1933     40.00      1.00     33.00     35.00     51.00
1934     46.00      1.00     56.00     47.00     50.00
1935     69.00      1.00     71.00     71.00     58.00
1936     40.00      1.00     56.00     57.00     61.00
1937     41.00      1.00     51.00     44.00     39.00
1938     57.00      1.00     66.00     65.00     46.00
1939     58.00      1.00     56.00     68.00     59.00
1940     57.00      1.00     66.00     73.00     55.00
1941     37.00      1.00     41.00     36.00     42.00
1942     55.00      1.00     46.00     43.00     55.00
1943     62.00      1.00     66.00     73.00     58.00
1944     64.00      1.00     56.00     52.00     58.00
1945     40.00      1.00     51.00     41.00     39.00
1946     50.00      1.00     51.00     60.00     50.00
1947     46.00      1.00     56.00     50.00     50.00
1948     53.00      1.00     56.00     50.00     39.00
1949     52.00      1.00     46.00     47.00     48.00
1950     45.00      1.00     46.00     47.00     34.00
1951     56.00      1.00     61.00     55.00     58.00
1952     45.00      1.00     56.00     50.00     44.00
1953     54.00      1.00     41.00     39.00     50.00
1954     56.00      1.00     46.00     50.00     47.00
1955     41.00      1.00     26.00     34.00     29.00
1956     54.00      1.00     56.00     57.00     50.00
1957     72.00      1.00     56.00     57.00     54.00
1958     56.00      1.00     51.00     68.00     50.00
1959     47.00      1.00     46.00     42.00     47.00
1960     49.00      1.00     66.00     61.00     44.00
1961     60.00      1.00     66.00     76.00     67.00
1962     54.00      1.00     46.00     47.00     58.00
1963     55.00      1.00     56.00     46.00     44.00
1964     33.00      1.00     41.00     39.00     42.00
1965     49.00      1.00     61.00     52.00     44.00
1966     43.00      1.00     51.00     28.00     44.00
1967     50.00      1.00     52.00     42.00     50.00
1968     52.00      1.00     51.00     47.00     39.00
1969     48.00      1.00     41.00     47.00     44.00
1970     58.00      1.00     66.00     52.00     53.00
1971     43.00      1.00     61.00     47.00     48.00
1972     41.00      1.00     31.00     50.00     55.00
1973     43.00      1.00     51.00     44.00     44.00
1974     46.00      1.00     41.00     47.00     40.00
1975     44.00      1.00     41.00     45.00     34.00
1976     43.00      1.00     46.00     47.00     42.00
1977     61.00      1.00     56.00     65.00     58.00
1978     40.00      1.00     51.00     43.00     50.00
1979     49.00      1.00     61.00     47.00     53.00
1980     56.00      1.00     66.00     57.00     58.00
1981     61.00      1.00     71.00     68.00     55.00
1982     50.00      1.00     61.00     52.00     54.00
1983     51.00      1.00     61.00     42.00     47.00
1984     42.00      1.00     41.00     42.00     42.00
1985     67.00      1.00     66.00     66.00     61.00
1986     53.00      1.00     61.00     47.00     53.00
1987     50.00      1.00     58.00     57.00     51.00
1988     51.00      1.00     31.00     47.00     63.00
1989     72.00      1.00     61.00     57.00     61.00
1990     48.00      1.00     61.00     52.00     55.00
1991     40.00      1.00     31.00     44.00     40.00
1992     53.00      1.00     61.00     50.00     61.00
1993     39.00      1.00     36.00     39.00     47.00
1994     63.00      1.00     41.00     57.00     55.00
1995     51.00      1.00     37.00     57.00     53.00
1996     45.00      1.00     43.00     42.00     50.00
1997     39.00      1.00     61.00     47.00     47.00
1998     42.00      1.00     39.00     42.00     31.00
1999     62.00      1.00     51.00     60.00     61.00
2000     44.00      1.00     51.00     44.00     35.00
2001     65.00      1.00     66.00     63.00     54.00
2002     63.00      1.00     71.00     65.00     55.00
2003     54.00      1.00     41.00     39.00     53.00
2004     45.00      1.00     36.00     50.00     58.00
2005     60.00      1.00     51.00     52.00     56.00
2006     49.00      1.00     51.00     60.00     50.00
2007     48.00      1.00     51.00     44.00     39.00
2008     57.00      1.00     61.00     52.00     63.00
2009     55.00      1.00     61.00     55.00     50.00
2010     66.00      1.00     56.00     50.00     66.00
2011     64.00      1.00     71.00     65.00     58.00
2012     55.00      1.00     51.00     52.00     53.00
2013     42.00      1.00     36.00     47.00     42.00
2014     56.00      1.00     61.00     63.00     55.00
2015     53.00      1.00     66.00     50.00     53.00
2016     41.00      1.00     41.00     42.00     42.00
2017     42.00      1.00     41.00     36.00     50.00
2018     53.00      1.00     56.00     50.00     55.00
2019     42.00      1.00     51.00     41.00     34.00
2020     60.00      1.00     56.00     47.00     50.00
2021     52.00      1.00     56.00     55.00     42.00
2022     38.00      1.00     46.00     42.00     36.00
2023     57.00      1.00     52.00     57.00     55.00
2024     58.00      1.00     61.00     55.00     58.00
2025     65.00      1.00     61.00     63.00     53.00
2026 end data.
2027
2028 regression
2029  /variables = math female socst read
2030  /statistics = coeff r anova ci (95)
2031  /dependent = science
2032  /method = enter
2033 ])
2034
2035 AT_CHECK([pspp -O format=csv conf.sps], [0], [dnl
2036 Table: Model Summary (science)
2037 ,R,R Square,Adjusted R Square,Std. Error of the Estimate
2038 ,.699,.489,.479,7.148
2039
2040 Table: ANOVA (science)
2041 ,,Sum of Squares,df,Mean Square,F,Sig.
2042 ,Regression,9543.721,4,2385.930,46.695,.000
2043 ,Residual,9963.779,195,51.096,,
2044 ,Total,19507.500,199,,,
2045
2046 Table: Coefficients (science)
2047 ,,Unstandardized Coefficients,,Standardized Coefficients,,,95% Confidence Interval for B,
2048 ,,B,Std. Error,Beta,t,Sig.,Lower Bound,Upper Bound
2049 ,(Constant),12.325,3.194,.000,3.859,.000,6.027,18.624
2050 ,math,.389,.074,.368,5.252,.000,.243,.535
2051 ,female,-2.010,1.023,-.101,-1.965,.051,-4.027,.007
2052 ,socst,.050,.062,.054,.801,.424,-.073,.173
2053 ,read,.335,.073,.347,4.607,.000,.192,.479
2054 ])
2055
2056
2057 AT_CLEANUP