diffPeakSummary(chipseq)
diffPeakSummary()所属R语言包:chipseq
A function to identify and produce summary statistics for differentially expressed peaks.
一个功能识别并产生汇总统计差异表达的高峰。
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Given two sets of peaks, this function combines them and summarizes the individual coverage vectors under the combined peak set.
由于两套峰,此功能相结合,并总结下个人保险向量合并峰集。
用法----------Usage----------
diffPeakSummary(ranges1, ranges2,
viewSummary = list(sums = viewSums, maxs = viewMaxs))
参数----------Arguments----------
参数:ranges1
First set of peaks (typically an RleViewsList).
第一组峰(通常是一个RleViewsList)。
参数:ranges2
Second set of peaks (typically an RleViewsList).
第二组峰(通常是一个RleViewsList)。
参数:viewSummary
A list of the per peak summary functions.
每峰汇总函数的列表。
值----------Value----------
A data.frame with one row for each peak in the combined data. The chromosome, start and stop nucleotide positions (+ strand) are given as are the summary statistics requested.
一个data.frame一排在合并后的数据为每个高峰。染色体,启动和停止核苷酸位置(股),以及所要求的汇总统计。
作者(S)----------Author(s)----------
D. Sarkar
举例----------Examples----------
data(cstest)
library(BSgenome.Mmusculus.UCSC.mm9)
seqlevels(cstest) <- seqlevels(Mmusculus)
seqlengths(cstest) <- seqlengths(Mmusculus)
## find peaks[#找到峰]
findPeaks <- function(reads) {
reads.ext <- resize(reads, width = 200)
slice(coverage(reads.ext), lower = 8)
}
peakSummary <- diffPeakSummary(findPeaks(cstest$gfp), findPeaks(cstest$ctcf))
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|