<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META content="text/html; charset=iso-8859-1" http-equiv=Content-Type>
<META name=GENERATOR content="MSHTML 8.00.6001.23543">
<STYLE></STYLE>
</HEAD>
<BODY bgColor=#ffffff>
<DIV><FONT face=Arial>Preciso traçar elipse de confiança no fim desta sequencia,
queria saber como entro com os dados, que está no último comando
abaixo:</FONT></DIV>
<DIV><FONT color=#0000ff face=Arial><FONT color=#000000>Data set:</FONT><A
href="https://www.dropbox.com/s/405ggda0h6scnie/marajoara320amostras.txt"><FONT
face="Times New Roman">https://www.dropbox.com/s/405ggda0h6scnie/marajoara320amostras.txt</FONT></A></FONT></DIV>
<DIV><FONT size=2 face=Arial></FONT> </DIV>
<DIV><FONT size=2
face=Arial>marajoara<-read.table('marajoara320amostras.txt')</FONT></DIV>
<DIV><FONT size=2 face=Arial></FONT> </DIV>
<DIV><FONT size=2
face=Arial>=======NORMALIZAÇÃO========<BR>marajoara.norm<-sapply(data.frame(marajoara),scale)</FONT></DIV>
<DIV><FONT size=2 face=Arial></FONT> </DIV>
<DIV><FONT size=2 face=Arial>=======CLUSTER=========</FONT></DIV>
<DIV> </DIV>
<DIV><FONT size=2
face=Arial>win.graph()<BR>marajoara.norm.h<-hclust(dist(marajoara.norm)^2,method="ward")<BR>plot(marajoara.norm.h,main="marajoara320amostras")</FONT></DIV>
<DIV> </DIV>
<DIV><FONT size=2 face=Arial>=======PCA========</FONT></DIV>
<DIV> </DIV>
<DIV><FONT size=2
face=Arial>win.graph()<BR>initial<-tapply(marajoara.norm,list(rep(cutree(marajoara.norm.h,3),ncol(marajoara.norm)),col(marajoara.norm)),mean)<BR>marajoara.norm.km1<-kmeans(marajoara.norm,initial)<BR>marajoara.norm.pca<-princomp(marajoara.norm)<BR>marajoara.norm.px<-predict(marajoara.norm.pca)<BR>library(MASS)<BR>eqscplot(marajoara.norm.px[,1:2],type="n",xlab="componente
1",ylab="componente
2")<BR>text(marajoara.norm.px[,1:2],labels=as.character(marajoara.norm.km1$cluster),col=as.numeric(marajoara.norm.km1$cluster))<BR>title("marajoara320amostras")</FONT></DIV>
<DIV> </DIV>
<DIV><FONT size=2 face=Arial># identificar amostras no gráfico de
PCA</FONT></DIV>
<DIV> </DIV>
<DIV><FONT size=2 face=Arial>identify(marajoara.norm.px[,1:2])</FONT></DIV>
<DIV> </DIV>
<DIV><FONT size=2 face=Arial># Ellipse no gráfico de PCA
(********DUVIDA*********)</FONT></DIV>
<DIV> </DIV>
<DIV><FONT size=2 face=Arial>library(ellipse)<BR>plot(ellipse(marajoara.norm.px,
centre = ???, t = ???), type = "???")</FONT></DIV></BODY></HTML>