DLR: Estimating the Capacity for Improvement in Diagnostic Risk...

In DTComPair: Comparison of Binary Diagnostic Tests in a Paired Study Design

Estimating the Capacity for Improvement in Diagnostic Risk Prediction with an additional marker based on the Diagnostic Likelihood Ratio (DLR)

Description

This function allows for estimating the log diagnostic likelihood ratio in a regression model approach. It can be used to assess the gain in diagnostic accuracy for a new binary or continuous diagnostic marker compared to established markers, to determine the impact of covariates on the risk prediction model, and to estimate the DLR for selected marker/covariate values.

Usage

DLR(basemodel, augmentedmodel, diseasestatus, dataset, clustervar = NULL, alpha=0.05)

Arguments

basemodel	pre-test/base model X, formula character string
augmentedmodel	post-test/ augmented model V, formula character string, this is usually the basemodel X including the additional diagnostic test of interest Y and interactions XY
diseasestatus	variable name containing disease status, assumed to be a 0/1 variable, for having condition of interest (1) or not (0), character string
dataset	dataframe, needs to be in wide format with one observation per subject
clustervar	optional, cluster variable name in dataset, character string
alpha	significance level alpha used for confidence intervals, the default is 0.05.

Details

This function is an implementation of the algorithm described in the appendix of Gu and Pepe (2009) using the GEE approach in order to get standard error estimates. The definition of I and Zero matrices is slightly more flexible than the ones described in section 3 in order to allow for models without interaction.

Value

Returns a list including

logPreTestModel	logistic regression model output for prior disease using base model X: P(D=1|X). All estimates are on a log scale.
logPostTestModel	logistic regression model output for posterior disease using augmented model V: P(D=1|X,Y),i.e. P(D=1|V). All estimates are on a log scale.
logDLRModel	regression model output for log DLR defined as difference between logPostTestModel and logPreTestModel. All estimates are on a log scale.
DLR	Positive/negative DLR for diagnostic marker Y, with all base covariates X set to 1. Results are only sensible for binary marker Y taking values 0/1.

Author(s)

Thomas Hielscher ([e-mail])

References

Gu, W. and Pepe, M. S. (2009). Estimating the capacity for improvement in risk prediction with a marker. Biostatistics, 10(1):172-86.

See Also

dlr.regtest

Examples

library(DTComPair)
data(Paired1)

# test y1 conditioned on null model: DLR+(Y1=1) and DLR-(Y1=0)

DLR("~ 1","~ y1","d",Paired1)

# test y1 conditioned on test y2 with interaction, DLR+(Y1=1|Y2=1) and DLR-(Y1=0|Y2=1)

DLR("~ y2","~ y2 * y1","d",Paired1) 

DTComPair documentation built on Sept. 24, 2024, 5:08 p.m.