[R-br] combinar Cluster e Ordination (PCA)

Jackson Rodrigues jacksonmrodrigues em gmail.com
Segunda Junho 30 10:13:33 BRT 2014


Olá todos,

Gostaria de uma ajuda para plotar Principal Components Analysis (PCA) e
clusters juntos.
Possuo dados ambientais de 25 locais espalhados em 5 diferentes
ecossistemas que quando agrupados em 5 clusters, locais de ecossitemas
diferentes se juntam no mesmo grupo. Logo, uma bagunça.
Assim, quero plotar os 2 juntos (PCA e Clusters) de maneira que, locais
pertencentes ao mesmo ecossistema possuam as mesmas CORES, e locais
agrupados no mesmo cluster, possuam mesma FORMA. Por exemplo:
1 local da Floresta de Araucária e 1 local da Floresta Decidual (logo com
cores diferentes azul e verde respectivamente) estão no mesmo cluster
(logo, representados pela mesma figura, triângulo).

Como fazer essa combinação?
Abaixo segue os comandos que estou usando para fazer minhas combinação, no
entanto. Nesse código, consigo juntar PCA e Cluster. no entanto sem
diferenciar as cores de cada local de acordo com eu ecossistema.

#PCA
mydata.pca<-rda(mydata)
#Cluster
mydata.w <- hclust(dist(Sqchord.mydata), "ward")
plot(mydata.w, hang=-1)
rect.hclust(mydata.w, 5)

# Cut the dendrogram to yield 5 groups
gr <- cutree(mydata.w, k=5)
grl <- levels(factor(gr))

# Get the site scores, scaling 1
sit.sc1 <- scores(mydata.pca, display="w", scaling=3)

# Plot the sites with cluster symbols and colours (scaling 3)
p <- plot(mydata.pca,  display="wa", scaling=3, type="n",
          main="Mydata PCA and clusters")
abline(v=0, lty="dotted")
abline(h=0, lty="dotted")
for (i in 1:length(grl))
{  points(sit.sc1[gr==i,], pch=(14+i), cex=2.5, col=i+1)}
text(sit.sc1, row.names(Pre_euro_veg.1.all), cex=0.7, pos=3)

# Add legend interactively
legend(locator(1), paste("Cluster", c(1:length(grl))),
pch=14+c(1:length(grl)),
       col=1+c(1:length(grl)), pt.cex=2)

Desde já, agradeço por qualquer ajuda.

Jackson M. Rodrigues
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140630/39df7a78/attachment.html>


Mais detalhes sobre a lista de discussão R-br