Compare commits

...

6 commits
0.1.3 ... main

2 changed files with 18 additions and 10 deletions

View file

@ -1,7 +1,5 @@
# guoyi.run
[![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.8431529.svg)](https://doi.org/10.5281/zenodo.8431529)
TNT Script for Maximum Parsimony Analysis (Phylogeny)
## Usage
@ -58,7 +56,9 @@ tnt run guoyi.run filename datatype weight 0/K cons resample prefix,
- resample should be sum of what you want
- relative bremer support (rbrs)=0.1, bremer support (brs)=0.2, jackknifing (jak)=1, bootstrap (boot)=2, symmetric resampling (sym)=4 i.e. rbrs+jak+boot+sym=7.1 (default)
- relative bremer support (rbrs)=0.1, bremer support (brs)=0.2, jackknifing (jak)=1, bootstrap (boot)=2, symmetric resampling (sym)=4 i.e. jak+boot+sym=7 (default)
- Notice: Bremer related support would change the tree topology, so it should be used with caution.
- prefix can be empty or a string

View file

@ -111,7 +111,7 @@ end
/*handle resmaple type*/
if (argnumber<6)
set dojak 1; set doboot 1; set dosym 1; set dorbrs 1;
set dojak 1; set doboot 1; set dosym 1; set dorbrs 0;
else
set rsmp %6;
loop 1 5
@ -356,7 +356,6 @@ set contree ntrees;
/*Store consensus tree to tree vault*/
hold /+0;
tv>/;
tchoose 0.'npars';
/*Get rbr/br/jak/boot/sym support and get consensus tree*/
ttags=;
@ -377,20 +376,21 @@ if ('dobremer')
/*calculate relative bremer support*/
if ('dorbrs')
bs ]!! 0;
bs ]!! 0.'npars';
end
/*calculate bremer support*/
if ('dobrs')
macfloat 1; /*set the br value float*/
bsupport !! 0;
bsupport !! 0.'npars';
end
end
/*Choose final tree*/
tchoose 0.'npars';
tv<;
tchoose/;
tchoose /;
/*adjust to zero following Pablo Goloboff*/
sub 0;
@ -426,6 +426,10 @@ taxname=;
/*Export consensus tree*/
tchoose /;
ttags=;
tp/;
ttags & %7.original.svg thickness 7 italics fontsize 15;
ttags-;
export - %7.original.tre;
taxname-;
export - %7.original_no.tre;
@ -545,7 +549,11 @@ if ('dohlf')
quote
| half strict consensus tree with |;
end end end
if ('dobremer')
quote
| WARNING: bremer-related will |
| WARNING: change the polytomy |;
end
if ('dorbrs')
quote
| relative bremer support |;