00001 /*************************************************************************** 00002 ** 00003 ** This file is part of QGpCoreTools. 00004 ** 00005 ** This library is free software; you can redistribute it and/or 00006 ** modify it under the terms of the GNU Lesser General Public 00007 ** License as published by the Free Software Foundation; either 00008 ** version 2.1 of the License, or (at your option) any later version. 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 Lesser General Public 00013 ** License for more details. 00014 ** 00015 ** You should have received a copy of the GNU Lesser General Public 00016 ** License along with this library; if not, write to the Free Software 00017 ** Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA 00018 ** 00019 ** See http://www.geopsy.org for more information. 00020 ** 00021 ** Created : 2003-12-15 00022 ** Authors : 00023 ** Marc Wathelet 00024 ** Marc Wathelet (ULg, Liège, Belgium) 00025 ** Marc Wathelet (LGIT, Grenoble, France) 00026 ** 00027 ***************************************************************************/ 00028 00029 #ifndef StatComparator_H 00030 #define StatComparator_H 00031 00032 #include "QGpCoreToolsDLLExport.h" 00033 #include "Trace.h" 00034 00035 namespace QGpCoreTools { 00036 00037 class QGPCORETOOLS_EXPORT StatComparator 00038 { 00039 public: 00040 StatComparator() {} 00041 ~StatComparator() {} 00042 00043 static inline double student(double t, double nu); 00044 static inline double fisher(double F, double nu1, double nu2); 00045 static double invStudent(double p, double nu); 00046 static double invFisher(double p, double nu1, double nu2); 00047 protected: 00048 static double betai(double a, double b, double x); 00049 static double gammln(double x); 00050 static double betacf(double a, double b, double x); 00051 }; 00052 00058 inline double StatComparator::student(double t, double nu) 00059 { 00060 TRACE; 00061 return betai(nu/2.0, 0.5, nu/(nu + t * t)); 00062 } 00063 00069 inline double StatComparator::fisher(double F, double nu1, double nu2) 00070 { 00071 TRACE; 00072 return betai(nu2/2.0, nu1/2.0, nu2/(nu2 + nu1 * F)); 00073 } 00074 00075 } // namespace QGpCoreTools 00076 00077 #endif // STATCOMPARATOR_H
