Natalia,
creio que a classe geodata não é compatível com indentify que é compatível com data frame(list)...
tente:
identify(geo[,2],geo[,3])
Att
Éder

Em 22 de agosto de 2012 18:59, Natalia Martins <nsmbarreto@gmail.com> escreveu:
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@polygons[[1]]@Polygons[[1]]@coords)
polygon(meso@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
 



_______________________________________________
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.