R/lme_mass_FDR.R

In Deep-MI/fslmer: Univariate And Mass-Univariate Linear Mixed-Effects Analysis For FreeSurfer Imaging Data

#' Conduct FDR thresholding
#'
#' @param p a vector of p-values
#' @param fdr threshold
#'
#' @return
#' The function returns an updated threshold for uncorrected p-values that keeps the FDR.
#'
#' @export
#'
#' @examples
#' \dontrun{pthFDR <- lme_mass_FDR(pvals, 0.05)}

lme_mass_FDR<-function(p,fdr)
{
    p<-sort(abs(p))
    nv<-length(p)

    nn<-matrix(c(1:nv),ncol=1)

    if (length(which(p<=(fdr*nn/nv)))>0)
    {
        imax<-max(which(p<=(fdr*nn/nv)))
        pthresh<-p[imax]
    }
    else
    {
        pthresh<-min(p)/10
    }

    return(pthresh)
}