Preciso traçar elipse de confiança no fim desta sequencia, queria saber como
entro com os dados, que está no último comando abaixo:
Data set:https://www.dropbox.com/s/405ggda0h6scnie/marajoara320amostras.txt
marajoara<-read.table('marajoara320amostras.txt')
=======NORMALIZAÇÃO========
marajoara.norm<-sapply(data.frame(marajoara),scale)
=======CLUSTER=========
win.graph()
marajoara.norm.h<-hclust(dist(marajoara.norm)^2,method="ward")
plot(marajoara.norm.h,main="marajoara320amostras")
=======PCA========
win.graph()
initial<-tapply(marajoara.norm,list(rep(cutree(marajoara.norm.h,3),ncol(marajoara.norm)),col(marajoara.norm)),mean)
marajoara.norm.km1<-kmeans(marajoara.norm,initial)
marajoara.norm.pca<-princomp(marajoara.norm)
marajoara.norm.px<-predict(marajoara.norm.pca)
library(MASS)
eqscplot(marajoara.norm.px[,1:2],type="n",xlab="componente 1",ylab="componente
2")
text(marajoara.norm.px[,1:2],labels=as.character(marajoara.norm.km1$cluster),col=as.numeric(marajoara.norm.km1$cluster))
title("marajoara320amostras")
# identificar amostras no gráfico de PCA
identify(marajoara.norm.px[,1:2])
# Ellipse no gráfico de PCA (********DUVIDA*********)
library(ellipse)
plot(ellipse(marajoara.norm.px, centre = ???, t = ???), type = "???")_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
