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:exemplo4 [2011/05/05 07:41]
paulojus criada
cursos:mct:exemplo4 [2011/05/27 01:03] (atual)
paulojus
Linha 2: Linha 2:
  
 <code R> <code R>
 +-
 ------------------------------------------------------------------ ------------------------------------------------------------------
 # dados de renda por município no ano de 2000 (.xls) # dados de renda por município no ano de 2000 (.xls)
  
 +## lendo direto do arquivo excel
 require(gdata) # contém a função read.xls() para importar planilhas require(gdata) # contém a função read.xls() para importar planilhas
 renda <- read.xls("​http://​www.leg.ufpr.br/​~paulojus/​mct/​renda2000municipios.xls",​ renda <- read.xls("​http://​www.leg.ufpr.br/​~paulojus/​mct/​renda2000municipios.xls",​
                   sheet=1, fileEncoding='​latin1',​ encoding="​latin1"​)                   sheet=1, fileEncoding='​latin1',​ encoding="​latin1"​)
 +names(renda) <- tolower(names(renda))
 +names(renda)[4] <- "​renda"​
 +str(renda)
 +
 +## lendo arquivo csv exportado do excel
 +renda <- read.csv("​http://​www.leg.ufpr.br/​~paulojus/​mct/​renda2000municipios.csv",​ fileEnc="​latin1",​
 +                  sep="​\t",​ dec=","​)
 +renda
 +head(renda)
 names(renda) <- tolower(names(renda)) names(renda) <- tolower(names(renda))
 names(renda)[4] <- "​renda"​ names(renda)[4] <- "​renda"​
Linha 15: Linha 26:
 # dados de alfabetizados por municÃípio no ano de 2000 (.xls) # dados de alfabetizados por municÃípio no ano de 2000 (.xls)
  
 +## lendo direto do arquivo excel
 alfab <- read.xls("​http://​www.leg.ufpr.br/​~paulojus/​mct/​alfabetizados2000municipios.xls",​ alfab <- read.xls("​http://​www.leg.ufpr.br/​~paulojus/​mct/​alfabetizados2000municipios.xls",​
                   sheet=1, fileEncoding='​latin1',​ encoding="​latin1"​)                   sheet=1, fileEncoding='​latin1',​ encoding="​latin1"​)
 +names(alfab) <- tolower(names(alfab))
 +names(alfab)[4] <- "​alfabetizados"​
 +str(alfab)
 +
 +## lendo arquivo csv exportado do excel
 +
 +## opcoes usuais para leitura/​importaçao dos dados: read.csv(), read.csv2() e read.table()
 +## dados em: http://​www.leg.ufpr.br/​~paulojus/​mct/​alfabetizados2000municipios.csv
 +file.show("​dados/​alfabetizados2000municipios.csv"​)
 +alfab <- read.csv("​http://​www.leg.ufpr.br/​~paulojus/​mct/​alfabetizados2000municipios.csv",​
 +                   ​fileEncoding='​latin1',​ sep="​\t",​ dec=","​)
 +alfab
 names(alfab) <- tolower(names(alfab)) names(alfab) <- tolower(names(alfab))
 names(alfab)[4] <- "​alfabetizados"​ names(alfab)[4] <- "​alfabetizados"​
Linha 32: Linha 56:
 str(teste) str(teste)
  
-teste <- merge(alfab[sample(1:​nrow(alfab),​ nrow(alfab)/​2),​],​ # metade ​aleatóriamente +teste <- merge(alfab[sample(1:​nrow(alfab),​ nrow(alfab)/​2),​],​ # metade ​aleatóriamente 
-               ​renda[sample(1:​nrow(renda),​ nrow(renda)/​2),​],​ # metade ​aleatóriamente+               ​renda[sample(1:​nrow(renda),​ nrow(renda)/​2),​],​ # metade ​aleatóriamente
                ​by.x="​codigo",​ by.y="​codigo"​)                ​by.x="​codigo",​ by.y="​codigo"​)
 str(teste) str(teste)
Linha 47: Linha 71:
  
 #​------------------------------------------------------------------------------------------ #​------------------------------------------------------------------------------------------
-qual a distribuição de frequência das variáveis+# distribuição de frequência das variáveis
  
 hist(ar$alfabetizados,​ freq=FALSE) hist(ar$alfabetizados,​ freq=FALSE)
Linha 103: Linha 127:
  
 #​------------------------------------------------------------------------------------------ #​------------------------------------------------------------------------------------------
-# Opção 2+# Opção 2:  
 + 
 +## graficos de dispsersao
  
 require(lattice) require(lattice)
 xyplot(alfabetizados~log(renda)|sigla,​ data=ar) xyplot(alfabetizados~log(renda)|sigla,​ data=ar)
 +
 +##
 +xyplot(alfabetizados~renda,​ data=ar)
 +xyplot(alfabetizados~log10(renda),​ data=ar)
 +xyplot(alfabetizados~renda,​ data=ar, scales=list(x=list(log=TRUE)))
 +
 +xyplot(alfabetizados~log10(renda)|sigla,​ data=ar)
 +xyplot(alfabetizados~log10(renda)|sigla,​ data=ar, type=c("​p","​r"​))
 +xyplot(alfabetizados~log10(renda)|sigla,​ data=ar,
 +       ​type=c("​p","​r"​),​ cex=0.2)
 +
 +
 +#​------------------------------------------------------------------------------------------
 +# box plot
 +bwplot(alfabetizados~sigla,​ data=ar)
 +bwplot(alfabetizados~sigla,​ data=ar, pch="​|"​)
 +bwplot(alfabetizados~sigla,​ data=ar, pch="​|",​
 +       ​scales=list(x=list(rot=90)))
 +bwplot(alfabetizados~sigla,​ data=ar, pch="​|",​
 +       ​fill="​red",​ scales=list(x=list(rot=90)))
  
 #​------------------------------------------------------------------------------------------ #​------------------------------------------------------------------------------------------
Linha 119: Linha 165:
 densityplot(~log(renda)|sigla,​ data=ar) densityplot(~log(renda)|sigla,​ data=ar)
 densityplot(~alfabetizados|sigla,​ data=ar) densityplot(~alfabetizados|sigla,​ data=ar)
 +
 +# criando 5 categorias de renda
 +
 +ar$r.cat <- equal.count(log(ar$renda),​ number=5, overlap=0.1)
 +str(ar)
 +
 +bwplot(alfabetizados~sigla|r.cat,​ data=ar, pch="​|"​)
 +bwplot(alfabetizados~r.cat|sigla,​ data=ar, pch="​|",​ fill="​gray70"​)
 +
 +histogram(~alfabetizados|r.cat,​ data=ar)
 +densityplot(~alfabetizados|r.cat,​ data=ar)
  
 #​------------------------------------------------------------------------------------------ #​------------------------------------------------------------------------------------------

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