R/hrSignificance.R

Defines functions hrSignificance

Documented in hrSignificance

#' Calculate HR significance using bootstrap results
#' @param hratios Hazard ratio vector, output of allHR()
#' @param bs_dfr Dataframe of bootstrapped hazard ratios
#' @return A per-separation point vector equal to the distance beyond or within
#' the bootstrapped thresholds
hrSignificance <- function(hratios, bs_dfr){
    if(all(bs_dfr == 0)){
        bsp <- NA
    } else {
        n <- ncol(bs_dfr)
        means <- apply(bs_dfr, 1, function(x) mean(x, na.rm = TRUE))
        obs_min_exp <-  abs(hratios - means)
        excd_up <- means + obs_min_exp
        excd_down <- means - obs_min_exp

        up_sums <- rowSums(bs_dfr > excd_up, na.rm = TRUE)
        down_sums <- rowSums(bs_dfr < excd_down, na.rm = TRUE)
        bsp <- round((up_sums + down_sums) / n, 10)
        bsp[is.na(hratios)] <- NA
    }
    bsp
}

Try the survivALL package in your browser

Any scripts or data that you put into this service are public.

survivALL documentation built on May 1, 2019, 11:02 p.m.