tenta <div>points(c(-70,-4),pch=20,col="blue")<br><br><div class="gmail_quote">Em 28 de outubro de 2012 09:08, D'Apolito <span dir="ltr"><<a href="mailto:carlosdapolito@gmail.com" target="_blank">carlosdapolito@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="text-indent:0px!important">Olá todos,</div></div><div style="text-indent:0px!important"><br></div><div style="text-indent:0px!important">
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.</div>


<div style="text-indent:0px!important"><br></div><div style="text-indent:0px!important">Att</div><div style="text-indent:0px!important">Carlos</div><div style="text-indent:0px!important"><br></div><div style="text-indent:0px!important">


<br></div><div style="text-indent:0px!important"><br></div><div style="text-indent:0px!important">o código é o seguinte:</div><div style="text-indent:0px!important"><br></div><div style="text-indent:0px!important"><div><br>


</div><div><div><div>elevation.loess = loess(z ~ x*y, data=data.s, degree = 2, span = .3, </div><div>family="gaussian")</div><div><br></div><div>elevation.fit = expand.grid(list(x = seq(65.86,71.6,0.01), </div>

<div>
y = seq(1.81,6.845,0.01)))</div><div><br></div><div>z = predict(elevation.loess, newdata = elevation.fit)</div><div><br></div><div>z1<-ifelse(z<0,0,z)</div></div></div><div><br></div><div><div>height<-(elevation.fit$Height = as.numeric(z1))</div>


</div><div><br></div><div>library(lattice)</div><div><br></div><div>levelplot(height~-x*-y, data = elevation.fit,</div><div> xlab = "Latitude", ylab = "Longitude",</div><div> main = "Surface elevation data",</div>


<div> col.regions = terrain.colors(100),cuts=40,contour=T, plot.new=T, points=T)</div><div><br></div><div><br></div><div>contourplot(height~-x*-y, data = elevation.fit,</div><div>cuts = 20,labels = TRUE,col.regions = terrain.colors(100),</div>


<div>contour = TRUE, pretty = TRUE,region = T)</div><div><br></div><div>points(-70,-4,pch=20,col="blue") # aqui é onde não funciona</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div>


<div><br></div></div><br>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>
<br style="font-family:Tahoma;font-size:13px"><font face="'Times New Roman'" size="4">Natália da Silva Martins<br>Bacharel em Estatística - Universidade Estadual de Maringá/ UEM<br>Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP<br>
Contato: (19) 8306-4743</font> <br><div><br></div><br>
</div>