poisson (version 1.0)

nhpp.plot: Plot simulated non-homogeneous Poisson processes

Description

Plot num.events simulated non-homogeneous Poisson processes, plus the mean and quantiles

Usage

nhpp.plot(rate, num.events, prob.func, num.sims = 100, t0 = 0, t1 = NULL,  num.points = 100, quantiles = c(0.025, 0.975), ...)

Arguments

rate

the rate at which events occur in the Poisson process, aka lambda

num.events

the number of event times to simulate in each process

prob.func

function that takes time as sole argument and returns value between 0 and 1

num.sims

number of simulated paths to plot

t0

start time

t1

end time

num.points

number of points to use in estimating mean and quantile processes

quantiles

plot these quantile processes

...

further arguments to be passed to or from methods

Value

x

Matrix of event times, one process per column

x.bar

Vector of mean process event times

x.q

Matrix of quantile event times, one process per column

See Also

nhpp.scenario, hpp.plot

Examples

intensity <- function(t) pmin(t/3, 1)
nhpp.plot(rate = 5, num.events = 20, num.sims = 100, main='My simulated NHPPs',
          prob.func=intensity)