16 int main(
int argc,
char **argv)
18 const char *infilename = argv[1];
20 vector<Sequence::Fasta> data;
22 Sequence::Alignment::GetData(data,infilename);
24 if ( Sequence::Alignment::IsAlignment(data) &&
25 Sequence::Alignment::validForPolyAnalysis(data.begin(),data.end()) )
27 const unsigned alignmentLength = data[0].length();
33 cout <<
"Tajima's D for the whole dataset is: " 34 << analyzeRegion.TajimasD()
43 for(
unsigned i = 0 ; i < windows.size() ; ++i)
47 cout <<
"Tajima's D for window " 50 << 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.
Declaration of Sequence::Fasta streams.
A container class for "sliding windows" along a polymorphism table.