Não foi possível enviar o arquivo. Será algum problema com as permissões?
Diferenças
Aqui você vê as diferenças entre duas revisões dessa página.
| Próxima revisão Ambos lados da revisão seguinte | |||
|
cursos:mct:exemplo3 [2011/05/05 07:03] paulojus criada |
cursos:mct:exemplo3 [2011/05/05 07:15] paulojus |
||
|---|---|---|---|
| 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> | ||
| + | |||
