Você está aqui: start » cursos » mct » exemplo4

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 anteriorRevisão anterior
cursos:mct:exemplo4 [2011/05/25 18:41] – paulojuscursos: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)