calibrationBeltIntersections: Calibration Belt Significant Deviations

In givitiR: The GiViTI Calibration Test and Belt

Description

calibrationBeltIntersections returns the intervals where the calibration belt significantly deviates from the bisector.

Usage

calibrationBeltIntersections(cbBound, seqP, minMax)

Arguments

cbBound	A data.frame object with the numeric variables "U" and "L", representing the upper and lower boundary of the calibration belt.
seqP	The vector of the the probabilities where the points of the calibration belt have been evaluated.
minMax	A list with two elements, named min and max, representing the minimum and maximum probabilities in the model under evaluation.

Value

A list with two components, overBisector and underBisector. Each component is a list containing all the intervals where the calibration belt is significantly over/under the bisector.

See Also

givitiCalibrationBelt and plot.givitiCalibrationBelt to compute and plot the calibaration belt, and givitiCalibrationTest to perform the associated calibration test.

Examples

e <- runif(1000)
logite <- logit(e)
eMod <- logistic(logit(e) +  (logit(e))^2)
o <- rbinom(1000, size = 1, prob = eMod)
data <- data.frame(e = e, o = o, logite = logite)

seqP <- seq(from = .01, to =.99, by = .01)
seqG <- logit(seqP)

minMax <- list(min = min(e), max = max(e))

fwLR <- polynomialLogRegrFw(data, .95, 4, 1)
cbBound <- calibrationBeltPoints(data, seqG, fwLR$m, fwLR$fit, .95, .90, "external")
calibrationBeltIntersections(cbBound, seqP, minMax)

Example output

$overBisector
$overBisector[[1]]
[1] 0.0009061052 0.3358427050

$overBisector[[2]]
[1] 0.6625993 0.9998250


$underBisector
list()

givitiR documentation built on May 2, 2019, 10:58 a.m.