Mercurial > hg > plosone_underreview

diff scripts_R/PlotOutliersCountry.R @ 80:103f7411c3ad branch-tests

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
changed interactive plot
author Maria Panteli <m.x.panteli@gmail.com>
date Tue, 26 Sep 2017 14:38:53 +0100
parents bde45ce0eeab
children d3e05cd49feb
line wrap: on
line diff
--- a/scripts_R/PlotOutliersCountry.R	Tue Sep 26 12:40:07 2017 +0100
+++ b/scripts_R/PlotOutliersCountry.R	Tue Sep 26 14:38:53 2017 +0100
@@ -16,7 +16,7 @@
 distMahal = as.dist(apply(data, 1, function(i) mahalanobis(data, i, cov = cov(data),tol=1e-18)))
 hc=hclust(distMahal, method="average")
 mypal = c("#000000", "#9B0000", "#9B0000", "#9B0000", "#9B0000")
-clus5 = cutree(hc, 5)
+clus5 = cutree(hc, 4)
 pdf('../data/results/hierarchical_cluster.pdf')
 par(mar=c(1,1,1,1))
 plot(as.phylo(hc),type="fan",tip.color=mypal[clus5], cex=.5, label.offset=.5)