********************************************************* * ODRPACK95 version 1.00 of 12-27-2005 (REAL (KIND=R8)) * ********************************************************* *** Initial summary for fit by method of ODR *** --- Problem Size: N = 10 (number with nonzero weight = 10) NQ = 1 M = 1 NP = 2 (number unfixed = 2) --- Control Values: JOB = 00000 = ABCDE, where A=0 ==> fit is not a restart. B=0 ==> deltas are initialized to zero. C=0 ==> covariance matrix will be computed using derivatives re-evaluated at the solution. D=0 ==> derivatives are estimated by forward differences. E=0 ==> method is explicit ODR. NDIGIT = 16 (estimated by ODRPACK95) TAUFAC = 1.00E+00 --- Stopping Criteria: SSTOL = 1.49E-08 (sum of squares stopping tolerance) PARTOL = 3.67E-11 (parameter stopping tolerance) MAXIT = 1000 (maximum number of iterations) --- Initial Weighted Sum of Squares = NaN Sum of Squared Weighted Deltas = 0.00000000E+00 Sum of Squared Weighted Epsilons = NaN *** Final summary for fit by method of ODR *** --- Stopping Conditions: INFO = 60000 = ABCDE, where a nonzero value for a given digit indicates an abnormal stopping condition. A=6 ==> numerical instabilities have been detected, possibly indicating a discontinuity in the derivatives or a poor poor choice of problem scale or weights. NITER = 1 (number of iterations) NFEV = 105 (number of function evaluations) IRANK = 0 (rank deficiency) RCOND = NaN (inverse condition number) ISTOP = 0 (returned by user from subroutine FCN) --- Final Weighted Sums of Squares = NaN Sum of Squared Weighted Deltas = 0.00000000E+00 Sum of Squared Weighted Epsilons = NaN --- Estimated BETA(J), J = 1, ..., NP: N.B. the standard errors of the estimated betas were not computed. (see info above.) Index Value --------------> 1 to 2 NaN NaN N.B. no parameters were fixed by the user or dropped at the last iteration because they caused the model to be rank deficient. --- Estimated EPSILON(I) and DELTA(I,*), I = 1, ..., N: I EPSILON(I,1) DELTA(I,1) 1 NaN 0.00000000E+00 2 NaN 0.00000000E+00 3 NaN 0.00000000E+00 4 NaN 0.00000000E+00 5 NaN 0.00000000E+00 6 NaN 0.00000000E+00 7 NaN 0.00000000E+00 8 NaN 0.00000000E+00 9 NaN 0.00000000E+00 10 NaN 0.00000000E+00 ********************************************************* * ODRPACK95 version 1.00 of 12-27-2005 (REAL (KIND=R8)) * ********************************************************* *** Initial summary for fit by method of ODR *** --- Problem Size: N = 10 (number with nonzero weight = 10) NQ = 1 M = 1 NP = 2 (number unfixed = 2) --- Control Values: JOB = 00000 = ABCDE, where A=0 ==> fit is not a restart. B=0 ==> deltas are initialized to zero. C=0 ==> covariance matrix will be computed using derivatives re-evaluated at the solution. D=0 ==> derivatives are estimated by forward differences. E=0 ==> method is explicit ODR. NDIGIT = 16 (estimated by ODRPACK95) TAUFAC = 1.00E+00 --- Stopping Criteria: SSTOL = 1.49E-08 (sum of squares stopping tolerance) PARTOL = 3.67E-11 (parameter stopping tolerance) MAXIT = 1000 (maximum number of iterations) --- Initial Weighted Sum of Squares = NaN Sum of Squared Weighted Deltas = 0.00000000E+00 Sum of Squared Weighted Epsilons = NaN *** Final summary for fit by method of ODR *** --- Stopping Conditions: INFO = 60000 = ABCDE, where a nonzero value for a given digit indicates an abnormal stopping condition. A=6 ==> numerical instabilities have been detected, possibly indicating a discontinuity in the derivatives or a poor poor choice of problem scale or weights. NITER = 1 (number of iterations) NFEV = 110 (number of function evaluations) IRANK = 0 (rank deficiency) RCOND = NaN (inverse condition number) ISTOP = 0 (returned by user from subroutine FCN) --- Final Weighted Sums of Squares = NaN Sum of Squared Weighted Deltas = 0.00000000E+00 Sum of Squared Weighted Epsilons = NaN --- Estimated BETA(J), J = 1, ..., NP: N.B. the standard errors of the estimated betas were not computed. (see info above.) Index Value --------------> 1 to 2 NaN NaN N.B. no parameters were fixed by the user or dropped at the last iteration because they caused the model to be rank deficient. --- Estimated EPSILON(I) and DELTA(I,*), I = 1, ..., N: I EPSILON(I,1) DELTA(I,1) 1 NaN 0.00000000E+00 2 NaN 0.00000000E+00 3 NaN 0.00000000E+00 4 NaN 0.00000000E+00 5 NaN 0.00000000E+00 6 NaN 0.00000000E+00 7 NaN 0.00000000E+00 8 NaN 0.00000000E+00 9 NaN 0.00000000E+00 10 NaN 0.00000000E+00 ********************************************************* * ODRPACK95 version 1.00 of 12-27-2005 (REAL (KIND=R8)) * ********************************************************* *** Initial summary for fit by method of ODR *** --- Problem Size: N = 10 (number with nonzero weight = 10) NQ = 1 M = 1 NP = 1 (number unfixed = 1) --- Control Values: JOB = 00000 = ABCDE, where A=0 ==> fit is not a restart. B=0 ==> deltas are initialized to zero. C=0 ==> covariance matrix will be computed using derivatives re-evaluated at the solution. D=0 ==> derivatives are estimated by forward differences. E=0 ==> method is explicit ODR. NDIGIT = 16 (estimated by ODRPACK95) TAUFAC = 1.00E+00 --- Stopping Criteria: SSTOL = 1.49E-08 (sum of squares stopping tolerance) PARTOL = 3.67E-11 (parameter stopping tolerance) MAXIT = 1000 (maximum number of iterations) --- Initial Weighted Sum of Squares = 1.31464741E+02 Sum of Squared Weighted Deltas = 0.00000000E+00 Sum of Squared Weighted Epsilons = 1.31464741E+02 *** Iteration reports for fit by method of ODR *** Cum. Act. Rel. Pred. Rel. It. No. FN Weighted Sum-of-Sqs Sum-of-Sqs G-N Num. Evals Sum-of-Sqs Reduction Reduction TAU/PNORM Step ---- ------ ----------- ----------- ----------- --------- ---- 1 8 1.16106E+01 9.1168E-01 9.7286E-01 4.645E-01 YES 2 11 3.88003E+00 6.6582E-01 6.7897E-01 1.730E-01 YES 3 14 3.76019E+00 3.0888E-02 3.1072E-02 2.794E-02 YES 4 17 3.76016E+00 8.0881E-06 8.1371E-06 4.604E-04 YES 5 20 3.76016E+00 3.4237E-10 3.5320E-10 3.686E-06 YES *** Final summary for fit by method of ODR *** --- Stopping Conditions: INFO = 1 ==> sum of squares convergence. NITER = 5 (number of iterations) NFEV = 22 (number of function evaluations) IRANK = 0 (rank deficiency) RCOND = 1.00E+00 (inverse condition number) ISTOP = 0 (returned by user from subroutine FCN) --- Final Weighted Sums of Squares = 3.76015539E+00 Sum of Squared Weighted Deltas = 1.71977286E+00 Sum of Squared Weighted Epsilons = 2.04038253E+00 --- Residual Standard Deviation = 6.46370670E-01 Degrees of Freedom = 9 --- Estimated BETA(J), J = 1, ..., NP: BETA LOWER UPPER S.D. ___ 95% Confidence ___ BETA Interval 1 4.09967822E-01 -1.80+308 1.80+308 1.20E-02 3.83E-01 to 4.37E-01 --- Estimated EPSILON(I) and DELTA(I,*), I = 1, ..., N: I EPSILON(I,1) DELTA(I,1) 1 6.37547286E-03 -1.45698265E-02 2 -1.14912603E-01 2.79417569E-03 3 1.62319283E-01 1.66708947E-03 4 -1.25177679E-01 -3.50463316E-03 5 4.98123683E-02 5.67737287E-03 6 1.63669850E-02 3.68250284E-03 7 -2.56066391E-02 -1.07081722E-02 8 -2.53298509E-02 -3.06023453E-02 9 3.97263837E-03 1.03745280E-01 10 5.40429771E-05 7.12863325E-01 ********************************************************* * ODRPACK95 version 1.00 of 12-27-2005 (REAL (KIND=R8)) * ********************************************************* *** Initial summary for fit by method of ODR *** --- Problem Size: N = 10 (number with nonzero weight = 10) NQ = 1 M = 1 NP = 1 (number unfixed = 1) --- Control Values: JOB = 00000 = ABCDE, where A=0 ==> fit is not a restart. B=0 ==> deltas are initialized to zero. C=0 ==> covariance matrix will be computed using derivatives re-evaluated at the solution. D=0 ==> derivatives are estimated by forward differences. E=0 ==> method is explicit ODR. NDIGIT = 16 (estimated by ODRPACK95) TAUFAC = 1.00E+00 --- Stopping Criteria: SSTOL = 1.49E-08 (sum of squares stopping tolerance) PARTOL = 3.67E-11 (parameter stopping tolerance) MAXIT = 1000 (maximum number of iterations) --- Initial Weighted Sum of Squares = 1.45361119E+02 Sum of Squared Weighted Deltas = 0.00000000E+00 Sum of Squared Weighted Epsilons = 1.45361119E+02 *** Iteration reports for fit by method of ODR *** Cum. Act. Rel. Pred. Rel. It. No. FN Weighted Sum-of-Sqs Sum-of-Sqs G-N Num. Evals Sum-of-Sqs Reduction Reduction TAU/PNORM Step ---- ------ ----------- ----------- ----------- --------- ---- 1 8 1.24190E+01 9.1456E-01 9.7598E-01 4.601E-01 YES 2 11 3.90185E+00 6.8581E-01 7.0020E-01 1.838E-01 YES 3 14 3.75939E+00 3.6512E-02 3.6732E-02 3.116E-02 YES 4 17 3.75934E+00 1.2071E-05 1.2145E-05 5.843E-04 YES 5 20 3.75934E+00 5.4746E-10 5.7337E-10 4.481E-06 YES *** Final summary for fit by method of ODR *** --- Stopping Conditions: INFO = 1 ==> sum of squares convergence. NITER = 5 (number of iterations) NFEV = 22 (number of function evaluations) IRANK = 0 (rank deficiency) RCOND = 1.00E+00 (inverse condition number) ISTOP = 0 (returned by user from subroutine FCN) --- Final Weighted Sums of Squares = 3.75934199E+00 Sum of Squared Weighted Deltas = 1.71789508E+00 Sum of Squared Weighted Epsilons = 2.04144691E+00 --- Residual Standard Deviation = 6.46300755E-01 Degrees of Freedom = 9 --- Estimated BETA(J), J = 1, ..., NP: BETA LOWER UPPER S.D. ___ 95% Confidence ___ BETA Interval 1 1.00038730E+00 -1.80+308 1.80+308 8.73E-03 9.81E-01 to 1.02E+00 --- Estimated EPSILON(I) and DELTA(I,*), I = 1, ..., N: I EPSILON(I,1) DELTA(I,1) 1 6.33525046E-03 -1.44492478E-02 2 -1.17910238E-01 2.86328356E-03 3 1.59524739E-01 1.62089019E-03 4 -1.27260221E-01 -3.55202417E-03 5 4.88188626E-02 5.55537178E-03 6 1.57946805E-02 3.55032777E-03 7 -2.59391250E-02 -1.08404622E-02 8 -2.54714942E-02 -3.07632755E-02 9 3.96335778E-03 1.03505689E-01 10 5.39993660E-05 7.12373201E-01 ********************************************************* * ODRPACK95 version 1.00 of 12-27-2005 (REAL (KIND=R8)) * ********************************************************* *** Initial summary for fit by method of ODR *** --- Problem Size: N = 10 (number with nonzero weight = 10) NQ = 1 M = 1 NP = 2 (number unfixed = 2) --- Control Values: JOB = 00000 = ABCDE, where A=0 ==> fit is not a restart. B=0 ==> deltas are initialized to zero. C=0 ==> covariance matrix will be computed using derivatives re-evaluated at the solution. D=0 ==> derivatives are estimated by forward differences. E=0 ==> method is explicit ODR. NDIGIT = 16 (estimated by ODRPACK95) TAUFAC = 1.00E+00 --- Stopping Criteria: SSTOL = 1.49E-08 (sum of squares stopping tolerance) PARTOL = 3.67E-11 (parameter stopping tolerance) MAXIT = 1000 (maximum number of iterations) --- Initial Weighted Sum of Squares = 1.44889752E+02 Sum of Squared Weighted Deltas = 0.00000000E+00 Sum of Squared Weighted Epsilons = 1.44889752E+02 *** Iteration reports for fit by method of ODR *** Cum. Act. Rel. Pred. Rel. It. No. FN Weighted Sum-of-Sqs Sum-of-Sqs G-N Num. Evals Sum-of-Sqs Reduction Reduction TAU/PNORM Step ---- ------ ----------- ----------- ----------- --------- ---- 1 9 1.26593E+01 9.1263E-01 9.7623E-01 4.318E-01 YES 2 13 3.89042E+00 6.9268E-01 7.0679E-01 1.649E-01 YES 3 17 3.74830E+00 3.6530E-02 3.6740E-02 2.811E-02 YES 4 21 3.74825E+00 1.4256E-05 1.3893E-05 6.014E-04 YES 5 25 3.74825E+00 9.1626E-08 7.8661E-08 8.161E-05 YES 6 29 3.74825E+00 2.6737E-09 2.1587E-09 1.333E-05 YES *** Final summary for fit by method of ODR *** --- Stopping Conditions: INFO = 1 ==> sum of squares convergence. NITER = 6 (number of iterations) NFEV = 32 (number of function evaluations) IRANK = 0 (rank deficiency) RCOND = 1.25E-01 (inverse condition number) ISTOP = 0 (returned by user from subroutine FCN) --- Final Weighted Sums of Squares = 3.74825109E+00 Sum of Squared Weighted Deltas = 1.69631727E+00 Sum of Squared Weighted Epsilons = 2.05193383E+00 --- Residual Standard Deviation = 6.84493526E-01 Degrees of Freedom = 8 --- Estimated BETA(J), J = 1, ..., NP: BETA LOWER UPPER S.D. ___ 95% Confidence ___ BETA Interval 1 4.17524609E-01 -1.80+308 1.80+308 4.84E-02 3.06E-01 to 5.29E-01 2 1.00566711E+00 -1.80+308 1.80+308 3.45E-02 9.26E-01 to 1.09E+00 --- Estimated EPSILON(I) and DELTA(I,*), I = 1, ..., N: I EPSILON(I,1) DELTA(I,1) 1 6.01323538E-03 -1.36036478E-02 2 -1.27351622E-01 3.16511408E-03 3 1.56333747E-01 1.42422929E-03 4 -1.25881315E-01 -3.44792864E-03 5 5.25058480E-02 5.93961619E-03 6 1.91767312E-02 4.29866450E-03 7 -2.32510314E-02 -9.70473307E-03 8 -2.39929864E-02 -2.89350990E-02 9 4.13411546E-03 1.07612829E-01 10 5.52238407E-05 7.25473428E-01