methPercent(charm)
methPercent()所属R语言包:charm
Estimate percentage DNA methylation from log-ratios
log比估计的百分比从DNA甲基化
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Estimate percentage DNA methylation from log-ratios
log比估计的百分比从DNA甲基化
用法----------Usage----------
methPercent(m, pmIndex, ngc, commonParams = TRUE)
参数----------Arguments----------
参数:m
a matrix of M-values (methylation log-ratios). One column per sample.
矩阵的m值(甲基数的比率)。每样一列。
参数:pmIndex
A vector of probe indices to use in the calculation. Usually set to the indices of the pm probes (excluding background and other non-specific controls) by using pmIndex=pmindex(dat)
在计算中使用的一个探针指数向量。通常设置为下午探针指数(不包括背景和其他非特定的控件)使用pmIndex = pmindex(DAT)
参数:ngc
a vector with GC-content of probes. Same length as nrow(m)
与GC含量探针的向量。作为NROW相同长度(M)
参数:commonParams
boolean indicating whether a common set of parameters should be used for all samples when converting M-values to percentage methylation.
布尔值,指示是否一套通用的参数应为所有样品转换m值的百分比甲基化时使用。
Details
详情----------Details----------
This function estimates percentage DNA methylation from normalized methylation log-ratios (M-values).
这个功能估计百分比归甲基log比率(M值)的DNA甲基化。
值----------Value----------
a matrix of percentage methylation estimates. Same dimensions as m
百分比甲基估计矩阵。相同尺寸为m
作者(S)----------Author(s)----------
Martin Aryee <aryee@jhu.edu>
举例----------Examples----------
if (require(charmData) & require(BSgenome.Hsapiens.UCSC.hg18)) {
phenodataDir <- system.file("extdata", package="charmData")
pd <- read.delim(file.path(phenodataDir, "phenodata.txt"))
pd <- subset(pd, sampleID=="441_liver")
dataDir <- system.file("data", package="charmData")
setwd(dataDir)
# Read in raw data[读取原始数据]
rawData <- readCharm(files=pd$filename, sampleKey=pd)
# Find non-CpG control probes[非的CpG控制探针]
ctrlIdx <- getControlIndex(rawData, subject=Hsapiens)
# Get normalized methylation log-ratios[获取归甲基log比]
m <- methp(rawData, controlIndex=ctrlIdx, returnM=TRUE)
# Estimate percentage methylation[估计百分比甲基化]
ngc <- countGC(rawData)
p <- methPercent(m, ngc=ngc)
}
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|