From cc0e0468465cfc20a25d06dec8871f01df8b60c1 Mon Sep 17 00:00:00 2001
From: Eduardo Junior <edujrrib@gmail.com>
Date: Tue, 28 Jun 2016 03:07:37 -0300
Subject: [PATCH] =?UTF-8?q?Adiciona=20fun=C3=A7=C3=A3o=20para=20c=C3=A1lcu?=
=?UTF-8?q?lo=20da=20constante=20Z?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit
---
docs/_setup.R | 17 ++++++++++++++++-
1 file changed, 16 insertions(+), 1 deletion(-)
diff --git a/docs/_setup.R b/docs/_setup.R
index 0f04474..4575413 100644
--- a/docs/_setup.R
+++ b/docs/_setup.R
@@ -45,7 +45,7 @@ opts_chunk$set(
##======================================================================
## Configura opções de output no documento
-options(digits = 3, OutDec = ",",
+options(digits = 3, OutDec = ",", scipen=-1,
xtable.caption.placement = "top",
xtable.include.rownames = FALSE,
xtable.booktabs = TRUE)
@@ -101,6 +101,21 @@ fonte <- function(texto, side = 1, line = -1, adj = 0,
invisible()
}
+##======================================================================
+## Calcula Z para um c(lambda, phi)
+funZ <- function(lambda, nu, maxit = 500, tol = 1e-5) {
+ z <- rep(NA, maxit)
+ j = 1
+ ##
+ z[j] <- exp(j * log(lambda) - nu * lfactorial(j))
+ ##
+ while (abs(z[j] - 0) > tol && j <= maxit) {
+ j = j + 1
+ z[j] <- exp(j * log(lambda) - nu * lfactorial(j))
+ }
+ return(cbind("j" = 0:j, "z" = c(1, z[!is.na(z)])))
+}
+
##======================================================================
## Para padronizar as analises de deviance de todos os modelos
## -
--
GitLab