Extraindo os raios e área de uma elipse ellipsoidhull() do pacote cluster
Caros Membros, Com fins ilustrativos e exploratório, tenho pontos correspondentes a dois tratamentos A e B espacialmente distribuídos no campo, estou delimitando o limite (zona de influência) de cada um como uma elipse, através da função ellipsoidhull() do pacote cluster conforme CRM abaixo. Gostaria de saber se alguém familiarizado com a função poderia me dizer como faço para extrair os raios e as áreas das elipses EpointsA e EpointsB? #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 ellipsoidhull require(cluster) PointsA<-unname(cbind(dados$x[1:6],dados$y[1:6])) EpointsA<-ellipsoidhull(PointsA) lines(predict(EpointsA), col="blue") PointsB<-unname(cbind(dados$x[7:12],dados$y[7:12])) EpointsB<-ellipsoidhull(PointsB) lines(predict(EpointsB), col="orange") # Abraços, -- ====================================================================== 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 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 - ResearcherID: A-5790-2016 Researchgate: www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635 Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/ ======================================================================
Os valores retornados por EpointsA$d2 e EpointsB$d2 " average squared radius. " cfme. a documentação não lhe são suficientes? 2018-05-23 11:14 GMT-03:00 ASANTOS via R-br <r-br@listas.c3sl.ufpr.br>:
Caros Membros,
Com fins ilustrativos e exploratório, tenho pontos correspondentes a dois tratamentos A e B espacialmente distribuídos no campo, estou delimitando o limite (zona de influência) de cada um como uma elipse, através da função ellipsoidhull() do pacote cluster conforme CRM abaixo. Gostaria de saber se alguém familiarizado com a função poderia me dizer como faço para extrair os raios e as áreas das elipses EpointsA e EpointsB?
#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 ellipsoidhull require(cluster) PointsA<-unname(cbind(dados$x[1:6],dados$y[1:6])) EpointsA<-ellipsoidhull(PointsA) lines(predict(EpointsA), col="blue") PointsB<-unname(cbind(dados$x[7:12],dados$y[7:12])) EpointsB<-ellipsoidhull(PointsB) lines(predict(EpointsB), col="orange") #
Abraços,
-- ====================================================================== 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 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 - ResearcherID: A-5790-2016 Researchgate: www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635 Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/ ======================================================================
_______________________________________________ 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 fornea cdigo mnimo reproduzvel.
Muito obrigado César, vi a documentação, mas mesmo assim submeti para ver com quem trabalha com a ferramenta, se essa abordagem baseada na área e raios, está correta, uma vez que a ferramenta não é empregada para essa finalidade que propus. Em qui, 24 de mai de 2018 08:03, Cesar Rabak <cesar.rabak@gmail.com> escreveu:
Os valores retornados por EpointsA$d2 e EpointsB$d2 " average squared radius. " cfme. a documentação não lhe são suficientes?
2018-05-23 11:14 GMT-03:00 ASANTOS via R-br <r-br@listas.c3sl.ufpr.br>:
Caros Membros,
Com fins ilustrativos e exploratório, tenho pontos correspondentes a dois tratamentos A e B espacialmente distribuídos no campo, estou delimitando o limite (zona de influência) de cada um como uma elipse, através da função ellipsoidhull() do pacote cluster conforme CRM abaixo. Gostaria de saber se alguém familiarizado com a função poderia me dizer como faço para extrair os raios e as áreas das elipses EpointsA e EpointsB?
#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 ellipsoidhull require(cluster) PointsA<-unname(cbind(dados$x[1:6],dados$y[1:6])) EpointsA<-ellipsoidhull(PointsA) lines(predict(EpointsA), col="blue") PointsB<-unname(cbind(dados$x[7:12],dados$y[7:12])) EpointsB<-ellipsoidhull(PointsB) lines(predict(EpointsB), col="orange") #
Abraços,
-- ====================================================================== 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 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 - ResearcherID: A-5790-2016 Researchgate: www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635 Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/ ======================================================================
_______________________________________________ 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 fornea cdigo mnimo reproduzvel.
OK Alexandre! De fato, como eu não tenho conhecimento do domínio do problema que você estuda, me soa "não intuitivo" que se use uma ferramenta para medir uma área de influência sem levar em conta alguma 'resposta' dos tipos de tratamento... Boa sorte na busca de um pacote ou funções no R para resolver seu problema! []s 2018-05-24 12:06 GMT-03:00 Alexandre Dos Santos < alexandre.santos@cas.ifmt.edu.br>:
Muito obrigado César, vi a documentação, mas mesmo assim submeti para ver com quem trabalha com a ferramenta, se essa abordagem baseada na área e raios, está correta, uma vez que a ferramenta não é empregada para essa finalidade que propus.
Em qui, 24 de mai de 2018 08:03, Cesar Rabak <cesar.rabak@gmail.com> escreveu:
Os valores retornados por EpointsA$d2 e EpointsB$d2 " average squared radius. " cfme. a documentação não lhe são suficientes?
2018-05-23 11:14 GMT-03:00 ASANTOS via R-br <r-br@listas.c3sl.ufpr.br>:
Caros Membros,
Com fins ilustrativos e exploratório, tenho pontos correspondentes a dois tratamentos A e B espacialmente distribuídos no campo, estou delimitando o limite (zona de influência) de cada um como uma elipse, através da função ellipsoidhull() do pacote cluster conforme CRM abaixo. Gostaria de saber se alguém familiarizado com a função poderia me dizer como faço para extrair os raios e as áreas das elipses EpointsA e EpointsB?
#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 ellipsoidhull require(cluster) PointsA<-unname(cbind(dados$x[1:6],dados$y[1:6])) EpointsA<-ellipsoidhull(PointsA) lines(predict(EpointsA), col="blue") PointsB<-unname(cbind(dados$x[7:12],dados$y[7:12])) EpointsB<-ellipsoidhull(PointsB) lines(predict(EpointsB), col="orange") #
Abraços,
-- ====================================================================== 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 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 - ResearcherID: A-5790-2016 Researchgate: www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635 Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/ ======================================================================
_______________________________________________ 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 fornea cdigo mnimo reproduzvel.
participantes (3)
-
Alexandre Dos Santos -
ASANTOS
-
Cesar Rabak