diff nb_clust_G.R @ 1:e94a25eed489 draft

planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Ecoregionalization_workflow commit 459ba1277acd7d8d4a02f90dbd7ff444bf8eac92

--- a/nb_clust_G.R	Wed Oct 18 09:59:06 2023 +0000
+++ b/nb_clust_G.R	Wed Jan 24 15:53:32 2024 +0000
@@ -21,7 +21,7 @@
 
 #load data 
 
-env.data <- read.table(enviro, header = TRUE, dec = ".", na.strings = "-9999.00") 
+env.data <- read.table(enviro, sep="\t", header = TRUE, dec = ".", na.strings = "-9999") 
 
 ##List of modelled taxa used for clustering
 tv <- read.table(taxa_list, dec=".", sep=" ", header=F, na.strings = "NA") 
@@ -33,7 +33,7 @@
 data.bio = NULL
 
 for (i in 1:length(data_split[[1]])) {
-data.bio1 <- read.table(data_split[[1]][i], dec=".", sep=" ", header=T, na.strings = "NA")
+data.bio1 <- read.table(data_split[[1]][i], dec=".", sep="\t", header=T, na.strings = "NA")
 data.bio <- rbind(data.bio,data.bio1)
 remove(data.bio1)
 }
@@ -69,5 +69,5 @@
 
 # Plot SIH Index Chart by Number of Clusters
 png("Indices_SIH.png")
-plot(2:max_k, sih_values[2:max_k], type = "b", xlab = "Nombre de clusters", ylab = "Indice SIH")
+plot(2:max_k, sih_values[2:max_k], type = "b", xlab = "Number of clusters", ylab = "SIH index")
 dev.off()