[R-br] Elipse circundados pontos e cálculo do centro de massa

ASANTOS alexandresantosbr em yahoo.com.br
Domingo Setembro 20 21:06:09 BRT 2015


Obrigado pelas sugestões Walmes, Mauro e Cleber,

            Respondendo a questão do Walmes, não é uma questão de elipse 
com intervalo de confiança mesmo, mas sim mera questão geométrica de 
escolher a elipse que envolva um conjunto de dados com a menor área, 
consegui avançar com as sugestões de vocês, mas agora preciso de ajuda 
para sistematizar a criação dos convexhull e elipses, no meu CRM tenho 2 
conjunto de dados, mas caso eu tivesse 80 conjuntos eu precisaria cria 
uma função para automatizar a criação desses convex's e elipses, alguém 
poderia me dizer qual um bom caminho, segue novo CRM:

#Criação do data frame
dados<-NULL

#Coordenadas x

dados$x<-c(0.350,0.487,0.637,0.775,0.825,0.087,0.237,0.400,0.575,0.737,0.062,0.212)

#Coordenadas y
dados$y<-c(0.025,0.087,0.050,0.025,0.125,0.187,0.150,0.162,0.212,0.237,0.362,0.337)

#Classificação
dados$Class<-c("A","A","A","A","A","A","B","B","B","B","B","B")

#Visualização
plot(dados$x,dados$y)
text(dados$x,dados$y+0.01,labels=dados$Class,col="red")
#

##Usando convexhull
require(grDevices)
pointsA<-as.data.frame(cbind(dados$x[1:6],dados$y[1:6]))
hpts <- chull(pointsA)
hpts <- c(hpts, hpts[1])
lines(pointsA[hpts, ], col="red")

#Usando ellipsoidhull
require(cluster)
PointsA<-unname(cbind(dados$x[1:6],dados$y[1:6]))
EpointsA<-ellipsoidhull(PointsA)
lines(predict(EpointsA), col="blue")
#

Obrigado,

-- 
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT                      CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM)   (+55) 65 9686-6970 (VIVO)
e-mails:alexandresantosbr em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================




Em 20/09/2015 17:30, Walmes Zeviani escreveu:
> Para traçar elipse que contenha os pontos você precisa definir melhor. 
> Cleber sugeriu as funções do pacote caar r para traçar elipses de 
> confiança, seria isso? Mas elipses de confiança seriam apropriadas 
> para um caso de variável aleatória bivariada como uma normal bivariada 
> ou estimadores de mínimos quadrados em problemas de regressão (\beta_0 
> e \beta_1) por exemplo. Essas coordenadas não se encaixem nessas 
> situações, certo? Existem infinitas elipses que podem ser traçadas 
> para envolver um conjunto de pontos. Talvez você queira algo como a 
> elipse de menor área que contenha um conjunto de pontos. Procure por 
> "draw ellipse over given points". Você pode usar também o convex hull, 
> o poligono convexo que cerce todas as observações. Talvez a partir do 
> convex hull seja mais fácil posicionar a elipse.
>
> À disposição.
> Walmes.
>>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> 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.

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150920/eb9b2228/attachment.html>


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