[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