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

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