report mixture model component in .iqtree file if necessary (requeste…

…d by Joran Martijn)
Cibiv · Dec 8, 2017 · 6be8320 · 6be8320
1 parent 32a5317
commit 6be8320
Showing 1 changed file with 8 additions and 0 deletions.
diff --git a/main/phyloanalysis.cpp b/main/phyloanalysis.cpp
@@ -400,6 +400,14 @@ void reportModel(ofstream &out, PhyloTree &tree) {
                 reportModel(out, tree.aln, m);
             }
 		}
+        if (tree.aln->seq_type != SEQ_POMO)
+		for (i = 0; i < nmix; i++) {
+            ModelMarkov *m = (ModelMarkov*)mmodel->getMixtureClass(i);
+            if (m->getFreqType() == FREQ_EQUAL || m->getFreqType() == FREQ_USER_DEFINED)
+                continue;
+            out << endl << "Model for mixture component "  << i+1 << ": " << (m)->name << endl;
+            reportModel(out, tree.aln, m);
+        }
 		out << endl;
 	} else {
 		out << "Model of substitution: " << tree.getModelName() << endl << endl;