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

Ambos lados da revisão anterior Revisão anterior
cursos:mct:exemplo4 [2011/05/25 18:41]
paulojus
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 ​aleatoriamente +teste <- merge(alfab[sample(1:​nrow(alfab),​ nrow(alfab)/​2),​],​ # metade ​aleatóriamente 
-               ​renda[sample(1:​nrow(renda),​ nrow(renda)/​2),​],​ # metade ​aleatoriamente+               ​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 54: Linha 78:
  
 par(mfrow=c(1,​2)) par(mfrow=c(1,​2))
-with(ar, {hist(alfabetizados,​ prob=T); rug(alfabetizados)lines(density(alfabetizados))}) +with(ar, {hist(alfabetizados,​ prob=T); rug(alfabetizados)lines(density(alfabetizados))}) 
-with(ar, {hist(log(renda),​ prob=T); rug(leg(renda))lines(density(log(renda)))})+with(ar, {hist(log(renda),​ prob=T); rug(leg(renda))lines(density(log(renda)))})
 par(mfrow=c(1,​1)) par(mfrow=c(1,​1))
  
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)