[R-br] como identificar pontos no mapa

Natalia Martins nsmbarreto em gmail.com
Quarta Agosto 22 18:59:44 BRT 2012


Boa noite pessoal,
estou tentando identificar uns pontos em um mapa.
Os comandos que utilizei são os seguintes:

geo = data.frame(dados$idestacao,dados$latitude,dados$longitude)
head(geo)
  dados.idestacao dados.latitude dados.longitude
1               1       -24.3209        -54.0111
2               1       -24.3209        -54.0111
3               1       -24.3209        -54.0111
4               1       -24.3209        -54.0111
5               1       -24.3209        -54.0111
6               1       -24.3209        -54.0111

dadosgeo <- as.geodata(geo,coords.col=c(3,2),data.col=1)
points(dadosgeo)
meso = readShapePoly("E:/natalia/41ME2500G.shp")
estado = readShapePoly("E:/natalia/41UF2500G.shp")
polygon(estado em polygons[[1]]@Polygons[[1]]@coords)
polygon(meso em polygons[[4]]@Polygons[[1]]@coords)
identify(dadosgeo)

mas quando peço a identificação aparece o seguinte erro:
'x' is a list, but does not have components 'x' and 'y'

Não entendi onde estou errando, no entanto quando pedi pra ver meus
dadosgeo vi que eles tem apenas  as  idestacao, as quais deveriam ser
identificada no mapa.

Muito obrigada.

Natália da Silva Martins
Bacharel em Estatística - Universidade Estadual de Maringá/ UEM
Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP
Contato: (19) 8306-4743
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120822/87f1b4f0/attachment.html>


Mais detalhes sobre a lista de discussão R-br