Method: CArray#fit_wlinear

Defined in:
lib/carray-gsl/core.rb

#fit_wlinear(x, w) ⇒ Object



259
260
261
262
263
264
265
266
267
268
# File 'lib/carray-gsl/core.rb', line 259

def fit_wlinear (x, w)
  CArray.attach(self, x, w) {
    c0, c1, c00, c01, c11, chi2, status = 
                       *GSL::Fit.wlinear(x.gv, w.gv, self.gv)
    dof   = self.elements - 2
    covar = CA_DOUBLE([[c00, c01],[c01, c11]])
    err   = Array.new(2){|i| Math::sqrt(chi2/dof*covar[i,i]) }
    return [c0, c1], err, chi2, dof, covar
  }
end