
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

Tem alguma coisa em http://www.leg.ufpr.br/doku.php/pessoais:elias:cursos:uem2012 Há outros materiais procure por spatial analysis R no google (tem muita coisa) On 07/14/2013 11:19 AM, Natalia Martins wrote:
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
_______________________________________________ R-br mailing list R-br@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.
participantes (2)
-
Elias Krainski
-
Natalia Martins