sts_ttest.h

Go to the documentation of this file.

//=============================================================================================================
 
#ifndef STS_TTEST_H
#define STS_TTEST_H
 
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
#include "sts_global.h"
#include "sts_types.h"
 
//=============================================================================================================
// EIGEN INCLUDES
//=============================================================================================================
 
#include <Eigen/Core>
 
//=============================================================================================================
// DEFINE NAMESPACE STSLIB
//=============================================================================================================
 
namespace STSLIB
{
 
//=============================================================================================================
struct STSSHARED_EXPORT StatsTtestResult
{
    Eigen::MatrixXd matTstat; 
    Eigen::MatrixXd matPval;  
    int degreesOfFreedom = 0; 
};
 
//=============================================================================================================
class STSSHARED_EXPORT StatsTtest
{
public:
    //=========================================================================================================
    static StatsTtestResult oneSample(const Eigen::MatrixXd& data,
                                      double mu = 0.0,
                                      StatsTailType tail = StatsTailType::Both);
 
    //=========================================================================================================
    static StatsTtestResult paired(const Eigen::MatrixXd& dataA,
                                   const Eigen::MatrixXd& dataB,
                                   StatsTailType tail = StatsTailType::Both);
 
    //=========================================================================================================
    static StatsTtestResult independent(const Eigen::MatrixXd& dataA,
                                        const Eigen::MatrixXd& dataB,
                                        StatsTailType tail = StatsTailType::Both);
 
    //=========================================================================================================
    static StatsTtestResult permutation(const Eigen::MatrixXd& data,
                                        int nPermutations = 10000,
                                        StatsTailType tail = StatsTailType::Both,
                                        unsigned int seed = 0);
 
    //=========================================================================================================
    static double tCdf(double t, int df);
 
    //=========================================================================================================
    static double regularizedBeta(double x, double a, double b);
 
private:
    //=========================================================================================================
    static double betaCf(double x, double a, double b);
 
    //=========================================================================================================
    static double logBeta(double a, double b);
 
    //=========================================================================================================
    static double tToPval(double t, int df, StatsTailType tail);
};
 
} // namespace STSLIB
 
#endif // STS_TTEST_H