Learn R Programming

HDMT (version 1.0.5)

A Multiple Testing Procedure for High-Dimensional Mediation Hypotheses

Description

A multiple-testing procedure for high-dimensional mediation hypotheses. Mediation analysis is of rising interest in epidemiology and clinical trials. Among existing methods for mediation analyses, the popular joint significance (JS) test yields an overly conservative type I error rate and therefore low power. In the R package 'HDMT' we implement a multiple-testing procedure that accurately controls the family-wise error rate (FWER) and the false discovery rate (FDR) when using JS for testing high-dimensional mediation hypotheses. The core of our procedure is based on estimating the proportions of three component null hypotheses and deriving the corresponding mixture distribution of null p-values. Results of the data examples include better-behaved quantile-quantile plots and improved detection of novel mediation relationships on the role of DNA methylation in genetic regulation of gene expression. With increasing interest in mediation by molecular intermediaries such as gene expression, the proposed method addresses an unmet methodological challenge. Methods used in the package refer to James Y. Dai, Janet L. Stanford & Michael LeBlanc (2020) .

Copy Link

Version

Install

install.packages('HDMT')

Monthly Downloads

321

Version

1.0.5

License

MIT + file LICENSE

Maintainer

James Dai

Last Published

January 29th, 2022

Functions in HDMT (1.0.5)

adjust_quantile

A function to compute the quantiles of the estimated mixture null distribution for pmax using either the approximation or exact method
correct_qqplot

A function to draw the corrected quantile-quantile plot for p-max using the expected quantiles
snp_input

An example dataset to demonstrate the usage of 'HDMT'
fdr_est

A function to compute the estimated pointwise FDR for every observed p-max
exercise_input

An example dataset to demonstrate the usage of 'HDMT'
null_estimation

A function to estimate the proportions of the three component nulls
fwer_est

A function used to compute Family wise error rate (FWER) cutoff for p-max at a designated level