[R-br] em plot.xy(xy.coords(x, y), type = type, ...) : objeto 'membros' não encontrado

[Eder Comunello]

Kin, boa tarde!

O erro listado é por atribuir ao parâmetro pch o objeto 'membros' que não
existe ou não está acessível (veja seu arquivo de dados e declarações).
Acredito que a ideia seria diferenciar os pontos (tipos) de acordo com o
valor em membros...

Acredito que o exemplo que segue faça algo parecido com o que você deseja,
embora usando somente a biblioteca sp().

Inclui os pontos que você forneceu.

###########

require(sp)
br.df <- read.table("http://dl.dropbox.com/u/117618178/brlite.dat",
sep=";", head=T)
class(br.df); head(br.df); str(br.df)

br.lin <- Lines(list(Line(br.df)), ID = "Brasil")
br.sl  <- SpatialLines(list(br.lin))
proj4string(br.sl) <- CRS('+proj=longlat +datum=WGS84')
summary(br.sl); class(br.sl)

est.x  <- c(-53.35,-55.60,-51.80,-51.05,-50.13,-52.41,-52.10,-52.51)
est.y  <- c(-33.51,-30.83,-29.45,-29.68,-29.05,-31.78,-32.03,-30.53)
est.id <- c(1:8)
est.df <- data.frame(est.x,est.y,est.id)

est.sp <- est.df
coordinates(est.sp) <- c('est.x','est.y')
proj4string(br.sl) <- CRS('+proj=longlat +datum=WGS84')
est.sp

#png(filename="c:/temp/mapa.png",width=600,height=800)
plot(br.sl, axes=T)
plot(est.sp, pch=21, bg=3, add=T)
#dev.off()

###########

-- 

================================================
Éder Comunello

Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ)
Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W]

Researcher at Embrapa Western Region Agriculture
Dourados, MS, Brazil [22 16.5'S, 54 49.0'W]
================================================
UTC-03:00



> Date: Tue, 30 Oct 2012 21:09:26 -0200
> From: kin Huts <kim-met em hotmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=kim-met@hotmail.com>
> >
> To: <r-br em listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>
> >
> Subject: [R-br] em plot.xy(xy.coords(x, y), type = type, ...) :
>         objeto 'membros' não encontrado
> Message-ID: <SNT141-W57AFF1FD679C8C33FDDE18F1620 em phx.gbl>
> Content-Type: text/plain; charset="iso-8859-1"
>
>
>
> Olá,
> Gostaria de uma ajuda para entender o erro encontrado e como corrigir esse
> script,Estou tentando gerar um mapa para o RS com a localização  das
> estações que estou utilizando.
>
>
> require(MASS)require(maps)
>
> # GERAÇÃO DO MAPA COM AS ESTAÇÕES
> # carregando divisão política brasileiramapaBrasil <-
> read.table("BASE_BR.txt",sep="")
> # carregando informações de localização das estações
> meteorológicaslatlonRS <- read.csv2("coordenadas.txt",h=F,sep="")
> # gerando mapa em PNGpng(filename="mapa_grupos.png",width=600,height=800)
> # criando grafico em branco p/ o estado do
> RSmap(xlim=c(-80,-35),ylim=c(-45,10), type="n")map.axes()
> # plotando mapalines(mapaBrasil[,2], mapaBrasil[,1])
> # plotando estaçõespoints( latlonRS[,2], latlonRS[,1], pch=membros)
> # fechando PNGdev.off()
>
>  dados que estou usando são:            V1     V21 -33.51 -53.352 -30.83
>  -55.63 -29.45  -51.84 -29.68 -51.055 -29.05 -50.136 -31.78 -52.417 -32.03
>  -52.18 -30.53 -52.51
> e o erro que aparece é:
>
> Carregando pacotes exigidos: MASSCarregando pacotes exigidos: mapsErro em
> plot.xy(xy.coords(x, y), type = type, ...) :   objeto 'membros' não
> encontrado
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121031/b3ce063b/attachment.html>