gene_methy_level_distdecay: Obtain gene's methylation level

In NWPU-903PR/m6AexpressBHM: m6Aexpress-BHM: Predicting m6A regulation of gene expression in multiple-groups context by a Bayesian Hierarchical Mixture model

Obtain gene's methylation level

Description

Calculate methylation level of each gene, which is weighted by the distance of peak center to stop-codon

Usage

gene_methy_level_distdecay(mapLTX_peakinfor,
                           size_factor,
                           peak_dist_stopcodon)

Arguments

mapLTX_peakinfor	A list containing peak sites information mapped to the longest transcript of gene.
size_factor	A numberic vector containing the sequencing depth size factors for each MeRIP-seq samples, which are save in the peak_calling output directory in a "size_factor.Rdata" file
peak_dist_stopcodon	A data.frame specifing the distance of peak center to stop cond and the site of each peak center.

Value

It will return a data.frame file containing the methylation level weighted by the distance of peak center to sop codon for each gene in each concerned sample.

Examples

## Not run: 
load("./Group1_peakcalling/size_factor.Rdata")
Group1_sizefactor <- size_factor
Group1_methylevel_distdecay <- gene_methy_level_distdecay(mapLTX_peakinfor=Group1_mappeak_LTX,
                                                          peak_dist_stopcodon=Group1_peakcenter2stopcondon,
                                                          size_factor=Group1_sizefactor)


## End(Not run)

NWPU-903PR/m6AexpressBHM documentation built on May 29, 2022, 11:07 p.m.