Nothing
lambda <- function(t, x, y, param, object)
{
if (!is.numeric(t) || !is.numeric(x) || !is.numeric(y))
stop(paste("Arguments", sQuote(t), ",", sQuote(x), "and", sQuote(y), "must be numeric."))
if (length(t) != length(x) || length(t) != length(y) || length(x) != length(y))
stop(paste("Arguments", sQuote(t), ",", sQuote(x), "and", sQuote(y), "must be of the same length."))
spatstat.geom::verifyclass(object, "catalog")
revents <- object$revents
storage.mode(revents) <- "double"
theta <- sqrt(param)
out <- numeric(length(t))
for (i in 1:length(t))
{
if (t[i] < revents[1,1])
out[i] <- 0
else
{
out[i] <- .Call("clambdax", as.double(t[i]), as.double(x[i]),
as.double(y[i]), as.double(theta),
revents, PACKAGE="ETAS")[[1]]
}
}
return(out)
}
lambdatemporal <- function(t, fit)
{
obj <- fit$object
cxxlambdtemp(t, fit$param, obj$revents, obj$rpoly,
obj$rtperiod, fit$integ0, fit$ndiv, fit$mver)
}
lambdaspatial <- function(x, y, fit)
{
obj <- fit$object
cxxlambspat(x, y, fit$param, obj$revents, obj$rpoly,
obj$rtperiod, fit$bwd, fit$mver)
}
timetransform <- function(fit)
{
obj <- fit$object
cxxtimetrans(fit$param, obj$revents, obj$rpoly, obj$rtperiod,
fit$integ, fit$ndiv, fit$mver)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.