Compare commits

...

2 commits
0.1.5 ... main

Author SHA1 Message Date
a5c85911d0 fix: typo 2025-04-05 18:12:42 +11:00
561f88be60 polish: just add warn to polytomy 2025-04-05 18:04:33 +11:00
2 changed files with 13 additions and 7 deletions

View file

@ -56,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,18 +376,19 @@ if ('dobremer')
/*calculate relative bremer support*/
if ('dorbrs')
bs *]!! 0.'npars';
bs ]!! 0.'npars';
end
/*calculate bremer support*/
if ('dobrs')
macfloat 1; /*set the br value float*/
bsupport *!! 0.'npars';
bsupport !! 0.'npars';
end
end
/*Choose final tree*/
tchoose 0.'npars';
tv<;
tchoose /;
@ -549,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 |;