FIT OF THE INDOMETHACIN DATA BY GLS 
 
 
OLS estimate  1.271474 1.040768 -1.232717 -1.506841 
 
 
Iteration  1 
Estimate of theta  0.6928 
GLS estimate of beta  1.241283 0.972565 -1.432457 -1.747787 
 
Iteration  2 
Estimate of theta  0.8192 
GLS estimate of beta  1.236734 0.970361 -1.428257 -1.741751 
 
Iteration  3 
Estimate of theta  0.8188 
GLS estimate of beta  1.236594 0.970094 -1.428944 -1.74254 
 
Iteration  4 
Estimate of theta  0.8189 
GLS estimate of beta  1.236587 0.970098 -1.428899 -1.742481 
 
Iteration  5 
Estimate of theta  0.8189 
GLS estimate of beta  1.236587 0.970097 -1.428904 -1.742486 
 
Iteration  6 
Estimate of theta  0.8189 
GLS estimate of beta  1.236587 0.970097 -1.428904 -1.742486 
 
Iteration  7 
Estimate of theta  0.8189 
GLS estimate of beta  1.236587 0.970097 -1.428904 -1.742486 
 
Iteration  8 
Estimate of theta  0.8189 
GLS estimate of beta  1.236587 0.970097 -1.428904 -1.742486 
 
Iteration  9 
Estimate of theta  0.8189 
GLS estimate of beta  1.236587 0.970097 -1.428904 -1.742486 
 
Iteration  10 
Estimate of theta  0.8189 
GLS estimate of beta  1.236587 0.970097 -1.428904 -1.742486 
 
Final estimate of sigma  0.131344 
 

Formula: concwt ~ weightfunc(time, b1, b2, b3, b4, wt)

Parameters:
   Estimate Std. Error t value Pr(>|t|)    
b1   1.2366     0.1768   6.993 0.000213 ***
b2   0.9701     0.1379   7.034 0.000205 ***
b3  -1.4289     0.2322  -6.155 0.000465 ***
b4  -1.7425     0.2632  -6.621 0.000298 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

Residual standard error: 0.1313 on 7 degrees of freedom

Number of iterations to convergence: 5 
Achieved convergence tolerance: 9.459e-06