Hélio e demais colegas, boa tarde!
Tenho aqui um exemplo que pode ser adaptado ao caso em questão:
### <code r>
require(geoR)
attach(parana)
plot(borders, xlab = "X Coord", ylab = "Y Coord ", type = "n", asp=1)
polygon(borders)
dNum <- 5
dBrk <- as.vector(quantile(data, probs = seq(0, 1, by = (1/dNum))))
dCut <- cut(data, breaks=dBrk, include.l=T, labels=F)
points(coords, pch=(1:dNum)[dCut], col=c(4,3,7,2))
# points(coords, pch=(15:18)[dCut], col=c(4,3,7,2))
### posicionamento relativo
legend('topright', 'coords', paste0("Q",1:dNum), pch=(1:dNum)[dCut],
col=c(4,3,7,2), cex=.5)
legend('bottomleft', 'coords', paste0("Q",1:dNum), pch=(1:dNum)[dCut],
col=c(4,3,7,2), horiz=T, cex=.5)
# c("bottomright", "bottom", "bottomleft", "left", "topleft", "top",
"topright", "right", "center")
### posicionamento absoluto
# locator(1)
legend(150, 540, paste0("Q",1:dNum), pch=(1:dNum)[dCut], col=c(4,3,7,2),
horiz=T, cex=.5)
detach()
### </code>
Estou na 'correria' com os experimentos de campo, então o código não está
devidamente comentado. Qualquer dúvida, retorne na lista.
Éder Comunello <c <[email protected]>[email protected]>
Dourados, MS - [22 16.5'S, 54 49'W]
_______________________________________________
R-br mailing list
[email protected]
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.
