ridiculas - Wiki do LEG

Não foi possível enviar o arquivo. Será algum problema com as permissões?

Você está aqui: start » ridiculas

Diferenças

Aqui você vê as diferenças entre duas revisões dessa página.

--- ridiculas [2011/05/31 22:51]
walmes
+++ ridiculas [2011/07/05 14:13]
walmes [section 2]
@@ Linha 8: / Linha 8: @@
 ----
-==== Análise de dados de proporção usando modelo lineares generalizados ====
+==== 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
+#------------------------------------------------------------------------------------------
+</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.\\
+palavras-chave: #tikz, #sweave, #latex.
+<code>
+\documentclass{article}
+\usepackage{Sweave}
+\usepackage{tikz}
+\SweaveOpts{keep.source=true}
+\title{Usando ti\textit{k}z no Sweave}
+\author{Walmes Zeviani\\ LEG/UFPR}
+\begin{document}
+\maketitle
+C\'{o}digo m\'{i}nimo reproduz\'{i}vel usando \texttt{tikzDevice} para exportar gr\'{a}ficos feitos no R
+para codifica\c{c}\~{a}o ti\textit{k}z.
+{\footnotesize
+<<results=hide>>=
+#------------------------------------------------------------------------------------------
+require(tikzDevice)
+set.seed(2011); x <- rnorm(200)
+tikz("plot.tex", w=5, h=3)
+hist(x, freq=FALSE, ylab="Densidade",
+     main="Histograma de uma amostra de $X \\sim N(\\mu=0, \\sigma^2=1)$")
+curve(dnorm(x), col=2, add=TRUE, lwd=2); rug(x)
+legend("topleft", col=2, lty=1, lwd=2, bty="n",
+       legend="$\\displaystyle \\frac{1}{\\sqrt{2\\pi\\sigma^2}}\\cdot e^{-\\frac{(x-\\mu)^2}{2\\sigma^2}}$")
+box(); dev.off()
+#------------------------------------------------------------------------------------------
+@
+}
+\input{plot.tex}
+\end{document}
+</code>
+----
+==== Gráfico com dois eixos coordenados ====
+Temporariamente sem descrição.\\
+palavras-chave: #série, #eixo.
+<code R>
+#------------------------------------------------------------------------------------------
+#                                                               por Ivan, Benilton e Walmes
+#------------------------------------------------------------------------------------------
+# dados de séries de dados indexadas no tempo (meses)
+lines <- 'meses  temp umidade    rad   chuva
+    Jan    26.49   86.58    795.88  0.36
+    Fev    26.65   88.49    710.24  0.34
+    Mar    27.19   86.16    772.99  0.21
+    Abr    26.28   89.75    574.88  0.67
+    Mai    26.62   89.22    614.02  0.31
+    Jun    26.13   87.83    680.08  0.26
+    Jul    25.83   86.57    675.97  0.15
+    Ago    27.05   83.14    756.44  0.07
+    Set    27.60   83.02    925.57  0.14
+    Out    27.44   85.16    927.71  0.17
+    Nov    26.56   88.18    788.87  0.19
+    Dez    25.87   90.63    703.94  0.33'
+da <- read.table(textConnection(lines), header=TRUE)
+str(da)
+#------------------------------------------------------------------------------------------
+# gráfico de duas séries no mesmo gráfico com um eixo coordenado para cada
+par(mar=c(5,5,4,5))
+plot(da$temp, type="b", pch=15, lwd=1.5,
+     xlab="Meses (Anos 2010)", main="Temperatura e umidade em 2010",
+     ylab=expression(Temperatura~group("(", degree*C, ")")),
+     ylim=c(24,30), xlim=c(1,12), axes=FALSE)
+axis(1, at=1:12, labels=da$meses)
+axis(2)
+par(new=TRUE)
+plot(da$umidade, type="b", pch=14, axes=FALSE, frame=TRUE, ann=FALSE)
+axis(4)
+mtext(text="Umidade (%)", 4, line=3)
+legend("topleft", bty="n", seg.len=3, lty=1, lwd=c(1.5,1), pch=c(15,14),
+       legend=c("Temperatura","Umidade"), merge=TRUE, trace=FALSE)
+#------------------------------------------------------------------------------------------
+# funções sugeridas
+apropos("month")
+demo(plotmath)
+#------------------------------------------------------------------------------------------
+</code>
+----
+==== Análise de dados de proporção usando modelo linear generalizado ====
 Temporariamente sem descrição.\\

Diferenças

Navegação

Busca

Ferramentas

QR Code