fix: themin=themax

README.md

@@ -103,3 +103,5 @@ Cite this script is mandatory, list me (Guoyi Zhang) in your Acknowledgements is
 # Further information
 
 To make tnt exported trees can be readable by WinClada (Nixon, 2021) and FigTree, please check tnt2winclada and tnt2figtree in this repository.
+
+runwincladtree.run can handle [wincladtree script](https://www.lillo.org.ar/phylogeny/tnt/scripts/wincladtree.run) with input file which must be $(PREFIX).winclada.ss or winclada.tree.

guoyi.run

@@ -463,13 +463,21 @@ set themin minsteps;
 set themax maxsteps;
 set TL length[0];
 set CI 'themin'/'TL'; /*CI=1 means no homoplasy*/
-set RI ('themax'-'TL')/('themax'-'themin'); /*RI=1 character fits perfetcly*/
+if ('themax' != 'themin')
+	set RI ('themax'-'TL')/('themax'-'themin'); /*RI=1 character fits perfetcly*/
+else
+	set RI 999;
+end
 
 /*Report CI/RI/TL */
 log %7.report.log;
 macfloat 3;
 quote Consistency Index (CI) is 'CI';
-quote Retention Index (RI) is 'RI';
+if ('RI' == 999)
+	quote Retention Index (RI) is NA;
+else
+	quote Retention Index (RI) is 'RI';
+end
 quote Tree Length (TL) is 'TL';
 log/;
 

runwincladtree.run

@@ -0,0 +1,7 @@
+macro=;
+taxname+5000;
+p %1;
+ttags=;
+apo>0;
+wincladtree %1.wincladtree.svg noalign shade rotate squaresize 15 taxangle 0;
+zzz;