nd observed responses, respectively.

 CALL: [BETA, SIGMA2] = stk_param_gls (MODEL, XI, ZI)

   also returns the associated unbiased estimate SIGMA2 of sigma^2, assu-
   ming that the actual covariance matrix of the Gaussian process part of
   the model is sigma^2 K, with K the covariance matrix built from MODEL.

   SIGMA2 is actually the "best" unbiased estimate of sigma^2 :

                 1
      SIGMA2 = ----- * || ZI - P BETA ||^2_{K^{-1}}
               n - r

   where n is the number of observatio