[R-br] levelplot e contourplot

D'Apolito carlosdapolito em gmail.com
Domingo Outubro 28 09:35:49 BRST 2012


obrigado, Natália, mas ainda não deu:

Erro em plot.xy(xy.coords(x, y), type = type, ...) :
  plot.new ainda não foi chamado







Em 28 de outubro de 2012 08:15, Natalia Martins <nsmbarreto em gmail.com>escreveu:

> tenta
> points(c(-70,-4),pch=20,col="blue")
>
> Em 28 de outubro de 2012 09:08, D'Apolito <carlosdapolito em gmail.com>escreveu:
>
>> Olá todos,
>>
>> Estou produzindo mapas de nível topográfico com o Lattice, usando as
>> funções levelplot() e contourplot(). Depois de pronto, tento colocar pontos
>> no mapa, usando points() e não funciona. Se alguém souber onde está o erro,
>> agradeço qualquer ajudar.
>>
>> Att
>> Carlos
>>
>>
>>
>> o código é o seguinte:
>>
>>
>> elevation.loess = loess(z ~ x*y, data=data.s, degree = 2, span = .3,
>> family="gaussian")
>>
>> elevation.fit = expand.grid(list(x = seq(65.86,71.6,0.01),
>>  y = seq(1.81,6.845,0.01)))
>>
>> z = predict(elevation.loess, newdata = elevation.fit)
>>
>> z1<-ifelse(z<0,0,z)
>>
>> height<-(elevation.fit$Height = as.numeric(z1))
>>
>> library(lattice)
>>
>> levelplot(height~-x*-y, data = elevation.fit,
>>  xlab = "Latitude", ylab = "Longitude",
>>  main = "Surface elevation data",
>>  col.regions = terrain.colors(100),cuts=40,contour=T, plot.new=T,
>> points=T)
>>
>>
>> contourplot(height~-x*-y, data = elevation.fit,
>> cuts = 20,labels = TRUE,col.regions = terrain.colors(100),
>> contour = TRUE, pretty = TRUE,region = T)
>>
>> points(-70,-4,pch=20,col="blue") # aqui é onde não funciona
>>
>>
>>
>>
>>
>>
>>
>>
>> _______________________________________________
>> 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.
>>
>
>
>
> --
>
> 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 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.
>



-- 
_____________________________________
*Carlos D'Apolito - PhD candidate - ID:1239480*
Earth Sciences
School of Geography & Environmental Sciences
Aston Webb Building
University of Birmingham
Birmingham
B15 2TT
UK

phone: 44 - 07 919 673 173
CXD280 em bham.ac.uk
CV: http://lattes.cnpq.br/4557754746424026
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121028/a1dbf02d/attachment.html>


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