Compare commits
6 commits
Author | SHA1 | Date | |
---|---|---|---|
a5c85911d0 | |||
561f88be60 | |||
6c05ec1de7 | |||
3c9cdd67e0 | |||
ec4fe0222b | |||
2609a4fbd4 |
2 changed files with 18 additions and 10 deletions
|
@ -1,7 +1,5 @@
|
|||
# guoyi.run
|
||||
|
||||
[](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
|
||||
|
||||
|
|
22
guoyi.run
22
guoyi.run
|
@ -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 |;
|
||||
|
|
Loading…
Add table
Reference in a new issue