Function rgsl::fit::wmul

pub fn wmul(
    x: &[f64], 
    xstride: usize, 
    w: &[f64], 
    wstride: usize, 
    y: &[f64], 
    ystride: usize, 
    n: usize, 
    c1: &mut f64, 
    cov11: &mut f64, 
    sumsq: &mut f64
) -> Value

This function computes the best-fit linear regression coefficient c1 of the model Y = c_1 X for the weighted datasets (x, y), two vectors of length n with strides xstride and ystride. The vector w, of length n and stride wstride, specifies the weight of each datapoint. The weight is the reciprocal of the variance for each datapoint in y.

The variance of the parameter c1 is computed using the weights and returned via the parameter cov11. The weighted sum of squares of the residuals from the best-fit line, \chi^2, is returned in chisq.