1 // leastsqs.h -- Implements a simple linear least squares best fit routine
3 // Written by Curtis Olson, started September 1997.
5 // Copyright (C) 1997 Curtis L. Olson - curt@infoplane.com
7 // This program is free software; you can redistribute it and/or modify
8 // it under the terms of the GNU General Public License as published by
9 // the Free Software Foundation; either version 2 of the License, or
10 // (at your option) any later version.
12 // This program is distributed in the hope that it will be useful,
13 // but WITHOUT ANY WARRANTY; without even the implied warranty of
14 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
15 // GNU General Public License for more details.
17 // You should have received a copy of the GNU General Public License
18 // along with this program; if not, write to the Free Software
19 // Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
22 // (Log is kept at end of this file)
31 # error This library requires C++
40 n*sum(xi*yi) - (sum(xi)*sum(yi))
41 b1 = --------------------------------
42 n*sum(xi^2) - (sum(xi))^2
45 b0 = sum(yi)/n - b1*(sum(xi)/n)
48 void least_squares(double *x, double *y, int n, double *m, double *b);
52 return the least squares error:
58 double least_squares_error(double *x, double *y, int n, double m, double b);
62 return the maximum least squares error:
66 double least_squares_max_error(double *x, double *y, int n, double m, double b);
73 // Revision 1.2 1998/04/21 17:03:42 curt
74 // Prepairing for C++ integration.
76 // Revision 1.1 1998/04/08 22:57:25 curt
77 // Adopted Gnu automake/autoconf system.
79 // Revision 1.1 1998/03/19 02:54:48 curt
80 // Reorganized into a class lib called fgDEM.
82 // Revision 1.1 1997/10/13 17:02:35 curt