Scripts/RF distances

From TNT

Jump to: navigation, search

RF Distances
Author: Pablo Goloboff.
Version: 1.0.
For: not provided.
Arguments: not provided.
Example: not provided.
Script
Scripts may not be displaying correctly here at present, use the script from the provide file in the archive.


Calculates Robinson-Foulds distance between trees macro=;

if ( argnumber < 2 ) errmsg Give trees to compare as argument; end 
set 0 collapse ; 
collapse none ;
tcomp * %1 %2 ;
tcomp * %2 %1 ;
set 1 ((tnodes[ntrees]+tnodes[(ntrees-1)])/(tnodes[%1]+tnodes[%2])); 
collapse '0' ;
set 0 ntrees - 1 ;
keep '0' ; 

quote Robinson-Foulds distance between trees %1-%2 is '/.3(1)' ;

proc/;
Personal tools