[R-br] Dúvida - elipse de confiança

Mauro Sznelwar sznelwar em uol.com.br
Sexta Janeiro 24 01:47:44 BRST 2014


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 = "???")
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140124/f2fee60c/attachment.html>


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