R/validate.lambda.summaryData.R

In zhangh12/ARTP3: Pathway and Gene-Level Association Test

validate.lambda.summaryData <- function(summary.files, lambda){
  
  if(is.null(lambda)){
    lambda <- rep(1.0, length(summary.files))
  }
  
  if(!is.vector(lambda)){
    msg <- 'lambda should be a numeric vector'
    stop(msg)
  }
  # each summary file should has one inflation factor
  if(length(summary.files) != length(lambda)){
    msg <- 'Each summary file should has one inflation factor'
    stop(msg)
  }
  
  lambda
  
}