diff --git a/01b_MesobromionCluster.R b/01b_MesobromionCluster.R
index 401220b55c2ea6129248d43af539e4c9fbfa6afd..03a024d282a33cbf8d661583953a1326b3007592 100644
--- a/01b_MesobromionCluster.R
+++ b/01b_MesobromionCluster.R
@@ -80,7 +80,8 @@ Mesobromion <- function(species.path, traits.path, output, myfunction="get.corXY
   
   traits <- traits %>% 
     column_to_rownames("species0") %>% 
-    rename_all(.funs=~gsub(pattern=".mean$", replacement="", x=.))  #%>% 
+    rename_all(.funs=~gsub(pattern=".mean$", replacement="", x=.))  %>% 
+    mutate_if(~is.character(.), .funs=~as.factor(.))
 #temporary    ### Use only a subset of traits
     #dplyr::select(LeafArea:Disp.unit.leng)