Elipse circundados pontos e cálculo do centro de massa
Boa tarde Listeiros, Tenho dois conjuntos de pontos A e B (dados$Class) e suas coordenadas (dados$x e dados$y), gostaria de saber se é possível desenhar uma elipse que circunde todos os pontos A e outra elipse que circunde todos os pontos B, posteriormente gostaria de calcular o centro de massa para cada conjunto de pontos A e B, isso é possível? Segue 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") # 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@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================
há algumas funções dentro do pacote 'car' para gerar elipses. 1 ) Ellipses 2 ) Data Ellipses 3) and Confidence Ellipses veja se alguma parte disso pode te ajudar cleber Em 20/09/2015 17:27, ASANTOS escreveu:
Boa tarde Listeiros,
Tenho dois conjuntos de pontos A e B (dados$Class) e suas coordenadas (dados$x e dados$y), gostaria de saber se é possível desenhar uma elipse que circunde todos os pontos A e outra elipse que circunde todos os pontos B, posteriormente gostaria de calcular o centro de massa para cada conjunto de pontos A e B, isso é possível? Segue 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") #
Obrigado,
--- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus
Para traçar elipse que contenha os pontos você precisa definir melhor. Cleber sugeriu as funções do pacote car 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.
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@yahoo.com.br alexandre.santos@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@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.
participantes (3)
-
ASANTOS -
Cleber N.Borges -
Walmes Zeviani