get.intersection: get.intersection Obtain the intersection of two distributions using the kernel method

Description

Obtain the intersection of two distributions using the kernel method. Warning: This function does not check the parameters ref and test.

Usage

get.intersection(ref, test, model = c("kernel", "binormal", "ordinal"), ...)

Arguments

ref

The reference standard. A column in a data frame or a vector indicating the classification by the reference test. The reference standard must be coded either as 0 (absence of the condition) or 1 (presence of the condition)

test

The index test or test under evaluation. A column in a dataset or vector indicating the test results on a continuous scale.

model

The model used for estimating the intersection(s). Default = 'kernel'.

...

passing arguments to the kernel density function, other than kernel='gaussian' (default).

Value

A vector of points of intersection, ordered on their density. The tail has the highest density.

References

Landsheer, J. A. (2016). Interval of Uncertainty: An Alternative Approach for the Determination of Decision Thresholds, with an Illustrative Application for the Prediction of Prostate Cancer. PloS One, 11(11), e0166007.

Examples

Run this code

# NOT RUN {
ref=c(rep(0,500), rep(1,500))
test=c(rnorm(500,0,1), rnorm(500,1,2))
(get.intersection(ref, test)) # two intersections! Generates warning in other functions!
# }