|
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 |