Mercurial > repos > ecology > ecoregion_cluster_estimate

comparison cluster_ceamarc.R @ 4:cf7aefee9977 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Ecoregionalization_workflow commit c41939f1cdc03331ec021d47495576a6b0c5fd14
author ecology
date Wed, 16 Oct 2024 11:45:55 +0000
parents e94a25eed489
children
comparison
equal deleted inserted replaced
3:001d7d101915 4:cf7aefee9977
39 plot(silhouette(test5)) 39 plot(silhouette(test5))
40 dev.off() 40 dev.off()
41 41
42 clus <- cbind(data.bio[1:nrow(test3), 1:2],test5$clustering) 42 clus <- cbind(data.bio[1:nrow(test3), 1:2],test5$clustering)
43 names(clus) <- c("lat", "long", "cluster") 43 names(clus) <- c("lat", "long", "cluster")
44 clus <- cbind(clus,test3,env.data[,3:19]) 44 clus <- cbind(clus,test3,env.data[,3:ncol(env.data)])
45 45
46 write.table(clus[1:3], file = "points_clus.tsv", quote = FALSE, sep="\t", row.names = FALSE) 46 write.table(clus[1:3], file = "points_clus.tsv", quote = FALSE, sep="\t", row.names = FALSE)
47 write.table(clus, file = "clus.tsv", quote = FALSE, sep="\t", row.names = FALSE) 47 write.table(clus, file = "clus.tsv", quote = FALSE, sep="\t", row.names = FALSE)
48 48
49 49