An R package for identification of differentially methylated regions (DMRs) from high density chip data.
- Order the CpG sites based on the chromosome number and chromosome coordinates. Calculate the adjacent CpG sites’ correlation in the same chromosome one by one. Furthermore, get the clustered CpG sites that share similar information based on the correlation and treat these clusters as candidate DMRs.
- Generalized beta regression is used to model the CpG sites in the same cluster; parameter estimation is performed via maximum likelihood function (ML).
- Conduct a log likelihood ratio test to the interesting phenotype after adjusting for the covariables, then summarize the result and report the DMRs
library(devtools)
install_github("chengzhouwu/gbdmr")
To start using the package just load it as any other package.
library(maxLik)
library(gbdmr)
library(doParallel)
library(foreach)
library(IlluminaHumanMethylation450kanno.ilmn12.hg19)
library(minfi)
library(geograbi) ## https://github.com/yousefi138
library(limma)
library(corrplot)
For running the algorithm one has to have three objects:
- a matrix with methylation values;
- a phenotype dataset (the IDs of the two data sets need to correspond one-to-one);
- specific correlation value from 0 to 1.
All the work is done by one command gbdmr
and for more information use ?gbdmr
in R.
A comprehensive usage protocol for gbdmr has been documented and is available on RPubs: https://rpubs.com/Chengzhou/1139894
This detailed guide provides step-by-step instructions for leveraging the functionalities of gbdmr effectively.
The output generated by gbdmr is versatile, allowing users to conduct various downstream analyses. Some of the notable analyses include enrichment analysis, mediation analysis, association studies, and prediction studies. The flexibility of gbdmr output empowers researchers to delve deeper into the biological and statistical aspects of their data, making it a valuable tool for comprehensive data exploration and interpretation.
Wu, Chengzhou, Xichen Mou, and Hongmei Zhang. "Gbdmr: identifying differentially methylated CpG regions in the human genome via generalized beta regressions." BMC bioinformatics 25.1 (2024): 97.