Não foi possível enviar o arquivo. Será algum problema com as permissões?
Diferenças

Diferenças

Aqui você vê as diferenças entre duas revisões dessa página.

Link para esta página de comparações

Próxima revisão
Revisão anterior
cursos:mct:exemplo3 [2011/05/05 07:03]
paulojus criada
cursos:mct:exemplo3 [2013/02/20 23:32] (atual)
paulojus a revisão anterior foi restaurada (2011/05/05 07:15)
Linha 8: Linha 8:
 </​code>​ </​code>​
  
- +Primeiro vamos obter um arquivo de mapa do [[http://​www.ibge.gov.br|IBGE]]:​ 
-O mapa foi carregado do site do IBGE em  +  O mapa foi carregado do site do IBGE em  
- <​Geociências>​ + <​Mapeamento das unidades territoriais>​ +  +   ​<​Geociências>​ + <​Mapeamento das unidades territoriais>​ + <​Produtos>​ + <Malha municipal digital 2007> (clicar no disquete de 2007)  
- <​Produtos>​ + <Malha municipal digital 2007> (clicar no disquete de 2007)  +   ​Na nova janela:
- Na nova janela+
    ​Escolha E500, E100 ou E2500    ​Escolha E500, E100 ou E2500
    ​Escolha uma projeção .Usado: Proj_Geografica    ​Escolha uma projeção .Usado: Proj_Geografica
Linha 21: Linha 20:
 Salve os arquivos com extensão .shp, .shx e .dbf Salve os arquivos com extensão .shp, .shx e .dbf
  
-### carregando ​pacotes maptools e sp+Carregando ​pacotes maptools e sp 
 +<code R>
 require(maptools) require(maptools)
 require(sp) require(sp)
 +</​code>​
  
-### lendo o mapa (disponivel em leg.ufpr.br/​~paulojus/​mct/​dados )+Lendo o mapa no  
 +(Site do IBGE pode ser demorado. Os [[http://leg.ufpr.br/​~paulojus/​mct/​dados|arquivos 12mu500gc* podem ser compiados daqui]]) 
 +<code R>
 ac <- readShapePoly("​12mu500gc"​) ​ ac <- readShapePoly("​12mu500gc"​) ​
  
Linha 43: Linha 46:
 dim(ac@data) dim(ac@data)
 ac@data[1:​3,​] ac@data[1:​3,​]
 +</​code>​
 +
 +  Obtencao dos dados do  [[http://​www.datasus.gov.br|datasus]]
 +   - clicar em <​Informacoes de Saude> ​
 +   - clicar em <​Demograficas e Socioeconomicas> ​  
 +   - Clicar em Abastecimento de agua - 1991, 2000
 +   - Selecione Brasil por Municipio ​
 +   - Montar a tabela assim:
 +   - Linha: Muncipio <Todas as categorias>​
 +     - Coluna: Abastecimento de Agua
 +     - Conteudo: Moradores
 +   - Clicar em Formato: <Tabela com bordas> e <​Mostra>​ (no final da pg.)
 +   - Após os dados serem mostrados, clicar em <Copia com o .CSV> e salvar o arquivo
  
-#### +Alternativamente aqui os [[http://www.leg.ufpr.br/​~paulojus/​mct/​dados/​A104241200_17_213_55.csv|dados de abastecimento de agua]], do arquivo 
-### obtencao dos dados de alfabetizacao +obtidos ​a partir do site do DATASUS ​
-### - no site do datasus ​http://www.datasus.gov.br +
-### - clicar em <​Informacoes de Saude>  +
-### - clicar em <​Demograficas e Socioeconomicas>​  +
-### - Clicar em Abastecimento de agua - 1991, 2000 +
-### - Selecione Brasil por Municipio  +
-### - Montar a tabela assim: +
-###    - Linha: Muncipio <Todas as categorias>​ +
-###    - Coluna: Abastecimento de Agua +
-###    - Conteudo: Moradores +
-### - Clicar em Formato: <Tabela com bordas> e <​Mostra>​ (no final da pg.) +
-### - Apos os dados serem mostrados, clicar em <Copia com o .CSV> e salvar o arquvio +
-### lendo dados de abastecimento de agua, do arquivo +
-### A104241200_17_213_55.csv salvo a partir do site do DATASUS ​ +
-### (disponivel em leg.ufpr.br/​~paulojus/​mct/​dados )+
  
 +<code R>
 agua <- read.csv2("​A104241200_17_213_55.csv",​ skip=3) agua <- read.csv2("​A104241200_17_213_55.csv",​ skip=3)
  
Linha 68: Linha 71:
 dim(agua) dim(agua)
 agua[1:3,] agua[1:3,]
 +</code
  
-### Necessário ​colocar dados na mesma ordem do mapa+Combinando as bases. Notar que é necessário ​colocar dados na mesma ordem do mapa.
  
 +<code R>
 ### extraindo o codigo de municipios dos dados ### extraindo o codigo de municipios dos dados
 codmund <- substr(as.character(agua[,​1]),​ 1, 6) codmund <- substr(as.character(agua[,​1]),​ 1, 6)
Linha 110: Linha 115:
 legend("​topright",​ leglabs(c(0,​12,​25,​1),​ "<",​ ">"​),​ legend("​topright",​ leglabs(c(0,​12,​25,​1),​ "<",​ ">"​),​
        ​fill=cores3,​ bty="​n",​ cex=2)        ​fill=cores3,​ bty="​n",​ cex=2)
 +</​code>​
 +
 +Uma análise:
  
 +<code R>
 ### carrega pacote spdep ### carrega pacote spdep
 require(spdep) require(spdep)
Linha 123: Linha 132:
 eb <- EBImoran.mc(odat[,​10],​ odat[,​2]+odat[,​3],​ nbw, 999) eb <- EBImoran.mc(odat[,​10],​ odat[,​2]+odat[,​3],​ nbw, 999)
 eb eb
- 
-## HELP !!! 
-help(EBImoran.mc) ​ ## ou simplesmente ?​EBImoran.mc 
-## tipos: "​text",​ "​html",​ "​postscript",​ "​ps",​ "​pdf"​. ​ 
 </​code>​ </​code>​
-### o mapa foi carregado do site do IBGE em  
-### <​Geociências>​ + <​Mapeamento das unidades territoriais>​ +  
-### <​Produtos>​ + <Malha municipal digital 2007> (clicar no disquete de 2007)  
-### Na nova janela 
-### Escolha E500, E100 ou E2500 
-### Escolha uma projeção .Usado: Proj_Geografica 
-### Escolha um formato. Usado: ArcView_shp 
-### Escolha um nível territorial. ​ 
-### Escolha a unidade territorial.  ​ 
  
-### salve os arquivos. Pelo menos aqueles com extensão .shp, .shx e .dbf+HELP !!!
  
-### carregando pacotes maptools e sp +<code R>
-require(maptools) +
-require(sp) +
- +
-### lendo o mapa (disponivel em leg.ufpr.br/​~paulojus/​mct/​dados ) +
-ac <- readShapePoly("​12mu500gc"​)  +
- +
-### visualizando +
-plot(ac) +
- +
-### atributos dos municipios (variáveis) +
-names(ac) +
- +
-### atributos/​elementos do objeto +
-names(attributes(ac)) +
- +
-### bounding box +
-bbox(ac) +
- +
-### inspecionando os atributos (dados) dos municípios +
-dim(ac@data) +
-ac@data[1:​3,​] +
- +
-#### +
-### obtencao dos dados de alfabetizacao +
-### - no site do datasus http://​www.datasus.gov.br +
-### - clicar em <​Informacoes de Saude +
-### - clicar em <​Demograficas e Socioeconomicas>​  +
-### - Clicar em Abastecimento de agua - 1991, 2000 +
-### - Selecione Brasil por Municipio  +
-### - Montar a tabela assim: +
-###    - Linha: Muncipio <Todas as categorias>​ +
-###    - Coluna: Abastecimento de Agua +
-###    - Conteudo: Moradores +
-### - Clicar em Formato: <Tabela com bordas> e <​Mostra>​ (no final da pg.) +
-### - Apos os dados serem mostrados, clicar em <Copia com o .CSV> e salvar o arquvio +
-### lendo dados de abastecimento de agua, do arquivo +
-### A104241200_17_213_55.csv salvo a partir do site do DATASUS  +
-### (disponivel em leg.ufpr.br/​~paulojus/​mct/​dados ) +
- +
-agua <- read.csv2("​A104241200_17_213_55.csv",​ skip=3) +
- +
-system("​file -bi A104241200_17_213_55.csv"​) +
-## encondings em diferentes sistemas! +
- +
-dim(agua) +
-agua[1:​3,​] +
- +
-### Necessário colocar dados na mesma ordem do mapa +
- +
-### extraindo o codigo de municipios dos dados +
-codmund <- substr(as.character(agua[,​1]),​ 1, 6) +
-codmund[1:​3] +
- +
-### extraindo o código do mapa com apenas 6 digitos +
-ac@data[1:​3,​] +
-cod6map <- substr(as.character(ac@data[,​1]),​ 1, 6) +
-head(cod6map,​ 3) +
- +
-### selecionando os dados do ACRE +
-### e preservando a ordem do mapa +
-odat <- merge(data.frame(cod=cod6map),​ +
-              data.frame(cod=codmund,​ agua[,-1]), sort=FALSE) +
-dim(odat) +
-names(odat) +
-table(odat$cod==cod6map) ### verifica ordem +
- +
-### adiciona os dados ao SpatialPolygons +
-ac@data <- data.frame(ac@data,​ odat[,​-1]) +
- +
-### proporcao de rede geral +
-ac@data$prg <- (odat[,​2]+odat[,​3])/​odat[,​10] +
-spplot(ac, "​prg"​) +
-spplot(ac, "​prg",​ col.regions=heat.colors(16)) +
- +
-### visualiza no mapa +
- +
-### criando 3 categorias +
-cat3 <- findInterval(ac$prg,​ c(0, .12, .25, 1)) +
-table(cat3) +
- +
-### três cores +
-cores3 <- c("​red",​ "​orange",​ "​yellow"​) +
- +
-### plotando o mapa temático +
-par(mar=c(0,​0,​0,​0)) +
-plot(ac, col=cores3[cat3]) +
-legend("​topright",​ leglabs(c(0,​12,​25,​1),​ "<",​ ">"​),​ +
-       ​fill=cores3,​ bty="​n",​ cex=2) +
- +
-### carrega pacote spdep +
-require(spdep) +
- +
-args(EBImoran.mc) +
-nb <- poly2nb(ac) +
-nb +
- +
-plot(nb, coordinates(ac),​ add=T) +
-nbw <- nb2listw(nb) +
- +
-eb <- EBImoran.mc(odat[,​10],​ odat[,​2]+odat[,​3],​ nbw, 999) +
-eb +
- +
-## HELP !!!+
 help(EBImoran.mc) ​ ## ou simplesmente ?​EBImoran.mc help(EBImoran.mc) ​ ## ou simplesmente ?​EBImoran.mc
 ## tipos: "​text",​ "​html",​ "​postscript",​ "​ps",​ "​pdf"​. ​ ## tipos: "​text",​ "​html",​ "​postscript",​ "​ps",​ "​pdf"​. ​
 </​code>​ </​code>​
 +

QR Code
QR Code cursos:mct:exemplo3 (generated for current page)