Preciso traçar elipse de confiança no fim desta sequencia,
queria saber como entro com os dados, que está no último comando
abaixo:
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 = "???")