obtain_peakinfor: Obtain peak sites information.

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

Obtain peak sites information.

Description

Obtain peak sites information from the output of peak_calling.

Usage

obtain_peakinfor(peak_infor_dir)

Arguments

peak_infor_dir

a character specifies the directory of peak calling results by peak_calling function.

Details

This function will output the m6A peak sites informaiton in a data.frame file, which includes the seqnames, start, end, strand, the reads count in all concerned samples of each peak site, the log2Foldchagne and padj for each peak.

Value

peaks_site_infors

A data.frame file, which includes the seqnames, start, end, strand, the reads count in all concerned samples of each peak site, the log2Foldchagne and padj for each peak.

Examples

## Not run: 
Group1_peakinfor <- obtain_peakinfor(peak_infor_dir="./Group1_peakcalling")

## End(Not run)

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