本文整理汇总了C++中Sample::ConstantSiteCorrection方法的典型用法代码示例。如果您正苦于以下问题:C++ Sample::ConstantSiteCorrection方法的具体用法?C++ Sample::ConstantSiteCorrection怎么用?C++ Sample::ConstantSiteCorrection使用的例子?那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类Sample
的用法示例。
在下文中一共展示了Sample::ConstantSiteCorrection方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: main
//.........这里部分代码省略.........
cerr << "postscript in\t\t" << SampleName << ".con.tre.ps\n";
}
cerr << '\n';
if ((! mParam2) || (!mParam2->SaveAll)) {
exit(1);
}
}
// delete mParam2;
Sample* sample = new Sample(ChainName,burnin,every,until,Path);
if (save) {
sample->ToFile(SampleName);
}
MCParameters* mParam = sample->GetParameters();
// initialising
int Nsite = mParam->Nsite;
int Ntaxa = mParam->Ntaxa;
cout << '\n';
cout << "Nsite : " << Nsite << '\n';
cout << "Ntaxa : " << Ntaxa << '\n';
cout << sample->GetSize() << " points to read\n";
cout << '\n';
cout.flush();
if (mParam->NormalApprox) {
if (clustermodes) {
cerr << "no mode clustering under normal approx\n";
cerr << '\n';
exit(1);
}
if (modes) {
cerr << "no mode analysis under normal approx\n";
cerr << '\n';
exit(1);
}
if (rates) {
cerr << "no rate analysis under normal approx\n";
cerr << '\n';
exit(1);
}
if (sitestat) {
cerr << "no site-specific profile analysis under normal approx\n";
cerr << '\n';
exit(1);
}
}
if (clock) {
sample->Dating(ps);
sample->Reset();
}
else if (sitelogl) {
sample->ReadSiteLogLikelihood();
}
else if (intlnl) {
sample->ReadSummedLogLikelihood();
}
else if (entropyN) {
sample->ReadFiniteTimeEntropy(entropytimemax,entropyN);
}
else if (catgtrpoint) {
sample->ReadCATGTRPoint();
}
else if (popeff) {
sample->ReadPopEff(ps);
}
else if (branchfreqs) {
sample->ReadBranchFreqs();
}
else if (coaff) {
sample->Coaff(coaff,step);
}
else if (rr) {
sample->MeanRR();
}
else if (benchmark) {
sample->ReadBench();
}
else if (clustermodes) {
if (! mParam->SaveAll) {
cerr << "error : biochemical profiles were not saved. cannot cluster them\n";
exit(1);
}
sample->ClusterModes(minsize, mindist, 0);
}
else if (grepmode) {
sample->GrepModes();
}
else if (truelength) {
sample->TrueLength(nrep,sitebranch);
}
else if (constcorrect) {
sample->ConstantSiteCorrection(rho);
}
else {
sample->Read(rates, modes, sitestat, cutoff, ncat, ps);
}
}