parameter <- list(xi = c(2, 0.2))
Lambda <- function(t, xi) (t / xi[2])^xi[1]
priorDensity <- function(xi) dgamma(xi, c(2, 0.2), 1)
start <- parameter
model <- set.to.class("NHPP", parameter, start = start, Lambda = Lambda,
priorDensity = priorDensity)
Run the code above in your browser using DataLab