From e51cff538833e0e0a765aeb47012f60f80a307a9 Mon Sep 17 00:00:00 2001
From: Walmes Zeviani <walmes@ufpr.br>
Date: Sun, 15 May 2016 20:24:32 -0300
Subject: [PATCH] =?UTF-8?q?Adiciona=20c=C3=B3digo=20para=20fazer=20xyplot?=
=?UTF-8?q?=203x3=20da=20distribui=C3=A7=C3=A3o.?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit
---
vignettes/v04_poisson_generalizada.Rmd | 30 ++++++++++++++++++++++++++
1 file changed, 30 insertions(+)
diff --git a/vignettes/v04_poisson_generalizada.Rmd b/vignettes/v04_poisson_generalizada.Rmd
index e897731..2f6e42a 100644
--- a/vignettes/v04_poisson_generalizada.Rmd
+++ b/vignettes/v04_poisson_generalizada.Rmd
@@ -84,6 +84,36 @@ gamma <- 0
fy <- dpgnz0(y = y, theta = theta, gamma = gamma)
plot(fy ~ y, type = "h", xlab = "y", ylab = "f(y)")
lines(y + 0.3, dpois(y, lambda = theta), type = "h", col = 2)
+
+grid <- expand.grid(lambda = c(2, 8, 15),
+ alpha = c(-0.05, 0, 0.05))
+y <- 0:35
+
+py <- mapply(FUN = dpgnz,
+ lambda = grid$lambda,
+ alpha = grid$alpha,
+ MoreArgs = list(y = y), SIMPLIFY = FALSE)
+grid <- cbind(grid[rep(1:nrow(grid), each = length(y)), ],
+ y = y,
+ py = unlist(py))
+
+useOuterStrips(xyplot(py ~ y | factor(lambda) + factor(alpha),
+ ylab = expression(f(y)),
+ xlab = expression(y),
+ data = grid, type = "h",
+ panel = function(x, y, ...) {
+ m <- sum(x * y)
+ panel.xyplot(x, y, ...)
+ panel.abline(v = m, lty = 2)
+ }),
+ strip = strip.custom(
+ strip.names = TRUE,
+ var.name = expression(lambda == ""),
+ sep = ""),
+ strip.left = strip.custom(
+ strip.names = TRUE,
+ var.name = expression(alpha == ""),
+ sep = ""))
```
## Recursos interativos com o `rpanel` ##
--
GitLab