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
|
/*=========================================================================
Program: Visualization Toolkit
Module: vtkContingencyStatistics.h
Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
All rights reserved.
See Copyright.txt or http://www.kitware.com/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notice for more information.
=========================================================================*/
/*-------------------------------------------------------------------------
Copyright 2011 Sandia Corporation.
Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
the U.S. Government retains certain rights in this software.
-------------------------------------------------------------------------*/
// .NAME vtkContingencyStatistics - A class for bivariate correlation contigency
// tables, conditional probabilities, and information entropy
//
// .SECTION Description
// Given a pair of columns of interest, this class provides the
// following functionalities, depending on the operation in which it is executed:
// * Learn: calculate contigency tables and corresponding discrete joint
// probability distribution.
// * Derive: calculate conditional probabilities, information entropies, and
// pointwise mutual information.
// * Assess: given two columns of interest with the same number of entries as
// input in port INPUT_DATA, and a corresponding bivariate probability distribution,
// * Test: calculate Chi-square independence statistic and, if VTK to R interface is available,
// retrieve corresponding p-value for independence testing.
//
// .SECTION Thanks
// Thanks to Philippe Pebay and David Thompson from Sandia National Laboratories
// for implementing this class.
// Updated by Philippe Pebay, Kitware SAS 2012
#ifndef vtkContingencyStatistics_h
#define vtkContingencyStatistics_h
#include "vtkFiltersStatisticsModule.h" // For export macro
#include "vtkStatisticsAlgorithm.h"
class vtkMultiBlockDataSet;
class vtkStringArray;
class vtkTable;
class vtkVariant;
class vtkIdTypeArray;
class vtkDoubleArray;
class VTKFILTERSSTATISTICS_EXPORT vtkContingencyStatistics : public vtkStatisticsAlgorithm
{
public:
vtkTypeMacro(vtkContingencyStatistics, vtkStatisticsAlgorithm);
void PrintSelf(ostream& os, vtkIndent indent);
static vtkContingencyStatistics* New();
// Description:
// Given a collection of models, calculate aggregate model
// NB: not implemented
virtual void Aggregate( vtkDataObjectCollection*,
vtkMultiBlockDataSet* ) { return; };
protected:
vtkContingencyStatistics();
~vtkContingencyStatistics();
// Description:
// Execute the calculations required by the Learn option.
virtual void Learn( vtkTable*,
vtkTable*,
vtkMultiBlockDataSet* );
// Description:
// Execute the calculations required by the Derive option.
virtual void Derive( vtkMultiBlockDataSet* );
// Description:
// Execute the calculations required by the Test option.
virtual void Test( vtkTable*,
vtkMultiBlockDataSet*,
vtkTable* );
// Description:
// Execute the calculations required by the Assess option.
virtual void Assess( vtkTable*,
vtkMultiBlockDataSet*,
vtkTable* );
//BTX
// Description:
// Calculate p-value. This will be overridden using the object factory with an
// R implementation if R is present.
virtual void CalculatePValues(vtkTable*);
//BTX
// Description:
// Provide the appropriate assessment functor.
// This one does nothing because the API is not sufficient for tables indexed
// by a separate summary table.
virtual void SelectAssessFunctor( vtkTable* outData,
vtkDataObject* inMeta,
vtkStringArray* rowNames,
AssessFunctor*& dfunc );
// Description:
// Provide the appropriate assessment functor.
// This one is the one that is actually used.
virtual void SelectAssessFunctor( vtkTable* outData,
vtkMultiBlockDataSet* inMeta,
vtkIdType pairKey,
vtkStringArray* rowNames,
AssessFunctor*& dfunc );
//ETX
private:
vtkContingencyStatistics(const vtkContingencyStatistics&); // Not implemented
void operator=(const vtkContingencyStatistics&); // Not implemented
};
#endif
|
