00001 //# Copyright (C) 2000,2001 00002 //# Associated Universities, Inc. Washington DC, USA. 00003 //# 00004 //# This library is free software; you can redistribute it and/or modify it 00005 //# under the terms of the GNU Library General Public License as published by 00006 //# the Free Software Foundation; either version 2 of the License, or (at your 00007 //# option) any later version. 00008 //# 00009 //# This library is distributed in the hope that it will be useful, but WITHOUT 00010 //# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or 00011 //# FITNESS FOR A PARTICULAR PURPOSE. See the GNU Library General Public 00012 //# License for more details. 00013 //# 00014 //# You should have received a copy of the GNU Library General Public License 00015 //# along with this library; if not, write to the Free Software Foundation, 00016 //# Inc., 675 Massachusetts Ave, Cambridge, MA 02139, USA. 00017 //# 00018 //# Correspondence concerning AIPS++ should be addressed as follows: 00019 //# Internet email: aips2-request@nrao.edu. 00020 //# Postal address: AIPS++ Project Office 00021 //# National Radio Astronomy Observatory 00022 //# 520 Edgemont Road 00023 //# Charlottesville, VA 22903-2475 USA 00024 //# 00025 //# $Id: HostInfoDarwin.h 21521 2014-12-10 08:06:42Z gervandiepen $ 00026 00027 #ifndef SCIMATH_ZSCORECALCULATOR_H 00028 #define SCIMATH_ZSCORECALCULATOR_H 00029 00030 #include <casacore/casa/aips.h> 00031 00032 #include <map> 00033 #include <set> 00034 #include <math.h> 00035 00036 namespace casacore { 00037 00038 /* 00039 * This class contains static methods related to z-scores. A z-score is the number of 00040 * standard deviations from the mean in a normal distribution. 00041 */ 00042 00043 class ZScoreCalculator { 00044 public: 00045 00046 // compute the maximum expected zscore given the number of points 00047 // in a sample. 00048 static Double getMaxZScore(uInt64 npts); 00049 00050 // Get the minimum number of points in a Gaussian distribution, such that the 00051 // probability that the maximum of the distribution having the specified zscore 00052 // is 0.5. <src>zscore</src> should be non-negative. 00053 static inline uInt64 zscoreToNpts(Double zscore) { 00054 return (uInt64)(0.5/erfc(zscore/sqrt(2))); 00055 } 00056 00057 private: 00058 static std::map<uInt64, Double> _nptsToMaxZScore; 00059 00060 }; 00061 00062 } 00063 00064 #endif