Skip to content
Snippets Groups Projects
Commit 99b7075b authored by Walmes Marques Zeviani's avatar Walmes Marques Zeviani
Browse files

Finaliza os slides da Gamma-Count.

parent c4631ba9
No related branches found
No related tags found
No related merge requests found
Pipeline #
...@@ -76,7 +76,7 @@ set_parent("slides-mrdcr.Rnw") ...@@ -76,7 +76,7 @@ set_parent("slides-mrdcr.Rnw")
\framebreak \framebreak
\item A distribuição acumulada do tempo até $\vartheta_{n}$ é \item A distribuição acumulada do tempo até $\vartheta_{n}$ é
\item $$F_n(T) = \Pr(\vartheta_n \leq T) = \int_{0}^{T} $$F_n(T) = \Pr(\vartheta_n \leq T) = \int_{0}^{T}
\frac{\beta^{n\alpha}}{\Gamma(n\alpha)}\cdot t^{n\alpha-1}\cdot \frac{\beta^{n\alpha}}{\Gamma(n\alpha)}\cdot t^{n\alpha-1}\cdot
\exp\{-\beta t\}\,\text{d}t.$$ \exp\{-\beta t\}\,\text{d}t.$$
\item Seja $[0,T)$ um intervalo e $N_{T}$ a v.a. número de eventos \item Seja $[0,T)$ um intervalo e $N_{T}$ a v.a. número de eventos
...@@ -102,7 +102,8 @@ set_parent("slides-mrdcr.Rnw") ...@@ -102,7 +102,8 @@ set_parent("slides-mrdcr.Rnw")
&= \left[ \int_{0}^{T} &= \left[ \int_{0}^{T}
\frac{\beta^{n\alpha}}{\Gamma(n\alpha)} \frac{\beta^{n\alpha}}{\Gamma(n\alpha)}
t^{n\alpha-1}\cdot t^{n\alpha-1}\cdot
\exp\{-\beta t\}\, \text{d}t \right] - \exp\{-\beta t\}\, \text{d}t \right] \\
&\quad -
\left[ \int_{0}^{T} \left[ \int_{0}^{T}
\frac{\beta^{(n+1)\alpha}}{\Gamma((n+1)\alpha)} \frac{\beta^{(n+1)\alpha}}{\Gamma((n+1)\alpha)}
t^{(n+1)\alpha-1}\cdot t^{(n+1)\alpha-1}\cdot
...@@ -111,6 +112,40 @@ set_parent("slides-mrdcr.Rnw") ...@@ -111,6 +112,40 @@ set_parent("slides-mrdcr.Rnw")
\end{itemize} \end{itemize}
\end{frame} \end{frame}
\begin{frame}[fragile]
<<fig.width = 9, fig.height = 4.5, out.width = "0.95\\textwidth">>=
library(MRDCr)
grid <- expand.grid(lambda = c(2, 8, 15),
alpha = c(0.5, 1, 2.5))
y <- 0:30
py <- mapply(FUN = dgcnt,
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 = ""))
@
\end{frame}
\begin{frame}[allowframebreaks] \begin{frame}[allowframebreaks]
\frametitle{Parametrização para modelo de regressão} \frametitle{Parametrização para modelo de regressão}
\begin{itemize} \begin{itemize}
...@@ -167,39 +202,6 @@ llgcnt ...@@ -167,39 +202,6 @@ llgcnt
@ @
\end{frame} \end{frame}
\begin{frame}[fragile]
<<fig.width = 9, fig.height = 4.5, out.width = "0.95\\textwidth">>=
grid <- expand.grid(lambda = c(2, 8, 15),
alpha = c(0.5, 1, 2.5))
y <- 0:30
py <- mapply(FUN = dgcnt,
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 = ""))
@
\end{frame}
\begin{frame}{Estudos de caso} \begin{frame}{Estudos de caso}
{\it Vignette} \href{run:../doc/v06_gamma_count.html}{\tt {\it Vignette} \href{run:../doc/v06_gamma_count.html}{\tt
gamma\_count.html} gamma\_count.html}
...@@ -211,8 +213,3 @@ useOuterStrips(xyplot(py ~ y | factor(lambda) + factor(alpha), ...@@ -211,8 +213,3 @@ useOuterStrips(xyplot(py ~ y | factor(lambda) + factor(alpha),
\item[\tt cambras]: Gols do Campeonato Brasileiro de 2010. \item[\tt cambras]: Gols do Campeonato Brasileiro de 2010.
\end{description} \end{description}
\end{frame} \end{frame}
\begin{frame}
\frametitle{Conclusões}
TODO
\end{frame}
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment