PLearn 0.1
KolmogorovSmirnovCommand.cc
Go to the documentation of this file.
00001 // -*- C++ -*-
00002 
00003 // KolmogorovSmirnovCommand.cc
00004 // 
00005 // Copyright (C) 2003 Pascal Vincent
00006 // 
00007 // Redistribution and use in source and binary forms, with or without
00008 // modification, are permitted provided that the following conditions are met:
00009 // 
00010 //  1. Redistributions of source code must retain the above copyright
00011 //     notice, this list of conditions and the following disclaimer.
00012 // 
00013 //  2. Redistributions in binary form must reproduce the above copyright
00014 //     notice, this list of conditions and the following disclaimer in the
00015 //     documentation and/or other materials provided with the distribution.
00016 // 
00017 //  3. The name of the authors may not be used to endorse or promote
00018 //     products derived from this software without specific prior written
00019 //     permission.
00020 // 
00021 // THIS SOFTWARE IS PROVIDED BY THE AUTHORS ``AS IS'' AND ANY EXPRESS OR
00022 // IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
00023 // OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN
00024 // NO EVENT SHALL THE AUTHORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
00025 // SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED
00026 // TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
00027 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
00028 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
00029 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
00030 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
00031 // 
00032 // This file is part of the PLearn library. For more information on the PLearn
00033 // library, go to the PLearn Web site at www.plearn.org
00034 
00035 /* *******************************************************      
00036  * $Id: KolmogorovSmirnovCommand.cc 3995 2005-08-25 13:58:23Z chapados $ 
00037  ******************************************************* */
00038 
00040 #include "KolmogorovSmirnovCommand.h"
00041 #include <plearn/db/getDataSet.h>
00042 #include <plearn/math/stats_utils.h>
00043 
00044 namespace PLearn {
00045 using namespace std;
00046 
00048 PLearnCommandRegistry KolmogorovSmirnovCommand::reg_(new KolmogorovSmirnovCommand);
00049 
00051 void KolmogorovSmirnovCommand::run(const vector<string>& args)
00052 {
00053     if(args.size()!=4 && args.size()!=5)
00054         PLERROR("ks-stat expects 4 or 5 arguments, check the help");
00055 
00056     VMat m1 = getDataSet(args[0]);
00057     int c1 = toint(args[1]);
00058     VMat m2 = getDataSet(args[2]);
00059     int c2 = toint(args[3]);
00060     int conv = 10;
00061     if(args.size()==5)
00062         conv = toint(args[4]);
00063 
00064     Vec v1 = m1.getColumn(c1);
00065     Vec v2 = m2.getColumn(c2);
00066 
00067     real D, ks_stat;
00068     KS_test(v1, v2, conv, D, ks_stat);
00069 
00070     cout << "Maximum absolute difference between the 2 cdfs is: " << D << endl;
00071     cout << "Result of Kolmogorov-Smirnov test is: " << ks_stat << endl;
00072 }
00073 
00074 } // end of namespace PLearn
00075 
00076 
00077 /*
00078   Local Variables:
00079   mode:c++
00080   c-basic-offset:4
00081   c-file-style:"stroustrup"
00082   c-file-offsets:((innamespace . 0)(inline-open . 0))
00083   indent-tabs-mode:nil
00084   fill-column:79
00085   End:
00086 */
00087 // vim: filetype=cpp:expandtab:shiftwidth=4:tabstop=8:softtabstop=4:encoding=utf-8:textwidth=79 :
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines