Polyfit

Greg Willden gregwillden at gmail.com
Thu Oct 12 16:15:34 CDT 2006


On 10/12/06, Charles R Harris <charlesr.harris at gmail.com> wrote:
>
> I'm guessing that the rcond number in the lstsq version (default 1e-10) is
> the difference. Generally the lstsq version should work better than the MPL
> version because at*a is not as well conditioned and vandermonde matrices are
> notoriously ill conditioned anyway for higher degree fits. It would help if
> you attached the data files in ascii form unless they happen to contain
> thousands of data items. Just the x will suffice, zip it up if you have to.
>


Here are the files.

Since the two algorithms behave differently and each has it place then can
both be included in numpy?
i.e. numpy.polyfit(x,y,N, mode='alg1')
numpy.polyfit(x,y,N, mode='alg2')

replacing alg1 and alg2 with meaningful names.

Thanks
Greg
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://projects.scipy.org/pipermail/numpy-discussion/attachments/20061012/743d056a/attachment.html 
-------------- next part --------------
2
4
0
16
1
3
14
17
12
13
8
7
5
11
23
19
24
9
18
15
44
22
20
25
27
10
41
45
26
40
21
50
29
38
28
37
53
48
57
78
54
61
39
64
65
58
43
35
55
47
60
66
46
68
52
51
42
75
93
56
87
74
59
71
84
49
72
81
73
63
94
67
36
80
30
70
106
69
62
77
76
90
98
103
6
172
79
83
82
31
85
97
105
91
121
89
95
101
116
135
92
102
96
88
120
137
136
34
389
499
489
565
479
555
412
390
443
442
456
365
404
593
457
364
375
441
367
360
453
376
482
576
550
445
409
658
363
435
381
477
549
396
460
544
601
466
384
557
391
461
511
642
547
558
542
610
608
639
438
591
551
378
599
552
397
569
578
437
643
361
385
399
538
600
395
561
602
649
454
464
447
386
535
554
577
596
686
635
459
514
628
448
583
407
408
603
379
626
410
500
570
449
582
434
401
417
501
462
566
366
394
631
485
382
563
369
784
781
824
838
783
826
825
854
867
763
801
799
879
795
759
852
872
831
802
769
873
773
762
779
772
789
836
881
839
827
796
882
851
786
797
755
726
841
896
865
890
837
857
775
853
767
785
756
750
754
887
861
894
770
749
771
840
844
782
780
877
859
774
832
871
730
722
901
794
884
846
760
848
729
876
724
720
757
830
798
850
818
764
875
787
829
899
885
874
960
815
864
777
723
820
891
747
751
725
821
776
823
866
800
833
788
758
878
1411
1350
1358
1353
1362
1378
1326
1348
1321
1355
1357
1349
1352
1346
1373
1430
1408
1351
1310
1418
1281
1336
1416
1428
1423
1436
1364
1432
1419
1377
1421
1359
1279
1401
1424
1323
1340
1322
1339
1407
1438
1427
1247
1435
1290
1311
1262
1343
1307
1437
1369
1315
1288
1356
1412
1309
1327
1276
1367
1325
1425
1420
1344
1361
1337
1368
1314
1415
1335
1278
1320
1280
1382
1282
1434
1328
1273
1332
1410
1347
1284
1338
1417
1354
1360
1329
1270
1305
1375
1414
1287
1291
1429
1379
1256
1372
1370
1289
1306
1275
1342
1274
1286
1313
1271
1248
1319
1268
1479
1440
1452
1449
1441
1542
1444
1447
1645
1569
1485
1537
1527
1522
1462
1461
1554
1651
1448
1652
1538
1529
1510
1451
1639
1570
1486
1737
1481
1623
1725
1453
1458
1459
1680
1484
1445
1671
1460
1524
1673
1666
1627
1540
1564
1626
1672
1499
1674
1561
1555
1519
1472
1632
1534
1548
1650
1470
1656
1637
1648
1717
1661
1521
1724
1755
1631
1735
1517
1667
1551
1728
1566
1575
1515
1539
1531
1541
1563
1558
1643
1454
1549
1758
1483
1761
1714
1638
1505
1654
1664
1475
1556
1473
1528
1516
1552
1443
1709
1553
1794
1571
1644
1520
1718
1742
1498
1736
2149
2052
2128
2141
2019
2108
2125
1998
2055
2087
2106
2123
2029
2036
1960
2121
2049
1945
2132
2058
2035
2144
2053
2109
1963
1949
1966
2137
2060
2153
2147
1956
2126
2115
2130
2129
1808
2021
2038
1804
2064
2104
1952
2107
1911
1958
2116
1810
2031
1940
2148
2117
2150
1814
2059
1916
2010
2026
1809
1964
2103
2032
2039
2118
2151
1973
2152
1951
2084
2011
2065
2105
2028
2112
2044
2040
2122
2033
1812
1941
1946
1987
1996
2127
2113
2048
1995
1950
2155
1959
2138
1803
1970
1961
2020
2008
1811
1968
1965
1947
1819
2037
2063
1817
2133
1967
2140
2047
2185
2161
2217
2241
2163
2164
2212
2220
2247
2162
2208
2248
2225
2255
2206
2191
2218
2165
2238
2237
2166
2231
2211
2192
2258
2233
2201
2221
2239
2251
2249
2219
2186
2197
2246
2210
2199
2184
2226
2243
2188
2240
2250
2204
2205
2193
2260
2160
2236
2187
2228
2182
2223
2215
2224
2235
2213
2198
2168
2202
2195
2183
2253
2209
2234
2288
2261
2256
2254
2200
2257
2242
2167
2280
2252
2245
2207
2263
2216
2277
2284
2244
2190
2259
2266
2267
2180
2227
2196
2229
2297
2189
2169
2296
2181
2279
2281
2179
2178
2194
2287
2222
2285
2232
2301
2264
2214
2289
2787
2714
2823
2591
2640
2717
2650
2612
2838
2630
2827
2588
2685
2684
2677
2720
2613
2872
2729
2804
2622
2841
2642
2694
2867
2693
2865
2601
2808
2870
2691
2607
2877
2620
2836
2835
2576
2596
2636
2805
2802
2603
2579
2797
2610
2577
2716
2832
2626
2676
2678
2873
2785
2632
2608
2692
2643
2726
2813
2614
2675
2791
2803
2644
2701
2801
2594
2737
2625
2826
2839
2649
2700
2799
2634
2695
2730
2592
2874
2667
2647
2683
2586
2722
2679
2669
2629
2871
2549
2733
2582
2631
2807
2606
2641
2727
2724
2652
2623
2783
2646
2821
2602
2702
2708
2589
2815
2842
3222
3230
3229
2898
3228
3214
3232
3208
3213
2931
3202
2886
2922
3205
3221
3235
3212
3204
2887
2893
3209
2907
2921
3231
2891
2914
3223
3224
2890
2911
3189
2896
3233
2909
2917
3220
3198
3207
3226
3227
2900
2889
3225
2923
2884
3203
3192
2910
2926
3210
2885
2905
2913
2895
2897
2918
2924
3200
2927
2934
3234
2912
3206
3187
3236
3196
2881
3238
3190
3219
2899
2936
3237
3199
2919
2915
2882
3215
3188
2902
2937
2906
2930
2901
2883
3197
2908
2904
3193
3201
3218
2925
2892
2916
3194
2939
3195
2888
3191
2880
2894
2935
2932
2944
2920
2928
2929
2938
3581
3576
3572
3571
3568
3575
3569
3566
3562
3574
3577
3565
3570
3559
3550
3582
3578
3585
3583
3560
3573
3563
3580
3588
3586
3584
3579
3564
3567
3561
3549
3587
3551
3590
3591
3589
3558
3515
3523
3511
3526
3552
3512
3547
3510
3508
3520
3503
3524
3509
3518
3557
3545
3548
3514
3592
3517
3519
3528
3513
3522
3488
3521
3487
3516
3544
3497
3553
3507
3525
3504
3546
3499
3490
3527
3492
3494
3501
3489
3500
3593
3491
3496
3311
3458
3306
3594
3495
3394
3543
3498
3277
3486
3502
3370
3505
3308
3460
3469
3459
3466
3483
3366
3482
3269
3323
3264
3485
-------------- next part --------------
3956.970544
4114.773479
4313.241696
4384.792216
4391.872310
4411.476715
4482.940096
4594.108179
4645.381743
4730.829231
4763.895598
4781.032135
4880.854414
5008.228140
5013.866252
5026.548238
5040.509458
5061.949397
5098.508781
5118.791696
5162.805554
5186.805236
5193.941710
5256.845108
5286.838670
5364.175004
5413.297293
5418.790951
5446.085519
5475.639165
5508.616723
5604.737760
5614.409364
5636.101761
5638.542797
5649.773243
5660.913105
5719.715027
5815.218454
5832.621843
5837.578731
5862.529421
5863.607227
5877.191757
5952.905334
5984.647641
6033.194005
6048.375830
6058.973718
6073.108082
6100.784624
6112.341023
6138.705305
6142.491247
6170.101238
6190.399915
6224.593338
6239.605558
6248.244072
6254.401594
6280.211422
6317.197093
6317.614065
6325.296754
6339.852487
6381.476997
6383.259524
6391.571145
6397.082360
6418.080290
6438.735101
6440.724082
6454.977105
6477.089430
6560.501879
6637.842362
6710.064249
6724.529239
6728.116095
6728.478821
6774.863744
6874.242169
6902.087390
6919.527972
6979.799498
6980.973237
6989.850615
7015.371907
7019.392350
7050.062660
7050.597597
7061.397784
7105.324281
7107.397379
7113.332893
7123.788680
7132.977337
7138.892938
7172.808848
7176.242427
7186.719367
7191.745735
7239.760194
7243.749478
7255.868073
7257.267870
7268.538062
7272.629743
7679.940589
7683.236804
7780.699438
7873.928129
7899.908299
7904.060156
7925.056466
7986.214500
7992.112331
7999.563794
8001.244799
8018.175375
8030.246998
8042.876999
8066.748960
8068.723221
8070.288644
8071.628413
8130.859172
8156.118329
8175.197378
8206.961155
8218.103232
8222.490958
8239.676758
8241.640380
8249.615006
8267.669676
8268.053374
8283.319127
8288.302563
8293.092398
8316.055561
8327.772501
8334.822552
8338.317644
8339.073386
8342.683772
8350.156751
8357.692104
8400.200111
8408.910930
8417.560026
8424.256975
8424.301520
8435.165529
8442.029865
8463.791695
8467.858939
8468.002339
8473.892687
8475.457054
8481.274131
8484.884301
8503.376563
8506.217207
8513.976462
8516.443233
8518.839920
8521.927566
8522.369277
8524.426356
8529.354446
8530.550718
8537.266324
8543.730529
8554.276036
8555.003337
8573.504421
8575.663769
8578.398385
8586.508857
8597.926278
8601.801143
8607.638604
8613.355557
8618.223021
8619.832545
8625.856397
8628.547664
8635.510731
8635.616818
8638.093330
8638.725205
8639.290086
8640.963579
8641.074684
8644.931541
8650.771065
8651.748195
8652.153167
8654.191277
8662.647644
8663.908731
8672.156947
8673.406207
8673.653387
8676.156890
8684.297953
8688.001563
8691.611803
8691.714375
8692.654871
8705.209794
8714.218292
8715.503995
8720.084210
8720.465880
8950.316961
9014.234275
9016.582013
9068.120557
9079.086925
9084.084293
9092.854202
9152.516280
9206.442055
9211.797882
9259.935823
9274.088555
9319.200628
9418.736103
9426.658011
9437.435612
9448.483510
9448.919497
9451.123274
9475.891179
9488.585768
9492.490822
9513.276633
9516.127121
9521.205075
9530.898966
9564.375640
9588.755348
9596.967298
9615.132594
9642.107437
9645.510988
9646.175388
9666.070621
9677.726960
9695.380196
9701.239340
9709.602677
9711.670031
9739.504151
9761.667030
9761.780136
9769.510284
9799.923976
9807.513067
9809.803762
9810.694795
9823.853204
9838.334226
9848.403779
9876.275178
9887.129520
9887.560756
9887.888122
9895.889715
9899.907916
9909.691312
9911.111820
9949.714270
9952.911904
9960.038486
9962.479402
9971.996514
9980.326000
9984.951467
9994.006658
10001.073320
10003.476996
10014.757757
10022.105369
10025.979253
10049.807479
10054.324996
10058.877706
10058.904792
10085.774604
10088.651241
10091.809589
10097.187658
10107.697162
10109.684042
10113.805433
10114.595009
10137.018824
10138.092392
10150.691072
10151.749667
10171.572498
10177.785950
10182.552424
10185.711350
10193.402181
10198.124474
10202.876980
10203.927967
10204.398934
10213.144657
10229.042155
10235.314449
10238.453722
10241.626409
10241.820711
10263.338102
10266.747743
10319.492303
10332.302989
10339.341659
10343.745803
8793.912077
9240.132446
9371.045485
9559.747528
9566.311321
9605.890226
9663.291742
9727.898364
9817.130726
9829.579358
9894.337029
9943.329879
9986.701787
9988.681733
10007.595675
10039.175362
10058.567031
10071.732150
10080.223263
10080.572538
10082.640448
10091.338544
10096.808553
10104.572156
10196.766192
10243.213998
10247.809490
10253.163862
10292.904202
10312.906320
10315.642277
10327.044105
10342.839692
10346.473991
10355.275639
10362.684038
10372.115530
10372.748272
10382.069791
10387.721477
10389.628136
10399.630692
10430.582283
10442.936940
10449.245767
10454.870137
10472.091864
10478.317998
10489.109327
10509.306150
10513.053502
10514.097318
10525.194284
10527.980647
10529.909410
10531.238327
10537.874168
10544.844175
10555.905185
10563.925159
10565.773122
10575.339600
10586.124716
10590.426276
10611.035077
10622.729994
10636.607153
10657.784468
10666.334794
10672.555469
10686.503366
10686.595292
10693.772247
10703.728069
10703.988774
10704.045868
10708.960756
10720.027629
10738.600374
10777.951535
10783.717838
10797.095152
10803.299357
10804.900797
10805.415021
10809.890716
10825.610894
10830.556433
10854.840855
10861.376313
10885.763384
10918.933837
10961.511982
10966.953981
10967.776642
10973.773434
10977.243562
10978.798310
10978.830169
10989.002409
11009.161106
11015.743459
11020.016124
11056.625898
11092.348893
11093.400822
11100.135715
11153.180487
9814.960894
9981.143054
10098.094627
10098.936901
10251.240937
10282.032995
10284.640897
10322.598305
10360.982393
10376.872740
10392.495710
10417.277021
10425.686178
10434.329513
10450.355048
10478.041847
10499.163586
10499.215862
10526.877232
10532.554246
10536.409254
10539.253460
10608.396814
10645.575298
10659.240945
10670.922154
10689.184454
10715.836202
10722.298297
10725.640874
10748.404608
10753.162544
10756.517743
10757.745888
10772.695671
10782.428719
10790.583373
10791.898039
10801.684370
10811.977492
10816.266996
10818.283291
10820.330119
10834.649198
10857.733858
10874.140510
10876.059901
10877.637723
10882.405065
10894.207293
10900.073600
10912.167792
10924.820170
10926.306574
10929.257567
10959.628651
10960.178273
10975.859148
10976.998078
10978.812143
10991.596620
10999.187685
11000.809537
11010.126635
11024.307375
11071.514474
11085.007968
11094.736694
11109.770745
11113.527504
11116.002506
11122.254903
11133.988453
11136.612348
11156.603051
11165.834508
11166.019724
11167.998566
11175.515681
11193.546241
11199.134090
11238.198692
11244.422619
11253.787513
11264.283982
11271.505442
11274.799495
11299.434399
11316.162315
11317.946622
11318.921151
11319.649335
11346.282844
11358.697447
11362.651974
11372.919362
11374.332828
11382.460151
11394.420472
11396.858534
11397.989178
11401.709776
11418.079306
11432.017976
11438.011867
11459.901142
11468.659319
11472.542199
10681.749907
10956.703906
11125.372074
11136.314500
11172.483691
11178.012600
11269.727750
11300.670178
11340.773673
11383.722922
11457.659862
11462.856526
11489.727744
11551.874237
11579.544500
11637.402145
11686.213732
11690.080383
11741.072475
11747.628633
11787.678534
11801.682051
11814.371802
11829.473150
11838.299686
11858.353106
11870.960722
11871.392282
11871.694919
11876.150858
11883.402109
11895.260550
11913.376985
11913.856733
11942.894620
11951.949830
11969.258031
11983.875254
11993.436462
11995.024242
12021.141064
12023.333762
12048.485972
12050.270337
12055.970075
12066.039635
12073.944508
12093.030117
12098.907173
12116.424190
12116.649970
12120.302248
12124.960350
12136.994784
12137.012910
12150.298122
12150.452929
12166.326678
12170.039350
12171.669367
12173.044335
12174.841150
12179.606648
12207.744517
12209.538641
12225.212482
12229.865951
12232.316882
12246.263509
12256.108638
12276.069698
12286.868822
12297.149059
12306.631699
12312.322071
12320.164640
12324.374077
12336.310756
12339.376427
12368.598153
12380.799010
12382.583485
12384.933702
12389.988805
12391.470566
12395.745999
12401.831477
12402.594363
12413.304896
12413.370202
12444.759649
12448.495831
12459.577196
12470.502275
12470.742632
12472.315225
12474.617021
12478.433235
12478.483165
12482.236078
12498.028231
12513.399284
12525.231137
12532.312092
12534.743170
12537.889753
12542.266646
12551.826599
11193.319247
11514.884367
11602.445732
11622.344038
11622.438404
11630.019170
11641.485017
11667.278665
11708.588012
11789.904842
11836.441818
11850.623879
11881.255290
11925.047593
11951.347144
11964.909208
11994.423115
12038.711459
12090.526264
12106.978903
12160.550605
12177.736209
12215.007376
12230.055459
12250.739418
12257.282265
12259.861311
12260.730253
12262.735060
12291.224191
12316.483640
12320.419844
12358.768148
12371.123837
12374.651193
12388.097106
12445.256592
12448.681215
12473.311313
12504.297769
12505.134285
12527.592355
12575.220793
12625.944644
12651.232120
12654.872920
12667.986158
12673.676463
12698.035009
12713.500228
12766.790350
12795.378324
12811.393790
12816.705173
12821.887367
12836.305776
12859.372179
12868.193924
12869.219546
12883.997847
12913.567482
12928.892730
12951.562632
12957.166348
13033.097686
13039.613435
13052.920297
13103.059009
13121.399000
13137.420222
13143.819581
13189.273853
13217.056695
13227.901100
13251.107604
13252.729284
13256.649463
13269.470154
13291.730794
13301.492326
13306.257808
13307.889142
13360.021884
13423.017549
13429.234543
13437.179872
13484.090893
13496.163341
13595.522328
13606.157527
13609.179758
13617.822523
13619.608609
13690.242120
13717.833712
13735.254554
13755.434901
13807.103472
13824.931396
13827.734231
13834.988156
13914.748586
13978.380023
14060.743382
14083.144765
14103.587172
14155.685069
14169.683047
12042.398439
12228.734395
12520.514566
12535.035045
12580.133137
12590.412616
12629.404482
12654.158926
12675.304854
12733.821411
12802.369761
12809.260088
12868.977918
12891.461146
12902.360057
12944.236947
12970.555105
12984.851197
12993.054529
13024.315480
13028.277354
13036.186899
13065.224983
13074.307427
13077.734413
13120.257899
13138.725590
13139.563867
13144.916995
13147.708541
13150.495330
13150.670526
13153.697961
13158.210056
13170.801478
13173.506940
13174.016751
13175.306920
13181.264005
13182.217312
13189.197422
13191.435130
13212.256657
13212.308700
13216.024874
13216.050142
13217.485968
13262.914785
13271.722595
13282.591589
13291.680677
13313.503687
13317.334527
13331.044319
13332.073092
13332.970644
13351.033630
13367.919842
13372.757479
13383.369447
13400.388651
13403.526696
13406.483139
13410.260061
13413.739986
13423.056454
13427.512376
13428.578560
13431.234935
13438.083882
13440.731763
13467.215303
13471.345252
13475.473964
13476.412131
13481.174967
13491.373425
13496.015255
13499.707775
13500.224718
13512.408367
13515.414742
13519.258531
13519.955607
13529.853371
13555.701103
13567.556673
13581.752972
13592.412676
13612.170183
13618.520577
13621.250454
13627.232255
13638.740682
13644.590236
13645.126659
13649.339579
13652.869745
13660.267656
13663.052077
13666.620509
13667.808315
13667.827996
13678.639363
13685.255280
13689.203630
13691.048025
13713.641976
11598.693262
11704.336246
11919.040078
11981.528008
12057.828308
12094.424616
12149.980041
12179.583269
12183.210365
12251.745007
12262.323134
12331.253864
12412.695494
12418.496422
12425.118997
12528.458621
12550.929713
12566.605855
12624.572874
12625.259561
12647.308193
12679.031675
12683.923804
12704.924039
12745.050992
12770.103628
12787.591668
12809.165792
12819.808330
12822.677224
12835.990692
12854.239067
12913.753810
12933.290259
12939.277223
12959.460061
12965.124068
12971.594291
12994.159714
13007.621915
13010.606602
13018.524252
13043.586740
13060.549764
13069.738691
13081.265005
13083.369723
13111.744055
13132.289855
13231.443965
13243.583985
13306.335621
13307.422095
13315.641892
13373.612750
13379.343564
13437.392770
13486.774587
13498.942593
13516.270938
13517.006573
13617.129826
13651.302768
13655.976384
13659.231345
13669.637967
13692.471342
13702.477965
13709.766221
13724.485954
13748.119265
13764.708728
13780.746631
13821.058818
13822.464110
13864.480461
13867.177653
13868.782455
13869.689934
13920.956099
13928.184482
13975.859048
13985.023508
13988.347825
13992.894227
14009.444963
14037.708338
14058.926075
14150.438533
14172.273453
14186.946281
14194.290515
14246.683826
14305.245169
14308.202607
14310.050716
14351.339722
14376.973672
14392.142221
14397.889970
14434.301413
14436.668387
14454.755269
14546.002108
14564.107635
14609.352294
14610.444061
14658.388376
6191.820693
6417.995072
6450.171771
6639.135507
6702.781151
6808.454342
6858.701101
6875.949295
6896.356613
6899.675204
6980.449533
7001.842011
7076.463311
7081.714631
7082.396339
7176.817676
7216.309297
7266.648313
7280.514286
7297.524297
7355.848621
7373.904267
7378.860143
7408.329335
7458.162090
7556.879500
7567.809539
7642.724706
7735.706825
8038.110570
8048.230497
8061.020685
8148.930937
8210.087296
8262.783633
8377.205167
8474.150854
8624.533634
8822.084093
8902.373733
9086.760508
9093.644637
9159.773432
9176.508060
9188.687554
9189.349289
9191.366199
9192.863741
9292.189993
9293.081584
9297.837340
9391.546496
9403.902790
9423.233460
9528.854294
9531.233339
9555.442442
9596.527664
9621.645865
9664.165188
9748.118478
9801.103502
9845.340170
9883.906937
9910.602414
9914.261309
9914.838154
9944.961769
9947.652699
10030.420720
10032.086427
10085.772122
10106.567350
10131.776082
10171.059604
10171.271504
10250.514131
10294.686428
10310.671898
10354.773039
10359.713965
10417.971882
10419.604148
10518.670282
10543.442886
10604.550121
10604.651083
10638.592841
10715.252457
10750.162225
10763.995006
10825.711128
10857.760357
10872.843194
10894.852725
10910.851029
10927.298566
10932.509272
10938.224926
11014.452664
11019.602796
11031.241511
11052.111308
11055.718676
11153.896765
11160.080716
11179.872948
11193.467570
-------------- next part --------------
-------------------------------------------------------------------------
Using Tomcat but need to do more? Need to support web services, security?
Get stuff done quickly with pre-integrated technology to make your job easier
Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo
http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
-------------- next part --------------
_______________________________________________
Numpy-discussion mailing list
Numpy-discussion at lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/numpy-discussion


More information about the Numpy-discussion mailing list