Plotar objeto polylist pacote spdep
Boa noite pessoal, Estou praticando o tutorial sobre o pacote spdep e modelos de regressão espacial disponível no link: http://wiki.dpi.inpe.br/lib/exe/fetch.php?id=materiais_para_leitura_e_tutori... O tutorial utiliza o dataset columbus, o qual possui um objeto polys que é da classe polylist e que possui 49 poligonos. Numa determinada etapa do tutorial é utilizado o comando plot para exibir os poligonos: ### <code r> library(spdep) data(columbus) plot(polys) plot(polys,bbs,border="grey") ### </code> entretanto quando executo aparecem os seguintes erros:
plot(polys) Error in xy.coords(x, y, xlabel, ylabel, log) : 'x' is a list, but does not have components 'x' and 'y'
plot(polys,bbs,border="grey") # não funciona Error in xy.coords(x, y, xlabel, ylabel, log) : 'x' and 'y' lengths differ
Alguém sabe o que pode estar acontecendo, e como corrigir? Desde já agradeço, *Hélder Gramacho * Recife-PE / *agrohelder@gmail.com <agrohelder@hotmail.com>*
Helder, bom dia! Eu não sei dizer com certeza, mas tenho a impressão que a classe 'polylist' se tornou obsoleta. Me parece que as funções de certos pacotes que tratavam da 'polylist', não o fazem mais (ex: {maptools}). Alternativamente, é possível converter para objetos do pacote {sp} ou usar as funções básicas (lines, polygons)... ### <code r> library(spdep) data(columbus) plot(namevec=polys) plot(polys,bbs,border="grey") ### </code> ### conversão básica para SpatialPolygons Pol1 <- sapply(polys, function(x) Polygon(as.data.frame(x))) SPol1 <- SpatialPolygons(list(Polygons(Pol1, "All"))) plot(SPol1, col=terrain.colors(8)) ### conversão com identificadores attributes(polys) Pol2 <- sapply(polys, function(x) Polygon(as.data.frame(x))); names(Pol2) <- attr(polys, 'region.id') Pols2 <- lapply(1:length(Pol2), function(x) Polygons(Pol2[x], attr(polys, ' region.id')[x])) SPol2 <- SpatialPolygons(Pols2) plot(SPol2, col=terrain.colors(8)) ### </code> Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W] Em 19 de junho de 2014 23:10, Helder Gramacho <agrohelder@gmail.com> escreveu:
Boa noite pessoal,
Estou praticando o tutorial sobre o pacote spdep e modelos de regressão espacial disponível no link:
http://wiki.dpi.inpe.br/lib/exe/fetch.php?id=materiais_para_leitura_e_tutori...
O tutorial utiliza o dataset columbus, o qual possui um objeto polys que é da classe polylist e que possui 49 poligonos.
Numa determinada etapa do tutorial é utilizado o comando plot para exibir os poligonos: ### <code r> library(spdep) data(columbus) plot(polys) plot(polys,bbs,border="grey") ### </code>
entretanto quando executo aparecem os seguintes erros:
plot(polys) Error in xy.coords(x, y, xlabel, ylabel, log) : 'x' is a list, but does not have components 'x' and 'y'
plot(polys,bbs,border="grey") # não funciona Error in xy.coords(x, y, xlabel, ylabel, log) : 'x' and 'y' lengths differ
Alguém sabe o que pode estar acontecendo, e como corrigir?
Desde já agradeço, *Hélder Gramacho * Recife-PE / *agrohelder@gmail.com <agrohelder@hotmail.com> *
_______________________________________________ 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.
Olá Éder, Boa tarde! muito obrigado pela ajuda. *Hélder Gramacho * Recife-PE / *agrohelder@gmail.com <agrohelder@hotmail.com>* Em 20 de junho de 2014 10:33, Éder Comunello <comunello.eder@gmail.com> escreveu:
Helder, bom dia!
Eu não sei dizer com certeza, mas tenho a impressão que a classe 'polylist' se tornou obsoleta. Me parece que as funções de certos pacotes que tratavam da 'polylist', não o fazem mais (ex: {maptools}).
Alternativamente, é possível converter para objetos do pacote {sp} ou usar as funções básicas (lines, polygons)...
### <code r> library(spdep) data(columbus) plot(namevec=polys) plot(polys,bbs,border="grey") ### </code>
### conversão básica para SpatialPolygons Pol1 <- sapply(polys, function(x) Polygon(as.data.frame(x))) SPol1 <- SpatialPolygons(list(Polygons(Pol1, "All"))) plot(SPol1, col=terrain.colors(8))
### conversão com identificadores attributes(polys) Pol2 <- sapply(polys, function(x) Polygon(as.data.frame(x))); names(Pol2) <- attr(polys, 'region.id') Pols2 <- lapply(1:length(Pol2), function(x) Polygons(Pol2[x], attr(polys, ' region.id')[x])) SPol2 <- SpatialPolygons(Pols2) plot(SPol2, col=terrain.colors(8))
### </code>
Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
Em 19 de junho de 2014 23:10, Helder Gramacho <agrohelder@gmail.com> escreveu:
Boa noite pessoal,
Estou praticando o tutorial sobre o pacote spdep e modelos de regressão espacial disponível no link:
http://wiki.dpi.inpe.br/lib/exe/fetch.php?id=materiais_para_leitura_e_tutori...
O tutorial utiliza o dataset columbus, o qual possui um objeto polys que é da classe polylist e que possui 49 poligonos.
Numa determinada etapa do tutorial é utilizado o comando plot para exibir os poligonos: ### <code r> library(spdep) data(columbus) plot(polys) plot(polys,bbs,border="grey") ### </code>
entretanto quando executo aparecem os seguintes erros:
plot(polys) Error in xy.coords(x, y, xlabel, ylabel, log) : 'x' is a list, but does not have components 'x' and 'y'
plot(polys,bbs,border="grey") # não funciona Error in xy.coords(x, y, xlabel, ylabel, log) : 'x' and 'y' lengths differ
Alguém sabe o que pode estar acontecendo, e como corrigir?
Desde já agradeço, *Hélder Gramacho * Recife-PE / *agrohelder@gmail.com <agrohelder@hotmail.com> *
_______________________________________________ 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.
_______________________________________________ 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 (2)
-
Helder Gramacho -
Éder Comunello