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.
| 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] 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) | ||
| #------------------------------------------------------------------------------------------ | #------------------------------------------------------------------------------------------ | ||
