[R-br] Plotar objeto polylist pacote spdep

[Helder Gramacho]

Olá Éder,

Boa tarde! muito obrigado pela ajuda.

*Hélder Gramacho *
Recife-PE /
*agrohelder em gmail.com <agrohelder em hotmail.com>*



Em 20 de junho de 2014 10:33, Éder Comunello <comunello.eder em 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 em gmail.com>omunello.eder em gmail.com>
> Dourados, MS - [22 16.5'S, 54 49'W]
>
>
> Em 19 de junho de 2014 23:10, Helder Gramacho <agrohelder em 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_tutoriais&cache=cache&media=ser301:aula_r_regressao_anselin.txt
>>
>> 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 em gmail.com <agrohelder em hotmail.com> *
>>
>>
>> _______________________________________________
>> 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.
>>
>
>
> _______________________________________________
> 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/20140620/4790557c/attachment.html>