void linreg(n,x,y,a,siga,b,sigb)int n;double x[n], y[n];double a[6], siga[6], b[6], sigb[6];
The dispersions for OLS are the general form provided by T. Isobe, E. D. Feigelson, M. G. Akritas and G. J. Babu, Ap.J. 364, 104.
mean OLS is 0.5*[OLS(X/Y) + OLS(Y/X)]
~/src/kernel/math linreg.c
xx-jul-90 Created Isobe 13-aug-90 Ported into NEMO PJT 28-oct-90 upgraded from Feigelstein’s official rel. 1.0 PJT