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