/usr/include/shogun/multiclass/RejectionStrategy.h is in libshogun-dev 3.2.0-7.3build4.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 | /*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 3 of the License, or
* (at your option) any later version.
*
* Written (W) 2012 Sergey Lisitsyn
* Copyright (C) 2012 Sergey Lisitsyn
*/
#ifndef _REJECTIONSTRATEGY_H___
#define _REJECTIONSTRATEGY_H___
namespace shogun
{
/** @brief base rejection strategy class */
class CRejectionStrategy : public CSGObject
{
public:
/** default constructor */
CRejectionStrategy() { };
/** destructor */
virtual ~CRejectionStrategy() { };
/** get name */
virtual const char* get_name() const
{
return "RejectionStrategy";
};
/** returns true if given output set leads to rejection */
virtual bool reject(SGVector<float64_t> outputs) const = 0;
};
/** @brief threshold based rejection strategy */
class CThresholdRejectionStrategy : public CRejectionStrategy
{
public:
/** constructor */
CThresholdRejectionStrategy() :
CRejectionStrategy(), m_threshold(0.0) { };
/** constructor */
CThresholdRejectionStrategy(float64_t threshold) :
CRejectionStrategy(), m_threshold(threshold) { };
virtual ~CThresholdRejectionStrategy() {};
/** get name */
virtual const char* get_name() const
{
return "ThresholdRejectionStrategy";
}
/** returns true if given output set leads to rejection */
virtual bool reject(SGVector<float64_t> outputs) const
{
for (int32_t i=0; i<outputs.vlen; i++)
{
if (outputs[i]>m_threshold)
return false;
}
return true;
}
protected:
/** threshold */
float64_t m_threshold;
};
static const float64_t Q_test_statistic_values[10][8] =
{
/* 10,20,30,40,50,60,70,80,90,100 */
{0.713,0.683,0.637,0.597,0.551,0.477,0.409,0.325},
{0.627,0.604,0.568,0.538,0.503,0.450,0.401,0.339},
{0.539,0.517,0.484,0.456,0.425,0.376,0.332,0.278},
{0.490,0.469,0.438,0.412,0.382,0.337,0.295,0.246},
{0.460,0.439,0.410,0.384,0.355,0.312,0.272,0.226},
{0.437,0.417,0.388,0.363,0.336,0.294,0.256,0.211},
{0.422,0.403,0.374,0.349,0.321,0.280,0.244,0.201},
{0.408,0.389,0.360,0.337,0.310,0.270,0.234,0.192},
{0.397,0.377,0.350,0.326,0.300,0.261,0.226,0.185},
{0.387,0.368,0.341,0.317,0.292,0.253,0.219,0.179}
};
/** @brief simplified version of Dixon's Q test outlier based
* rejection strategy. Statistic values are taken from
* http://www.vias.org/tmdatanaleng/cc_outlier_tests_dixon.html
* */
class CDixonQTestRejectionStrategy : public CRejectionStrategy
{
public:
/** constructor */
CDixonQTestRejectionStrategy() :
CRejectionStrategy()
{
s_index = 3;
}
/** constructor
* @param significance_level either 0.001,0.002,0.005,
* 0.01,0.02,0.05,0.1 or 0.2
*/
CDixonQTestRejectionStrategy(float64_t significance_level) :
CRejectionStrategy()
{
if (significance_level==0.001)
s_index = 0;
else if (significance_level==0.002)
s_index = 1;
else if (significance_level==0.005)
s_index = 2;
else if (significance_level==0.01)
s_index = 3;
else if (significance_level==0.02)
s_index = 4;
else if (significance_level==0.05)
s_index = 5;
else if (significance_level==0.1)
s_index = 6;
else if (significance_level==0.2)
s_index = 7;
else SG_ERROR("Given significance level is not supported")
}
virtual ~CDixonQTestRejectionStrategy()
{
}
/** get name */
virtual const char* get_name() const
{
return "DixonQTestRejectionStrategy";
}
/** returns true if given output set leads to rejection */
virtual bool reject(SGVector<float64_t> outputs) const
{
int32_t N = outputs.vlen;
if (N<10 || N>100)
SG_ERROR("Given number of classes is not supported.")
int32_t Ni = N/10 - 1;
SGVector<float64_t> outputs_local = outputs.clone();
outputs_local.qsort();
float64_t Q = 0.0;
if (N==10)
Q = (outputs[N-1]-outputs[N-2])/(outputs[N-1]-outputs[0]);
if (N>=20)
Q = (outputs[N-1]-outputs[N-4])/(outputs[N-1]-outputs[2]);
if (Q>Q_test_statistic_values[Ni][s_index])
return false;
return true;
}
private:
int32_t s_index;
};
}
#endif
|