R/pval.R

In SharonLutz/UmediationThread: Runs the Umediation Function

mediate_pval <- function(x, xhat){
  if(is.null(xhat)){
    stop("xhat is null")
  }
  if(is.na(xhat)){
    stop("xhat is NA")
  }
  ## Compute p-values
  if(xhat == 0){
    out <- 1
  } else {
    out <- 2 * min(sum(x > 0), sum(x < 0)) / length(x)
  }
  return(min(out, 1))
}