COMMENT --------------------------------------- Betahat - Thurber---------------------------------------
INPATH C:\Projects\Code\Beta\Benchmarks\Nonlinear
OPEN Thurber_NIST.PRN
USE 1 37
NLS Y   (b1+b2*x+b3*x**2+b4*x**3)/(1+b5*x+b6*x**2+b7*x**3)  / INPUTB B1=1000 B2=1000 B3=400 B4=40   B5=0.7 B6=0.3 B7=0.03  
NLS Y    (b1+b2*x+b3*x**2+b4*x**3)/(1+b5*x+b6*x**2+b7*x**3)  / INPUTB B1=1300 B2=1500 B3=500 B4=75   B5=1 B6=0.4 B7=0.05 


COMMENT --------------------------------------- NIST - Thurber ---------------------------------------
'Source is http://www.itl.nist.gov/div898/strd/nls/data/thurber.shtml
COMMENT           Starting Values                  Certified Values
COMMENT 
COMMENT         Start 1     Start 2           Parameter     Standard Deviation
COMMENT   b1 =   1000        1300          1.2881396800E+03  4.6647963344E+00
COMMENT   b2 =   1000        1500          1.4910792535E+03  3.9571156086E+01
COMMENT   b3 =    400         500          5.8323836877E+02  2.8698696102E+01
COMMENT   b4 =     40          75          7.5416644291E+01  5.5675370270E+00
COMMENT   b5 =      0.7         1          9.6629502864E-01  3.1333340687E-02
COMMENT   b6 =      0.3         0.4        3.9797285797E-01  1.4984928198E-02
COMMENT   b7 =      0.03        0.05       4.9727297349E-02  6.5842344623E-03
COMMENT 
COMMENT Residual Sum of Squares:                    5.6427082397E+03
COMMENT Residual Standard Deviation:                1.3714600784E+01
COMMENT Degrees of Freedom:                                30
COMMENT Number of Observations:                            37