From f5b9a26100786f6327f8ba9d43112069e1739580 Mon Sep 17 00:00:00 2001
From: Francesco Sabatini <francesco.sabatini@idiv.de>
Date: Wed, 8 Apr 2020 19:08:30 +0200
Subject: [PATCH] fixed bug
---
01_Mesobromion.R | 4 ++--
01b_MesobromionCluster.R | 4 +---
2 files changed, 3 insertions(+), 5 deletions(-)
diff --git a/01_Mesobromion.R b/01_Mesobromion.R
index d6fd595..238aae5 100644
--- a/01_Mesobromion.R
+++ b/01_Mesobromion.R
@@ -293,9 +293,9 @@ ggplot(data=mydata) +
traits.sign <- traits %>%
+ rownames_to_column("species0") %>%
dplyr::select(traits.sign.alone)
-write_csv(traits.sign, "_data/Mesobromion/traits.out.10perc.sign.txt")
-
+write_delim(traits.sign, path="_data/Mesobromion/traits.out.10perc.sign.txt", delim="\t")
break()
diff --git a/01b_MesobromionCluster.R b/01b_MesobromionCluster.R
index 720a14e..8b1a9fa 100644
--- a/01b_MesobromionCluster.R
+++ b/01b_MesobromionCluster.R
@@ -59,15 +59,13 @@ Mesobromion <- function(species.path, traits.path, output, myfunction="get.corXY
# cor.perm <- matrix(NA, nrow=length(allcomb.t), ncol=nperm,
# dimnames = list(paste("t", names.t, sep=""), paste("p",1:nperm, sep="")))
- speciesb <- species
- traitsb <- traits
corXY.output <- NULL
for(i in 1:length(allcomb.t)) {
tt <- unlist(allcomb.t[i])
#bootstrap species matrix
# set.seed(1984)
corXY.output <- rbind(corXY.output,
- get.corXY.bootstrap(comm=speciesb, trait=traitsb.perm, trait.sel=tt, stat="RV", bootstrap=nperm))
+ get.corXY.bootstrap(comm=species, trait=traits, trait.sel=tt, stat="RV", bootstrap=nperm))
save(corXY.output, file = paste(output, "_", chunk.i, ".RData", sep=""))
}
}
--
GitLab