eval_hrz {highriskzone}

Evaluation of the high-risk zone.

Description

Evaluation of the high-risk zone, which is only possible with simulated or thinned data or if the locations of the unobserved events have been revealed..

Usage

eval_hrz(hrz, unobspp, obspp = NULL)

Arguments

hrz

High-risk zone of class owin based on a binary mask (see area.owin)

unobspp

Unobserved spatial point process

obspp

Observed spatial point process

Value

An object of class "hrzeval", which is a list of

numbermiss

number of unobserved events outside the high-risk zone

numberunobserved

number of events in the unobserved point pattern

missingfrac

fraction of unobserved events outside the high-risk zone (numbermiss/numberunobserved)

arearegion

area of the high-risk zone

numberobs

number of events in the observed point pattern

out

subset of the unobserved events, which are outside the high-risk zone

insd

subset of the unobserved events, which are inside the high-risk zone

See Also

inside.owin, area.owin

Examples

 data(craterB)
 # thin data
 set.seed(100)
 thdata <- thin(craterB, nxprob=0.1)
 
 # determine hrz for the "observed events"
 hrz <- det_hrz(thdata$observed, type = "dist", criterion = "area", cutoff = 1500000, nxprob = 0.1)
 
 # evaluate the hrz
 evaluation <- eval_hrz(hrz = hrz$zone, unobspp = thdata$unobserved, obspp = thdata$observed)
 evaluation$missingfrac
 
 op <- par(mar=c(1, 4, 1, 6) , xpd=TRUE)
 plot(evaluation, hrz = hrz, obspp = thdata$observed, plothrz = TRUE, plotobs = TRUE, 
 insidecol = "magenta", outsidecol = "magenta", obscol = "blue", insidepch = 1, 
 outsidepch = 19, main = "Evaluation visualized")
 legend(2400, 2456.4061, c("observed", "unobs inside", "unobs outside"), 
 col = c("blue", "magenta", "magenta"), yjust=1, pch=c(1, 1, 19), cex=0.8)
 par(op)