Dúvida - elipse de confiança

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 = "???")
participantes (1)
-
Mauro Sznelwar