IPPP (version 1.1)

IPPPconddens: Conditional Probability Density Given the Location of one Event

Description

Given knowledge of the location of one event, this function determines the probability density function (pdf) of the distribution of the n-th point above/below the known event.

Usage

IPPPconddens(x, pointlocation, xrate, yrate, nthpoint = 1, mode = "forward")

Arguments

x

Points at which the values of the pdf are determined

pointlocation

Location of the known event

xrate

Vector of (strictly increasing) real numbers

yrate

Vector of strictly positive real numbers of the same length as xrate. The vectors xrate and yrate form the rate function r in the sense that r=approxfun(xrate,yrate)

nthpoint

OPTIONAL, default is 1. Setting this to a value i will determine the pdf of the ith event above/below the event at the position pointlocation.

mode

OPTIONAL, default is "forward". Determines whether the pdf for points above(right) the known event is determined (mode="forward"), or whether the pdf for points below(left) the known event is determined (mode="backward")

Value

Returns a list containing two entries:

x

Duplicate of the input of the same name

densval

Vector consisting of the values of the conditional probability density, evaluated at x

Details

The value of the rate function r below min(xrate) is set to r(min(xrate)), and for values above max(xrate) it is set to r(max(xrate)). Both r(min(xrate)) and r(max(xrate)) need to be strictly positive for the results to be correct.

References

Hohmann, Niklas. "Conditional Densities and Simulations of Inhomogeneous Poisson Point Processes: The R package "IPPP"" arXiv 2019. <arXiv:1901.10754>

See Also

IPPPnthpointdens for the probability density of the location of the of the n-th event, given that a fixed number of events occur.

vignette("IPPP") for an overview of the features of the IPPP package and some background.

Examples

# NOT RUN {
#define rate function
xrate=seq(0,2*pi,length.out=1000)
yrate=sin(xrate)+1.01
plot(xrate,yrate,type="l",main="Rate Function")
#define known event:
pointlocation=3
lines(c(3,3),c(0,3),lwd=3)
#values where the density is calulated:
x=xrate
#simulate the pdf of the event above the known event:
r1=IPPPconddens(x,pointlocation,xrate,yrate)
lines(r1$x,r1$densval,col="red",lwd=3)
#simulate the pdf of the event below the known event:
r2=IPPPconddens(x,pointlocation,xrate,yrate,mode="backward")
lines(r2$x,r2$densval,col="blue",lwd=3)
#simulate the pdf of the third event below the known event:
r2=IPPPconddens(x,pointlocation,xrate,yrate,mode="backward",nthpoint=3)
lines(r2$x,r2$densval,col="green",lwd=3)
# }