00001 /*************************************************************************** 00002 ** 00003 ** This file is part of GpCoreTools. 00004 ** 00005 ** This file may be distributed and/or modified under the terms of the 00006 ** GNU General Public License version 2 or 3 as published by the Free 00007 ** Software Foundation and appearing in the file LICENSE.GPL included 00008 ** in the packaging of this file. 00009 ** 00010 ** This file is distributed in the hope that it will be useful, but WITHOUT 00011 ** ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or 00012 ** FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for 00013 ** more details. 00014 ** 00015 ** You should have received a copy of the GNU General Public License 00016 ** along with this program. If not, see <http://www.gnu.org/licenses/>. 00017 ** 00018 ** See http://www.geopsy.org for more information. 00019 ** 00020 ** Created: 2011-04-21 00021 ** Authors: 00022 ** Marc Wathelet (LGIT, Grenoble, France) 00023 ** 00024 ***************************************************************************/ 00025 00026 #ifndef LINEARREGRESSION_H 00027 #define LINEARREGRESSION_H 00028 00029 #include <vector> 00030 00031 #include "GpCoreToolsDLLExport.h" 00032 00033 namespace GpCoreTools { 00034 00035 class GPCORETOOLS_EXPORT LinearRegression 00036 { 00037 public: 00038 LinearRegression() {reset();} 00039 00040 void add(double x, double y); 00041 void reset(); 00042 00043 bool calculate(); 00044 float stddev() const; 00045 int badValue(double deviation) const; 00046 void remove(int i); 00047 00048 00049 unsigned int count() const {return _values.size();} 00050 double a() const {return _a;} 00051 double b() const {return _b;} 00052 double at(double x) const {return _a*(x-_x0)+_b;} 00053 private: 00054 double _sumX, _sumY, _sumXY, _sumX2, _x0; 00055 double _a, _b; 00056 std::vector<std::pair<double, double> > _values; 00057 }; 00058 00059 } // namespace GpCoreTools 00060 00061 #endif // LINEARREGRESSION_H
