sts_correction.h

Go to the documentation of this file.

//=============================================================================================================
 
#ifndef STS_CORRECTION_H
#define STS_CORRECTION_H
 
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
#include "sts_global.h"
 
//=============================================================================================================
// EIGEN INCLUDES
//=============================================================================================================
 
#include <Eigen/Core>
 
//=============================================================================================================
// DEFINE NAMESPACE STSLIB
//=============================================================================================================
 
namespace STSLIB
{
 
//=============================================================================================================
class STSSHARED_EXPORT StatsMcCorrection
{
public:
    //=========================================================================================================
    using RejectMask = Eigen::Array<bool, Eigen::Dynamic, Eigen::Dynamic>;

    enum class FdrMethod
    {
        Independent,         
        NegativelyCorrelated 
    };
 
    //=========================================================================================================
    static Eigen::MatrixXd bonferroni(const Eigen::MatrixXd& pValues, double alpha = 0.05, RejectMask* reject = nullptr);
 
    //=========================================================================================================
    static Eigen::MatrixXd holmBonferroni(const Eigen::MatrixXd& pValues);
 
    //=========================================================================================================
    static Eigen::MatrixXd fdr(const Eigen::MatrixXd& pValues,
                               double alpha = 0.05,
                               FdrMethod method = FdrMethod::Independent,
                               RejectMask* reject = nullptr);
};
 
} // namespace STSLIB
 
#endif // STS_CORRECTION_H