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@data)
#para ler a primeira linha como cabeçalho
brasil@data$NOME
<- as.character(brasil@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@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@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