Learn R Programming

survivALL (version 0.9.3)

Continuous Biomarker Assessment by Exhaustive Survival Analysis

Description

In routine practice, biomarker performance is calculated by splitting a patient cohort at some arbitrary level, often by median gene expression. The logic behind this is to divide patients into <80><9c>high<80><9d> or <80><9c>low<80><9d> expression groups that in turn correlate with either good or poor prognosis. However, this median-split approach assumes that the data set composition adheres to a strict 1:1 proportion of high vs. low expression, that for every one <80><9c>low<80><9d> there is an equivalent <80><9c>high<80><9d>. In reality, data sets are often heterogeneous in their composition (Perou, CM et al., 2000 )- i.e. this 1:1 relationship is unlikely to exist and the true relationship unknown. Given this limitation, it remains difficult to determine where the most significant separation should be made. For example, estrogen receptor (ER) status determined by immunohistochemistry is standard practice in predicting hormone therapy response, where ER is found in an ~1:3 ratio (-:+) in the population (Selli, C et al., 2016 ). We would expect therefore, upon dividing patients by ER expression, 25% to be classified <80><9c>low<80><9d> and 75% <80><9c>high<80><9d>, and an otherwise 50-50 split to incorrectly classify 25% of our patient cohort, rendering our survival estimate under powered. 'survivALL' is a data-driven approach to calculate the relative survival estimates for all possible points of separation - i.e. at all possible ratios of <80><9c>high<80><9d> vs. <80><9c>low<80><9d> - allowing a measure<80><99>s relationship with survival to be more reliably determined and quantified. We see this as a solution to a flaw in common research practice, namely the failure of a true biomarker as part of a meta-analysis.

Copy Link

Version

Install

install.packages('survivALL')

Monthly Downloads

15

Version

0.9.3

License

MIT + file LICENSE

Maintainer

Dominic Pearce

Last Published

April 25th, 2018

Functions in survivALL (0.9.3)

checkContSig

Calculate association between survival and a continuous measure#' @inheritParams allPvals
compSelect

Given a dataframe of phenotypic information, use a variable (i.e. a single column) to define a patient subset of given proportion
plotALL

Calculate and combine hazard ratio, pvalue, threshold and area-between-curve data and plot
removeOutliers

Calculate outliers in a numeric vector and then convert those values to NA
areaBetweenCurves

Calculate the area above the bootstrapped thresholds but below the HR distribution and vice versa for each point of separation
survivALL

Calculate and combine hazard ratio, pvalue, threshold and area-between-curve data as a single dataframe
allHR

For all possible separation points for a cohort ordered by a continuous measurement, calculate hazard ratio
allPvals

For all possible separation points for a cohort ordered by a continuous measurement, perform a uni- or multivariate log-rank test
hrSignificance

Calculate HR significance using bootstrap results
nki_subset

NKI breast cancer patients subset with complete t.dmfs and e.dmfs information
bootstrapThresholds

Calculate per-separation point hazard ratio thresholds