21 int main(
int argc,
char **argv)
23 const char * infilename = argv[1];
24 std::vector<Sequence::Fasta> data;
25 Sequence::Alignment::GetData(data,infilename);
27 if ( Sequence::Alignment::IsAlignment(data) &&
28 Sequence::Alignment::validForPolyAnalysis(data.begin(),data.end()) )
31 if (! SNPtable.empty())
34 std::cout <<
"Tajima's D for the region is: "<< analyzeRegion.TajimasD() << std::endl;
37 SNPtable.send(),1u,1u);
39 while(itr < windows.
cend())
43 std::cout <<
"D for window " 46 << analyzeWindow.TajimasD()
Polymorphism tables for sequence data.
Molecular population genetic analysis.
declaration of Sequence::PolySNP, a class to analyze SNP data
Declaration of namespace Sequence::Alignment.
Sequence::PolySites, generates polymorphism tables from data.
const_iterator cend() const
Declaration of Sequence::Fasta streams.
A container class for "sliding windows" along a polymorphism table.
std::vector< range >::const_iterator const_iterator
const_iterator cbegin() const