diff --git a/01_Mesobromion.R b/01_Mesobromion.R index d6fd595e08f3ac6c76470e59bea1b4efa4003fa2..238aae563bfddee8309384980fdbb0c02763fd6d 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 720a14e60b8bd6e760e6e1b2e5cf14c0a09f4e06..8b1a9fa0d6d2014e6ce2b526f21c4d476db90759 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="")) } }