

Trace output for unweighted fit of Tyuma data to MWC model
0.9212432 :  3e+05 1e+01 1e-01 
0.920535 :  2.746657e+05 1.001420e+01 1.022784e-01 
0.9196883 :  2.520595e+05 1.002836e+01 1.045538e-01 
0.9187126 :  2.318338e+05 1.004247e+01 1.068257e-01 
0.9176167 :  2.136910e+05 1.005654e+01 1.090939e-01 
0.9164085 :  1.973768e+05 1.007055e+01 1.113578e-01 
0.9150957 :  1.826716e+05 1.008451e+01 1.136173e-01 
0.913685 :  1.693866e+05 1.009841e+01 1.158719e-01 
0.912183 :  1.573580e+05 1.011225e+01 1.181214e-01 
0.9105958 :  1.464437e+05 1.012603e+01 1.203655e-01 
0.9089288 :  1.365204e+05 1.013975e+01 1.226040e-01 
0.9071873 :  1.274800e+05 1.015340e+01 1.248364e-01 
0.9069935 :  1.109767e+05 1.018056e+01 1.292890e-01 
0.9062772 :  9.717435e+04 1.020758e+01 1.337358e-01 
0.9050893 :  8.555802e+04 1.023442e+01 1.381734e-01 
0.9034757 :  7.572291e+04 1.026107e+01 1.425992e-01 
0.901478 :  6.734859e+04 1.028750e+01 1.470105e-01 
0.8991337 :  6.017969e+04 1.031370e+01 1.514049e-01 
0.8964765 :  5.401128e+04 1.033964e+01 1.557803e-01 
0.8935371 :  4.867798e+04 1.036530e+01 1.601346e-01 
0.8903432 :  4.404546e+04 1.039068e+01 1.644662e-01 
0.8869196 :  4.000403e+04 1.041576e+01 1.687732e-01 
0.886309 :  3.292313e+04 1.046529e+01 1.773351e-01 
0.8842008 :  2.744011e+04 1.051402e+01 1.858655e-01 
0.8807987 :  2.314035e+04 1.056179e+01 1.943458e-01 
0.8762832 :  1.972831e+04 1.060848e+01 2.027592e-01 
0.8708134 :  1.699058e+04 1.065396e+01 2.110916e-01 
0.8645286 :  1.477114e+04 1.069813e+01 2.193299e-01 
0.864196 :  1.113790e+04 1.078368e+01 2.355964e-01 
0.8592174 :  8649.3113143   10.8652466    0.2517340 
0.8504692 :  6900.5965208   10.9418562    0.2676208 
0.838732 :  5642.3190501   11.0127633    0.2831550 
0.8368308 :  3792.2214049   11.1421226    0.3133491 
0.8210401 :  2736.3395818   11.2507221    0.3428412 
0.7952516 :  2106.1371987   11.3337888    0.3708535 
0.7664592 :  1325.2987549   11.4454241    0.4227199 
0.7210365 :  666.0564523  11.4701195   0.5140121 
0.4930507 :  438.0523477  10.7861059   0.6476654 
0.3832181 :  1258.4566128    8.8469888    0.5245916 
0.3666529 :  1809.4626065    8.9143125    0.4837355 
0.3512292 :  2595.6576142    9.0220805    0.4462004 
0.3363632 :  3704.1678461    9.1636027    0.4121669 
0.3216882 :  5247.7606678    9.3339066    0.3815696 
0.3069969 :  7369.5954331    9.5292558    0.3542057 
0.2921959 :  10247.852917     9.746805     0.329807 
0.2772704 :  1.409984e+04 9.984350e+00 3.080836e-01 
0.2756696 :  2.427002e+04 1.049596e+01 2.694176e-01 
0.2509090 :  3.801872e+04 1.103888e+01 2.460407e-01 
0.2368678 :  7.597743e+04 1.221490e+01 2.103398e-01 
0.1138779 :  1.501072e+05 1.470782e+01 2.015080e-01 
0.09469324 :  2.471624e+05 1.623903e+01 1.827920e-01 
0.05070876 :  4.406931e+05 1.931413e+01 1.697299e-01 
0.005523758 :  8.012726e+05 2.525724e+01 1.652797e-01 
0.001875109 :  1.193920e+06 2.877658e+01 1.540837e-01 
0.000555952 :  1.319724e+06 2.932066e+01 1.538629e-01 
0.000555296 :  1.376899e+06 2.931119e+01 1.522332e-01 
0.0005551089 :  1.468793e+06 2.929652e+01 1.497692e-01 
0.0005522296 :  1.573780e+06 2.928074e+01 1.472670e-01 
0.000550839 :  1.585761e+06 2.927717e+01 1.470958e-01 
0.0005508384 :  1.592925e+06 2.927396e+01 1.469292e-01 
0.0005508383 :  1.593138e+06 2.927386e+01 1.469247e-01 


Trace output for weighted fit of Tyuma data to MWC model
82.8735 :  3e+05 1e+01 1e-01 
80.114 :  1.459458e+05 1.066133e+01 1.203105e-01 
70.1819 :  9.272611e+04 1.130205e+01 1.400593e-01 
50.43831 :  5.812285e+04 1.259158e+01 1.727625e-01 
19.92077 :  6.878221e+04 1.531741e+01 2.090566e-01 
10.13172 :  2.281734e+05 2.220330e+01 1.688506e-01 
0.7534203 :  2.977821e+05 2.824413e+01 2.098295e-01 
0.4307134 :  4.511020e+05 3.106440e+01 1.904118e-01 
0.2223312 :  7.332111e+05 3.474013e+01 1.732287e-01 
0.07319022 :  8.281407e+05 3.530510e+01 1.745496e-01 
0.07253278 :  8.802077e+05 3.537409e+01 1.723081e-01 
0.07249 :  8.846200e+05 3.537140e+01 1.722023e-01 
0.07248993 :  8.852499e+05 3.537075e+01 1.721759e-01 
0.07248993 :  8.852778e+05 3.537072e+01 1.721748e-01 


Trace output for weighted fit of Tyuma data to Adair model
0.2585905 :  0.0240 0.0743 0.0858 7.3700 
0.04729995 :  0.02426245 0.07630348 0.09108982 7.76367763 
0.04652856 :  0.02416987 0.07702789 0.08968054 7.89120788 
0.04652726 :  0.02415977 0.07718744 0.08942258 7.90412260 
0.04652726 :  0.02415899 0.07719909 0.08940569 7.90484162 
0.04652726 :  0.02415893 0.07719989 0.08940451 7.90488997 


Summary for unweighted fit of Tyuma data to MWC model

Formula: y ~ (L/Kt * p * (1 + p/Kt)^3 + p/Kr * (1 + p/Kr)^3)/(L * (1 + 
    p/Kt)^4 + (1 + p/Kr)^4)

Parameters:
    Estimate Std. Error t value Pr(>|t|)    
L  1.593e+06  2.340e+06   0.681   0.5052    
Kt 2.927e+01  1.659e+00  17.649 2.29e-12 ***
Kr 1.469e-01  5.479e-02   2.682   0.0158 *  
---
Signif. codes:  0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 

Residual standard error: 0.005692 on 17 degrees of freedom

Correlation of Parameter Estimates:
         L     Kt
Kt -0.7759       
Kr -0.9999 0.7808



Summary for weighted fit of Tyuma data to MWC model

Formula: rep(0, length(y)) ~ wt.mwc(y, p, L, Kt, Kr)

Parameters:
    Estimate Std. Error t value Pr(>|t|)    
L  8.853e+05  2.560e+05   3.458  0.00301 ** 
Kt 3.537e+01  1.137e+00  31.104  < 2e-16 ***
Kr 1.722e-01  1.209e-02  14.241 7.03e-11 ***
---
Signif. codes:  0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 

Residual standard error: 0.0653 on 17 degrees of freedom

Correlation of Parameter Estimates:
         L     Kt
Kt -0.5254       
Kr -0.9949 0.5716



Summary for weighted fit of Tyuma data to Adair model

Formula: rep(0, length(y)) ~ wt.adair(y, p, K1, K2, K3, K4)

Parameters:
   Estimate Std. Error t value Pr(>|t|)    
K1 0.024159   0.001549  15.594 4.26e-11 ***
K2 0.077200   0.018607   4.149 0.000755 ***
K3 0.089405   0.031573   2.832 0.012029 *  
K4 7.904890   1.425783   5.544 4.44e-05 ***
---
Signif. codes:  0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 

Residual standard error: 0.05393 on 16 degrees of freedom

Correlation of Parameter Estimates:
        K1      K2      K3
K2 -0.8937                
K3  0.7821 -0.9562        
K4 -0.7134  0.8980 -0.9843





Unweighted variance = sum((y - ycalc)^2)/df, 
for unweighted and weighted fits

         mwc       mwc.wt        ad.wt 
3.240225e-05 7.070510e-05 5.309127e-05 



Weighted variance = sum(((y - ycalc)/(y*(1-y)))^2)/df, 
for unweighted and weighted fits

        mwc      mwc.wt       ad.wt 
0.016560060 0.004264114 0.002907954 
