[R-br] Mapas tematicos
Natalia Martins
nsmbarreto em gmail.com
Domingo Julho 14 11:19:15 BRT 2013
Prezados, bom dia.
Gostaria de saber como posso trabalhar com mapas tematicos no R, pois
trabalhava com o software geoda e neste conseguia criar mapas tematicos,
calcular o indice de moran, LISA, entre outros.
Encontrei na internet comandos para calcular o indice de moran, mas a
construção dos mapas tematicos não consegui.
Segue os comando que estou utilizando:
## lendo dados tipo shapefile
brasil <- readShapePoly("G:/natalia/BRASIL.shp", IDvar="UF")
plot(brasil)
head(brasil em data) #para ler a primeira linha como cabeçalho
brasil em data$NOME <- as.character(brasil em data$UF)
dados <- read.table("G:/natalia/dados.csv", sep=";", head=T, enc="latin1")
head(dados)
UF Estado Analfabetismo Homicidio Renda
1 AC Acre 15.19 21.55768 663.2902
2 AL Alagoas 22.52 59.18695 378.0110
3 AM Amapa 7.89 30.48163 521.1194
4 AP Amazonas 9.60 26.99392 485.6758
5 BA Bahia 15.39 36.51580 463.8209
6 CE Ceara 17.19 25.34000 427.8397
#junta o shape com os dados pela varivel comum ind
ind <- match(brasil em data$UF, dados$UF)
dados <- dados[ind,]
row.names(dados) <- brasil$UF
#linca dados com formato espacial
brasil <- spCbind(brasil, dados)
names(brasil)
[1] "UF" "ESTADO" "REGIAO" "NOME" "UF.1"
[6] "Estado" "Analfabetismo" "Homicidio" "Renda"
head(brasil em data)
UF ESTADO REGIAO NOME UF.1 Estado Analfabetismo Homicidio Renda
AC AC Acre NO AC AC Acre 15.19 21.55768 663.2902
AL AL Alagoas NE AL AL Alagoas 22.52 59.18695 378.0110
AM AM Amazonas NO AM AM Amapa 7.89 30.48163 521.1194
AP AP Amap NO AP AP Amazonas 9.60 26.99392 485.6758
BA BA Bahia NE BA BA Bahia 15.39 36.51580 463.8209
CE CE Cear NE CE CE Ceara 17.19 25.34000 427.8397
Att, Natalia
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130714/01a7d509/attachment.html>
Mais detalhes sobre a lista de discussão R-br