Nothing
expoTreeSurvival <- function(pars,times,ttypes,shifts=NULL,vflag=0,
return.full=FALSE,rescale=TRUE,
root.lineages=0)
{
if (! is.matrix(pars)) pars <- matrix(as.numeric(pars),nrow=1)
matDim <- dim(pars)
nshifts <- length(shifts) + sum(ttypes==3)
ttypes[ttypes != 3] <- 99
if (matDim[2] < 5) {
cat("Minimum 5 columns required in pars.\n")
return(-Inf)
}
if (matDim[1] < nshifts+1) {
cat("Not enough parameters for shifts.\n")
return(-Inf)
}
if (matDim[1] > nshifts+1) pars <- pars[1:(nshifts+1),]
if (! is.null(shifts)) {
times <- c(times,shifts)
ttypes <- c(ttypes,rep(3,length(shifts)))
}
o <- order(times)
times <- times[o]
ttypes <- ttypes[o]
p <- .Call("expoTreeEval",parameters=pars,
times=as.numeric(times),
ttypes=as.integer(ttypes),
survival=as.integer(c(1,vflag,rescale,root.lineages)))
# p = (p0,p1,p2,...,pN)
if (return.full) {
return(p)
} else {
return(p[2])
}
}
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.