Nothing
STIKhat <- function(xyt, s.region, t.region, dist, times, lambda, correction = TRUE, infectious = TRUE)
{
require(splancs)
if (missing(s.region)) s.region <- sbox(xyt[,1:2],xfrac=0.01,yfrac=0.01)
if (missing(t.region)) t.region <- range(xyt[,3],na.rm=T)
pts <- xyt[,1:2]
xytimes <- xyt[,3]
ptsx <- pts[, 1]
ptsy <- pts[, 2]
ptst <- xytimes
npt <- length(ptsx)
ndist <- length(dist)
dist <- sort(dist)
ntimes <- length(times)
times <- sort(times)
bsupt <- max(t.region)
binft <- min(t.region)
area <- areapl(s.region)*(bsupt-binft)
np <- length(s.region[, 1])
polyx <- c(s.region[, 1], s.region[1, 1])
polyy <- c(s.region[, 2], s.region[1, 2])
hkhat <- array(0, dim = c(ndist,ntimes))
if(missing(lambda))
{
misl <- 1
lambda <- rep(1,npt)
}
else misl <- 0
if (length(lambda)==1) lambda <- rep(lambda,npt)
if (correction==T){edg <- 1} else edg <- 0
if (infectious==T){infd <- 1} else infd <- 0
storage.mode(hkhat) <- "double"
if (area>10) lambda <- lambda*area
nev <- rep(0,ntimes)
klist <- .Fortran("stikfunction", as.double(ptsx),
as.double(ptsy), as.double(ptst),
as.integer(npt), as.double(polyx),
as.double(polyy), as.integer(np),
as.double(dist), as.integer(ndist),
as.double(times), as.integer(ntimes),
as.double(bsupt), as.double(binft),
as.double(lambda), as.integer(infd),
as.integer(edg), (hkhat))
khat <- klist[[17]]
if (misl==1)
{
if (area>10)
hkhat <- ((area^3)/(npt*(npt-1)))*khat
else
hkhat <- (area/(npt*(npt-1)))*khat
}
else
{
hkhat <- khat
if (area>10)
hkhat <- hkhat*area
else
hkhat <- hkhat/area
}
if(dist[1]==0) hkhat[1,]=0
if(times[1]==0) hkhat[,1]=0
Khpp <- matrix(0,ncol=length(times),nrow=length(dist))
for(i in 1:length(dist)){Khpp[i,]<-pi*(dist[i]^2)*times}
if (infectious==TRUE) Ktheo <- Khpp
else Ktheo <- 2*Khpp
lhat <- hkhat-Ktheo
invisible(return(list(Khat=hkhat,Ktheo=Ktheo,dist=dist,times=times,infectious=infectious)))
}
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.