inst/examples/dsc-shrink/methods/mixfdr.wrapper.R

In stephens999/dscr: Dynamic Statistical Comparisons in R

#' @title wrapper for mixfdr for shrinkage DSC
#'
#' @description Runs mixfdr to compute effect estimates from z scores (computed from betahat and sebetahat)  
#' @details None
#' 
#' @param input a list with elements betahat and sebetahat
#' @param args a list containing other additional arguments to ash
#' 
#' @return output a list with the following elements
#' \item{beta_est}{vector containing point estimates for beta}
#' \item{pi0_est}{scalar estimate of pi0}
#'
library(mixfdr)

mixfdr.wrapper=function(input,args=NULL){
  res = try(do.call(mixFdr, args= c(list(x=input$betahat/input$sebetahat, noiseSD=1, plots=FALSE),args)))
  if(inherits(res,"try-error")){res=list(effectSize=rep(NA,length(input$betahat)),pi0_est=NA)}
  return(list(beta_est=res$effectSize*input$sebetahat,pi0_est=res$pi[1]))
}