ARGATROBAN PK DATA 
 
ML fit with estimation of theta 
 

**Iteration 1
LME step: Loglik: -2911.289 , nlm iterations: 34 
reStruct  parameters:
     indiv1      indiv2      indiv3 
   5.868500    6.083282 -548.297866 
varStruct  parameters:
     power 
0.07704409 

PNLS step: RSS =  1339347 
 fixed effects:-5.43331  -1.98126  
 iterations: 7 

Convergence:
     fixed   reStruct  varStruct 
 0.1042991 13.6429449  5.4897905 

**Iteration 2
LME step: Loglik: -2865.409 , nlm iterations: 25 
reStruct  parameters:
    indiv1     indiv2     indiv3 
  4.270221   5.964866 -61.487347 
varStruct  parameters:
    power 
0.2029238 

PNLS step: RSS =  325751.6 
 fixed effects:-5.43269  -1.92542  
 iterations: 6 

Convergence:
     fixed   reStruct  varStruct 
0.02899958 1.76831917 0.62033002 

**Iteration 3
LME step: Loglik: -2862.572 , nlm iterations: 11 
reStruct  parameters:
    indiv1     indiv2     indiv3 
  4.074090   5.758618 -89.603591 
varStruct  parameters:
    power 
0.2385338 

PNLS step: RSS =  210084.3 
 fixed effects:-5.43257  -1.91873  
 iterations: 4 

Convergence:
      fixed    reStruct   varStruct 
0.003485758 0.057633632 0.149286759 

**Iteration 4
LME step: Loglik: -2862.257 , nlm iterations: 6 
reStruct  parameters:
    indiv1     indiv2     indiv3 
  4.047004   5.717158 -84.722059 
varStruct  parameters:
    power 
0.2427381 

PNLS step: RSS =  199416.2 
 fixed effects:-5.43255  -1.91813  
 iterations: 2 

Convergence:
       fixed     reStruct    varStruct 
0.0003164304 0.0062071673 0.0173203627 

**Iteration 5
LME step: Loglik: -2862.221 , nlm iterations: 7 
reStruct  parameters:
    indiv1     indiv2     indiv3 
  4.044165   5.712790 -84.199422 
varStruct  parameters:
    power 
0.2431754 

PNLS step: RSS =  198334.8 
 fixed effects:-5.43255  -1.91795  
 iterations: 2 

Convergence:
       fixed     reStruct    varStruct 
0.0000910878 0.0010075849 0.0017982779 

**Iteration 6
LME step: Loglik: -2862.212 , nlm iterations: 4 
reStruct  parameters:
    indiv1     indiv2     indiv3 
  4.044138   5.712324 -84.114745 
varStruct  parameters:
    power 
0.2432005 

PNLS step: RSS =  198260.1 
 fixed effects:-5.43255  -1.91795  
 iterations: 1 

Convergence:
       fixed     reStruct    varStruct 
0.0000000000 0.0004478655 0.0001034273 

**Iteration 7
LME step: Loglik: -2862.212 , nlm iterations: 7 
reStruct  parameters:
    indiv1     indiv2     indiv3 
  4.043925   5.712194 -84.077090 
varStruct  parameters:
    power 
0.2432239 

PNLS step: RSS =  198205.8 
 fixed effects:-5.43255  -1.91795  
 iterations: 1 

Convergence:
       fixed     reStruct    varStruct 
0.000000e+00 3.126598e-05 9.626624e-05 

**Iteration 8
LME step: Loglik: -2862.212 , nlm iterations: 9 
reStruct  parameters:
    indiv1     indiv2     indiv3 
  4.043798   5.712062 -84.079602 
varStruct  parameters:
    power 
0.2432412 

PNLS step: RSS =  198165.8 
 fixed effects:-5.43255  -1.91795  
 iterations: 1 

Convergence:
       fixed     reStruct    varStruct 
0.000000e+00 2.305688e-05 7.092974e-05 

**Iteration 9
LME step: Loglik: -2862.212 , nlm iterations: 9 
reStruct  parameters:
    indiv1     indiv2     indiv3 
  4.043705   5.711955 -84.080817 
varStruct  parameters:
    power 
0.2432541 

PNLS step: RSS =  198135.9 
 fixed effects:-5.43255  -1.91795  
 iterations: 1 

Convergence:
       fixed     reStruct    varStruct 
0.000000e+00 1.727592e-05 5.299145e-05 

**Iteration 10
LME step: Loglik: -2862.212 , nlm iterations: 5 
reStruct  parameters:
    indiv1     indiv2     indiv3 
  4.043478   5.711703 -84.081229 
varStruct  parameters:
    power 
0.2432857 

PNLS step: RSS =  198062.7 
 fixed effects:-5.43255  -1.91795  
 iterations: 1 

Convergence:
       fixed     reStruct    varStruct 
0.000000e+00 4.225415e-05 1.297261e-04 

**Iteration 11
LME step: Loglik: -2862.212 , nlm iterations: 2 
reStruct  parameters:
    indiv1     indiv2     indiv3 
  4.043463   5.711668 -84.082725 
varStruct  parameters:
    power 
0.2432873 

PNLS step: RSS =  198058.8 
 fixed effects:-5.43255  -1.91795  
 iterations: 1 

Convergence:
       fixed     reStruct    varStruct 
0.000000e+00 1.308211e-05 6.771862e-06 

**Iteration 12
LME step: Loglik: -2862.212 , nlm iterations: 2 
reStruct  parameters:
    indiv1     indiv2     indiv3 
  4.043454   5.711647 -84.081625 
varStruct  parameters:
    power 
0.2432885 

PNLS step: RSS =  198055.9 
 fixed effects:-5.43255  -1.91795  
 iterations: 1 

Convergence:
       fixed     reStruct    varStruct 
0.000000e+00 7.524887e-06 5.036608e-06 

 
 
Nonlinear mixed-effects model fit by maximum likelihood
  Model: ys ~ meanfunc(xs, b1, b2, ds) 
 Data: sdat 
       AIC      BIC    logLik
  5738.425 5767.568 -2862.212

Random effects:
 Formula: list(b1 ~ 1, b2 ~ 1)
 Level: indiv
 Structure: General positive-definite, Log-Cholesky parametrization
         StdDev      Corr 
b1        0.37168388 b1   
b2        0.06753306 0.268
Residual 20.41959281      

Variance function:
 Structure: Power of variance covariate
 Formula: ~fitted(.) 
 Parameter estimates:
    power 
0.2432885 
Fixed effects: list(b1 ~ 1, b2 ~ 1) 
       Value  Std.Error  DF   t-value p-value
b1 -5.432545 0.06230331 437 -87.19513       0
b2 -1.917953 0.02513017 437 -76.32073       0
 Correlation: 
   b1   
b2 0.156

Standardized Within-Group Residuals:
        Min          Q1         Med          Q3         Max 
-2.48366029 -0.42424578  0.03283586  0.57935560  9.38411265 

Number of Observations: 475
Number of Groups: 37 
Estimate of sigma 20.41959