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 Ambos lados da revisão seguinte | ||
ridiculas [2011/06/01 17:19] walmes [Gráficos do R exportados pelo dispositivo ti//k//z] |
ridiculas [2011/06/07 19:25] walmes |
||
---|---|---|---|
Linha 5: | Linha 5: | ||
**//Ridículas//** é a página do LEG dedicada à fornecer //dicas curtas// sobre R, e.g. condução de análises, operação com dados e confecção de gráficos. As dicas estão organizadas pelo título, seguido de descrição, palavras-chave e CMR (código mínimo reproduzível). Se você deseja contribuir com a nossa página de Ridículas, envie e-mail para ''walmes@ufpr.br''. | **//Ridículas//** é a página do LEG dedicada à fornecer //dicas curtas// sobre R, e.g. condução de análises, operação com dados e confecção de gráficos. As dicas estão organizadas pelo título, seguido de descrição, palavras-chave e CMR (código mínimo reproduzível). Se você deseja contribuir com a nossa página de Ridículas, envie e-mail para ''walmes@ufpr.br''. | ||
+ | |||
+ | ---- | ||
+ | |||
+ | ==== 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> | ||
---- | ---- |