libsequence
1.9.5
|
results of calculations of the HKA test More...
#include <Sequence/HKA.hpp>
Public Types | |
enum | chisq_tuple_element : std::int8_t { POLY, DIV, POLYA, POLYB } |
typedef std::tuple< double, double, double, double > | chisq_tuple |
Public Member Functions | |
HKAresults (const std::vector< double > &_thetas, const std::vector< chisq_tuple > &_chisquareds, const double &_fhat, const double &_That, const double &_xsq, const double &_xsqA, const double &_xsqB) | |
Public Attributes | |
const std::vector< double > | thetas |
std::vector< chisq_tuple > | chisquareds |
const double | fhat |
const double | That |
const double | xsq |
const double | xsqA |
const double | xsqB |
results of calculations of the HKA test
Store the list of estimates of theta for each locus, as well as f (ratio of Ne of spp B relative to sppA), and T, the divergence time between species (on the scale of the Ne of spp A).
typedef std::tuple<double,double,double,double> Sequence::HKAresults::chisq_tuple |
|
strong |
|
explicit |
_thetas | a vector containing the estimats of theta for each locus |
_chisquareds | a vector of the per-locus contributions to the test statistic |
_fhat | an estimate of the relative Ne of species b to species a |
_That | an estimate of the divergence time, on the scale of Ne in species a. |
_xsq | the Chi-squared statistic for the HKA test |
_xsqA | The HKA chi-squared statistic, only considering species A |
_xsqB | The HKA chi-squared statistic, only considering species B |
std::vector< chisq_tuple > Sequence::HKAresults::chisquareds |
The order of elements is the same as the loci that were placed in the std::vector< HKAdata > passed to calcHKA. Members of the tuple are accessed via the std::get<int> template function. For example, if you have an HKAresults object called hkares, and you want to retrieve the deviation in divergence of the first locus:
The template argument to std::get can be taken from chisq_tuple_elements so that it's readable.
const double Sequence::HKAresults::fhat |
const double Sequence::HKAresults::That |
const std::vector< double > Sequence::HKAresults::thetas |
const double Sequence::HKAresults::xsq |
const double Sequence::HKAresults::xsqA |
const double Sequence::HKAresults::xsqB |