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 Próxima revisão | Revisão anterior Próxima revisão Ambos lados da revisão seguinte | ||
|
ridiculas [2011/06/01 17:18] walmes |
ridiculas [2011/07/05 14:18] walmes [Como fazer a justaposição de vários data.frames] |
||
|---|---|---|---|
| Linha 8: | Linha 8: | ||
| ---- | ---- | ||
| - | ==== Gráficos do R exportados pelo dispositivo ti//k//z ==== | + | ==== Como fazer a justaposição de vários data.frames ==== |
| + | |||
| + | Temporariamente sem descrição.\\ | ||
| + | palavras-chave: #merge, #Reduce. | ||
| + | |||
| + | <code R> | ||
| + | #------------------------------------------------------------------------------------------ | ||
| + | # por Walmes | ||
| + | #------------------------------------------------------------------------------------------ | ||
| + | |||
| + | id <- 1:30 # número único que identifica os registros | ||
| + | n <- 20 # número de registros por data.frame | ||
| + | |||
| + | a1 <- data.frame(id=sample(id, n), v1=rnorm(n)) # resposta 1 | ||
| + | a2 <- data.frame(id=sample(id, n), v2=rpois(n,10)) # resposta 2 | ||
| + | a3 <- data.frame(id=sample(id, n), v3=runif(n)) # resposta 3 | ||
| + | |||
| + | merge(a1, a2, by="id") # justapõe 2 data.frames de cada vez | ||
| + | |||
| + | a0 <- list(a1, a2, a3) # cria uma lista com todos os data.frames | ||
| + | |||
| + | |||
| + | Reduce(function(x, y) merge(x, y, by="id"), a0, accumulate=FALSE) # justapõe todos | ||
| + | |||
| + | #------------------------------------------------------------------------------------------ | ||
| + | # font: http://rwiki.sciviews.org/doku.php?id=tips:data-frames:merge | ||
| + | # https://stat.ethz.ch/pipermail/r-help/2008-April/160836.html | ||
| + | # http://econometricsense.blogspot.com/2011/01/merging-multiple-data-frames-in-r.html | ||
| + | # http://www.youtube.com/watch?v=E4uR5I1uLFM | ||
| + | |||
| + | #------------------------------------------------------------------------------------------ | ||
| + | </code> | ||
| + | |||
| + | ---- | ||
| + | |||
| + | |||
| + | ==== Gráfico de valores observados e curva de valores preditos para modelo linear generalizado ==== | ||
| + | |||
| + | Temporariamente sem descrição.\\ | ||
| + | palavras-chave: #glm, #poisson, #predict. | ||
| + | |||
| + | <code R> | ||
| + | #------------------------------------------------------------------------------------------ | ||
| + | # por Walmes | ||
| + | #------------------------------------------------------------------------------------------ | ||
| + | |||
| + | da <- expand.grid(trat=gl(2,1), tempo=1:20) | ||
| + | y <- rpois(nrow(da), lambda=da$tempo/5) | ||
| + | |||
| + | g0 <- glm(y~trat/(tempo+I(tempo^2)), data=da, family=poisson) | ||
| + | |||
| + | new <- expand.grid(trat=gl(2,1), tempo=seq(1,20,l=50)) | ||
| + | new$p0 <- predict(g0, newdata=new, type="response") | ||
| + | |||
| + | plot(y~tempo, da, col=da$trat) | ||
| + | with(subset(new, trat=="1"), lines(p0~tempo, col=1)) | ||
| + | with(subset(new, trat=="2"), lines(p0~tempo, col=2)) | ||
| + | |||
| + | #------------------------------------------------------------------------------------------ | ||
| + | </code> | ||
| + | |||
| + | ---- | ||
| + | |||
| + | |||
| + | ==== Gráficos de barras com intervalos de confiança para as médias ==== | ||
| + | |||
| + | Temporariamente sem descrição.\\ | ||
| + | palavras-chave: #intervalo, #erro_padrão, #barras. | ||
| + | |||
| + | <code R> | ||
| + | #------------------------------------------------------------------------------------------ | ||
| + | # por Walmes | ||
| + | #------------------------------------------------------------------------------------------ | ||
| + | # cmr para colocar intervalos de confiança para as medias num gráfico de barras | ||
| + | |||
| + | da <- data.frame(trat=gl(5,8)) | ||
| + | da$y <- as.numeric(da$trat)+rnorm(nrow(da)) | ||
| + | |||
| + | m0 <- lm(y~trat, da) | ||
| + | |||
| + | new <- data.frame(trat=levels(da$trat)) | ||
| + | new$pred <- predict(m0, newdata=new, interval="confidence") | ||
| + | str(new) | ||
| + | |||
| + | ylim <- c(0, max(new$pred)*1.05) | ||
| + | bp <- barplot(new$pred[,1], ylim=ylim) | ||
| + | arrows(bp, new$pred[,2], bp, new$pred[,3], code=3, angle=90) | ||
| + | box() | ||
| + | |||
| + | #------------------------------------------------------------------------------------------ | ||
| + | # outras referências | ||
| + | |||
| + | browseURL("http://addictedtor.free.fr/graphiques/graphcode.php?graph=54") | ||
| + | browseURL("http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=72") | ||
| + | |||
| + | #------------------------------------------------------------------------------------------ | ||
| + | </code> | ||
| + | |||
| + | ---- | ||
| + | |||
| + | ==== Gráficos do R exportados pelo dispositivo tikz ==== | ||
| Temporariamente sem descrição.\\ | Temporariamente sem descrição.\\ | ||
