map_peak_longTX: Map peak sites to the longest transcript

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

Map peak sites to the longest transcript

Description

Map m6A peak sites to the longest transcirpt for each gene

Usage

map_peak_longTX(filepath,
                annotation_file,
                peak_sites_infor)

Arguments

filepath	A character specifies the directory of the peak calling result obtained by peak_calling function.
annotation_file	A character specifies the directory toward a gene annotation GFF/GTF file.
peak_sites_infor	A data.frame file gives the m6A peak sites information obtained by obtain_peakinfor function.

Value

A list containing two components:

mapped_peakGRList	peak sites mapped to the longest transicript in GRangesList format
mapped_peankinfor	peak sites information mapped to the longest transcirpt of each gene in data.frame format

Examples

## Not run: 
Group1_mappeak_LTX <- map_peak_longTX(filepath="./Group1_peakcalling",
                                     annotation_file=GENE_ANNO_GTF,
                                     peak_sites_infor=Group1_peakinfor)

## End(Not run)

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