accumulation: Creates a vector with accumulation counts of transcription...

View source: R/accumulation.R

accumulationR Documentation

Creates a vector with accumulation counts of transcription factors for each chromosome base.

Description

From a dataset with transcription factor (TF) binding regions, this function creates a vector (accvector) in which, for each chromosome base, the accumulation of the TFs present in the input dataset is calculated. Three types of accumulation are considered: TF accumulation, region accumulation and base accumulation. TF accumulation: for each base, it is the number of different TFs present in the neighborhood of the considered base. The neighborhood is defined by a window with half-width w centered on the considered base. Region accumulation: for each base, it is the number of regions containing TFs in the neighborhood of the considered base. If in the neighborhood of a base there are two input binding regions of the same TF, the accumulation value in that base is equal to 2 (differently from the TF accumulation, whose value in the same case is equal to 1). Base accumulation: for each base, it is the total number of bases belonging to input regions containing TFs in the neighborhood of the considered base. With w=0, a single base approach is applied (no base neighborhood is considered). In this case, if in the input dataset overlapping regions for the same TF and chromosome do not exist, the results of TF, region and base accumulation are equal.

Usage

accumulation(data, acctype, chr, w)

Arguments

data

a GRanges object containing coordinates of TF binding regions and their TF name.

acctype

a string with the name of the accumulation type: "TF", "region", "base".

chr

a string with the name of the chromosome (e.g., "chr1"). With chr = "all" all the chromosomes in the input GRanges object are considered.

w

an integer, half-width of the window that defines the neighborhood of each base.

Value

A list of four elements:

accvector

a Rle (or SimpleRleList if chr = "all") object containing the accumulation for each base of the selected chromosome.

acctype

a string with the accumulation type used.

chr

a string with the chromosome name associated with the accumulation vector.

w

an integer with the half-width of the window used to calculate the accumulation vector.

Examples


# loading dataset
data("Ishikawa")
# to calculate TF accumulation for the chromosome 21 for w=0
TF_acc_21_w_0 <- accumulation(Ishikawa,"TF","chr21",0)


DEIB-GECO/TFHAZ documentation built on July 29, 2023, 5:45 p.m.