R/step_cdf.R

In jean997/sherlockAsh: CAUSE: Causal Analysis Using Summary Effect Estimates

#'@export
step_cdf <- function(x, start, stop, p){
    ix <- which(start < x &  stop <= x)
    c <- sum(p[ix])
    if(any(start <x &  stop > x)){
       if(sum(start <x &  stop > x) > 1) stop("error")
       ix <- which(start <x &  stop > x)
       prop <- (x-start[ix])/(stop[ix]-start[ix])
       c <- c + p[ix]*prop
    }
    return(c)
}