diff --git a/RF_old.R b/RF_old.R index 3ac0e75..9ef054c 100644 --- a/RF_old.R +++ b/RF_old.R @@ -8,7 +8,6 @@ library("readtext") #plot(unrooted_tree) #rooted_tree=reroot(unrooted_tree,n) #plot(rooted_tree) -#is.rooted(rooted_tree) #rooted_tree=drop.tip(rooted_tree,n) #plot(rooted_tree) #tree=rooted_tree @@ -138,7 +137,7 @@ qmT=function(R,n,m){ } return(qmt) } - +#this fonction computes the RF distribution polynomial=function(tree,n){ Coef=numeric() R=RF_old(tree,n) @@ -151,4 +150,3 @@ polynomial=function(tree,n){ -