********************************************************* * 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 = 100 (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 = 11 (number with nonzero weight = 11) 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 = 100 (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 11 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 = 11 (number with nonzero weight = 11) 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 = 100 (maximum number of iterations) --- Initial Weighted Sum of Squares = 7.27357095E+00 Sum of Squared Weighted Deltas = 0.00000000E+00 Sum of Squared Weighted Epsilons = 7.27357095E+00 *** 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 4.55422E-01 9.3739E-01 9.4466E-01 2.314E-01 YES 2 11 4.28383E-01 5.9370E-02 6.0041E-02 2.134E-02 YES 3 14 4.28329E-01 1.2670E-04 1.6083E-04 2.294E-03 YES 4 17 4.28326E-01 5.8743E-06 7.6181E-06 5.046E-04 YES 5 20 4.28326E-01 3.0929E-07 4.0827E-07 1.155E-04 YES 6 23 4.28326E-01 1.7750E-08 2.2019E-08 2.678E-05 YES 7 26 4.28326E-01 7.5104E-10 1.0451E-09 5.822E-06 YES *** Final summary for fit by method of ODR *** --- Stopping Conditions: INFO = 1 ==> sum of squares convergence. NITER = 7 (number of iterations) NFEV = 28 (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 = 4.28326170E-01 Sum of Squared Weighted Deltas = 2.92820636E-01 Sum of Squared Weighted Epsilons = 1.35505534E-01 --- Residual Standard Deviation = 2.06960424E-01 Degrees of Freedom = 10 --- Estimated BETA(J), J = 1, ..., NP: BETA LOWER UPPER S.D. ___ 95% Confidence ___ BETA Interval 1 1.15120566E+04 -1.80+308 1.80+308 1.30E+02 1.12E+04 to 1.18E+04 --- Estimated EPSILON(I) and DELTA(I,*), I = 1, ..., N: I EPSILON(I,1) DELTA(I,1) 1 4.77474686E-02 5.06054161E-02 2 3.53962313E-01 6.28986617E-02 3 1.32100288E+00 3.54572906E-02 4 1.47861049E+00 4.47224638E-02 5 3.64650678E+00 1.52708979E-02 6 4.62053884E+00 2.01987437E-02 7 2.91140389E+00 1.27215101E-02 8 9.26627469E-01 3.60035723E-03 9 -7.06053668E+00 -1.21146653E-02 10 -1.21492054E+01 4.13664867E-02 11 -8.34434473E-01 2.82802633E-02 ********************************************************* * ODRPACK95 version 1.00 of 12-27-2005 (REAL (KIND=R8)) * ********************************************************* *** Initial summary for fit by method of ODR *** --- Problem Size: N = 11 (number with nonzero weight = 11) 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 = 100 (maximum number of iterations) --- Initial Weighted Sum of Squares = 4.98191074E+00 Sum of Squared Weighted Deltas = 0.00000000E+00 Sum of Squared Weighted Epsilons = 4.98191074E+00 *** 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 6.47497E-02 9.8700E-01 9.9620E-01 1.867E-01 YES 2 13 1.43667E-02 7.7812E-01 7.7838E-01 1.382E-02 YES 3 17 1.43658E-02 5.8242E-05 5.8256E-05 1.056E-04 YES 4 21 1.43658E-02 9.4061E-11 8.7794E-11 1.742E-07 YES *** Final summary for fit by method of ODR *** --- Stopping Conditions: INFO = 1 ==> sum of squares convergence. NITER = 4 (number of iterations) NFEV = 24 (number of function evaluations) IRANK = 0 (rank deficiency) RCOND = 9.09E-02 (inverse condition number) ISTOP = 0 (returned by user from subroutine FCN) --- Final Weighted Sums of Squares = 1.43658280E-02 Sum of Squared Weighted Deltas = 1.24033960E-02 Sum of Squared Weighted Epsilons = 1.96243207E-03 --- Residual Standard Deviation = 3.99525107E-02 Degrees of Freedom = 9 --- Estimated BETA(J), J = 1, ..., NP: BETA LOWER UPPER S.D. ___ 95% Confidence ___ BETA Interval 1 1.02931489E+04 -1.80+308 1.80+308 7.59E+01 1.01E+04 to 1.05E+04 2 2.91697314E+00 -1.80+308 1.80+308 7.02E-03 2.90E+00 to 2.93E+00 --- Estimated EPSILON(I) and DELTA(I,*), I = 1, ..., N: I EPSILON(I,1) DELTA(I,1) 1 -5.33296258E-02 -5.76070992E-02 2 -3.60321968E-02 -6.50654366E-03 3 -2.15187489E-01 -5.80953023E-03 4 2.15836674E-01 6.63949749E-03 5 -6.12163837E-01 -2.58104111E-03 6 8.31668356E-01 3.70160613E-03 7 5.09150057E-02 2.29969450E-04 8 -1.32534505E-01 -5.50702527E-04 9 6.83763280E-01 1.67384983E-03 10 1.00346421E-02 -4.46431364E-05 11 -2.70063290E-01 9.81439552E-03