diff --git a/01_Mesobromion.R b/01_Mesobromion.R
index 796a0b2d787b3380a2a04cf4044cf688f0203a62..d6fd595e08f3ac6c76470e59bea1b4efa4003fa2 100644
--- a/01_Mesobromion.R
+++ b/01_Mesobromion.R
@@ -191,7 +191,7 @@ sign <- data.frame(Trait.comb=names(sign), sign=sign)
 corXY %>% 
   arrange(Test, desc(Coef))
 
-aa <- data.frame(Trait.comb=paste0("t", 1:81), 
+aa <- data.frame(Trait.comb=paste0("t", 1:ncol(traits)), 
                  trait.name=colnames(traits))#[-which(colnames(traits) %in% c("species", "species0"))])
 
 corXY.ci <- data.frame(Trait.comb=rownames(corXY.bootstr), 
@@ -264,7 +264,7 @@ traits.sign.alone <- mydata.SES %>%
   filter(ntraits==1) %>% 
   filter(sign.SES==T) %>% 
   pull(trait1)
-  
+
 mydata <- mydata.SES %>% 
   filter(ntraits==1)%>% 
   mutate(seq=1:nrow(.))
@@ -291,6 +291,12 @@ ggplot(data=mydata) +
 
 
 
+
+traits.sign <- traits %>% 
+  dplyr::select(traits.sign.alone)
+write_csv(traits.sign, "_data/Mesobromion/traits.out.10perc.sign.txt")
+
+
 break()
 
 
diff --git a/01b_MesobromionCluster.R b/01b_MesobromionCluster.R
index 60d753469de80b178d5474c349768370650f02f5..720a14e60b8bd6e760e6e1b2e5cf14c0a09f4e06 100644
--- a/01b_MesobromionCluster.R
+++ b/01b_MesobromionCluster.R
@@ -18,7 +18,7 @@ require(doParallel)
 source("99_HIDDEN_functions.R")
 
 
-Mesobromion <- function(species.path, traits.path, output, myfunction="get.corXY", max.inter.t, chunkn, chunk.i, nperm=199){
+Mesobromion <- function(species.path, traits.path, output, myfunction="get.corXY.bootstr", max.inter.t, chunkn, chunk.i, nperm=199){
 
   myfunction <- get(myfunction)
   ## calculate corr between species composition matrix and traits