diff --git a/01b_MesobromionCluster.R b/01b_MesobromionCluster.R index 3b60cdb79b7519e840a791e6a8e1dee84d1c1e22..7fef7c6d37173629f59b39eeb6520c540fee4bd1 100644 --- a/01b_MesobromionCluster.R +++ b/01b_MesobromionCluster.R @@ -216,7 +216,7 @@ Mesobromion <- function(species.path, traits.path, output, myfunction="get.corXY if(b==1) {cor.obs <- rbind(cor.obs, tmp)} cor.bootstr[i,b] <- tmp$Coef cor.perm[i,b] <- myfunction(comm=speciesb, trait=traitsb.perm, trait.sel=tt, stat="RV")$Coef - print(paste("trait", tt, "perm", b)) + print(paste("trait", paste0(tt, collapse="_"), "perm", b)) } print(i) } diff --git a/submit_01b.sh b/submit_01b.sh index 55cd08702aad53e5a66ce292f5ba24a18bf087a2..b165c0347cbaf6151e23bf494392dd110df9b668 100644 --- a/submit_01b.sh +++ b/submit_01b.sh @@ -22,10 +22,10 @@ output=/data/splot/HIDDEN/output/HIDDEN Rscript \ cli_01b.r \ - --max.inter.t 2 \ + --max.inter.t 3 \ --chunk.i $SGE_TASK_ID \ --chunkn 399 \ - --nperm 99 \ + --nperm 199 \ /data/splot/HIDDEN/species.out.10perc.txt \ /data/splot/HIDDEN/traits.out.10perc.txt \ "$output" \