################################################################################
### Part of the surveillance package, http://surveillance.r-forge.r-project.org
### Free software under the terms of the GNU General Public License, version 2,
### a copy of which is available at http://www.r-project.org/Licenses/.
###
### Temporal interaction functions for twinstim's epidemic component.
### Specific implementations are in seperate files (e.g.: exponential, step).
###
### Copyright (C) 2009-2014 Sebastian Meyer
### $Revision: 733 $
### $Date: 2014-01-31 17:46:47 +0100 (Fre, 31 Jan 2014) $
################################################################################
#####################
### "Constructor" ###
#####################
tiaf <- function (g, G, deriv, Deriv, npars, validpars = NULL)
{
npars <- as.integer(npars)
if (length(npars) != 1 || npars < 0L) {
stop("'tiaf'/'npars' must be a single nonnegative number")
}
haspars <- npars > 0L
g <- .checknargs3(g, "tiaf$g")
G <- .checknargs3(G, "tiaf$G")
if (!haspars || missing(deriv)) deriv <- NULL
if (!haspars || missing(Deriv)) Deriv <- NULL
if (!is.null(deriv)) deriv <- .checknargs3(deriv, "tiaf$deriv")
if (!is.null(Deriv)) Deriv <- .checknargs3(Deriv, "tiaf$Deriv")
validpars <- if (!haspars || is.null(validpars))
NULL else match.fun(validpars)
list(g = g, G = G,
deriv = deriv, Deriv = Deriv,
npars = npars, validpars = validpars)
}
#################################
### Constant temporal interaction
#################################
tiaf.constant <- function ()
{
res <- list(
g = as.function(alist(t=, pars=, types=, rep.int(1, length(t))), envir = .GlobalEnv),
G = as.function(alist(t=, pars=, types=, t), envir = .GlobalEnv),
npars = 0L
)
attr(res, "constant") <- TRUE
res
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.