findpeakcenter: Obtain peak center

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

Obtain peak center

Description

Obtain each peak center site based on transcriptome coordinate.

Usage

findpeakcenter(annotation_file,
              maplongTX_peak)

Arguments

annotation_file	A character specifies the directory toward a gene annotation GFF/GTF file.
maplongTX_peak	A list includes the information of peak sites mapped to the longest transcript of each gene.

Value

It will return a data.frame file including the sites of peak center for each peak.

Examples

## Not run: 
Group1_peakcenter <- findpeakcenter(annotation_file=GENE_ANNO_GTF,
                                    maplongTX_peak=Group1_mappeak_LTX)

## End(Not run)

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