annotateChipSeqPeaks: Function to match each peaks to associated genes in terms of...

In uc-bd2k/CLEAN: Clustering Enrichment Analysis

Function to match each peaks to associated genes in terms of RefID

Description

This function is intended to search all peaks within up and down distances for each gene. Users could specify their own platforms or the function will download platform from UCSC as default.

Usage

annotateChipSeqPeaks(chip.seq,transcriptDB = NULL, species = "Hs", distanceRange = c(-1e+06,1e+06), file = NULL)

Arguments

chip.seq	A matrix list all CHIP-Seq peaks in which the chromosome ID for all peaks are listed in the first column, peaks' start positions are listed in the second column, peaks' end position are listed in the third column, peaks' ID are listed in the forth column and the fifth column lists peaks' score.
file	File name with path where the final matrix will deposit. If it is null the final matrix will be returned from the function.
transcriptDB	Platform table name with path saved in the local matchine.
distanceRange	Chip.seq peaks found within distances from Transcription Start Site.

Details

This function is intended to search all peaks within up and down distances for each gene.

Value

This function returns a matrix in which the first column contains the RefSeqID, the second column contains Entrez ID , the third column contains chromosome ID, the fourth column contains the absolute distance from each peak to transcription start site and the fifth column contains chip-seq score.

Author(s)

Zhen Hu

References

Publication forthcoming.

Examples

data(nMyc)
RefTable<-createRefGenome("mm8_refGene",genome.id = "mm8", tablename = "refGene")
ChipSeq<-annotateChipSeqPeaks(n.Myc,transcriptDB=RefTable,species="Mm",distanceRange=c(-1e+06,1e+06),file=NULL)
res<-chipSeqWeightedSum(ChipSeq,verbose=T)

uc-bd2k/CLEAN documentation built on Oct. 19, 2024, 11:45 p.m.