segments.p(DNAcopy)
segments.p()所属R语言包:DNAcopy
p-values for the change-points
P-值的变化点
译者:生物统计家园网 机器人LoveR
描述----------Description----------
This program computes pseudo p-values and confidence intervals for the change-points found by the circular binary segmentation (CBS) algorithm.
此方案计算伪p值和循环二元分割算法(CBS)的变点的置信区间。
用法----------Usage----------
segments.p(x, ngrid=100, tol=1e-6, alpha=0.05, search.range=100, nperm=1000)
参数----------Arguments----------
参数:x
an object of class DNAcopy
一个类DNAcopy的对象
参数:ngrid
number grid points for the integral
积分网格点
参数:tol
tolerance level for calculating nu
公差等级为计算NU
参数:alpha
Confidence level is 1-alpha
置信水平为1-α
参数:search.range
statistic is maximized over nu +/- search.range
统计数字是最大化了NU + / - search.range
参数:nperm
number of permutations for confidence interval
置信区间排列的数
Details
详情----------Details----------
The p-values are obtained by applying Siegmund's approximation for the maximal statistic from binary segmenting consecutive segments within a chromosome. This p-value is only to give the relative importance of the change-points as the CBS is different from the algorithm used here.
通过申请从二进制的最大分割染色体内连续段的统计西格蒙德的近似P-值。 P-值是只给相对重要性的变化点在CBS从这里所用的算法是不同的。
The confidence intervals are obtained by a permutation algorithm. The data are permuted to the left and right of the identified change-point and the location of the maximal binary segmentation statistic computed. The confidence interval is given by the quantiles of the permutation distribution of the locations.
置换算法得到的置信区间。数据置换,以确定变化点和最大的二元分割,计算统计的位置的左侧和右侧。置信区间的位置排列分布的分位数。
The statistical properties of this confidence interval is unknown. It is used to give an idea of the uncertainity on the location of the change-points as the CBS is different from the algorithm used here.
这个置信区间的统计特性是未知之数。 CBS的是从这里使用的算法不同,它是用来给一个变化点的位置上的不确定度的想法。
值----------Value----------
a data frame with ten columns. The maximal statistic from binary segmentation, the p-values and lower and upper alpha/2 confidence limits (as genomic positions) are added to the six columns from the segment command.
与十列的数据框。从二元分割,P-值和α/ 2置信区间(如基因组的位置)上下最大的统计添加到六个段命令列。
NOTE: THE p VALUES ARE APPROXIMATE TAIL PROBABILITIES. ANY VALUE GREATER THAN 0.1 CAN HAVE LARGE ERROR. p > 1 ARE REPLACED WITH 1.
注:P值是近似的尾概率。任何值大于0.1,可以有较大的误差。 P> 1被替换为1。
作者(S)----------Author(s)----------
Venkatraman E. Seshan
举例----------Examples----------
# test code on an easy data set[一个简单的数据集上的测试代码]
set.seed(25)
genomdat <- rnorm(500, sd=0.1) +
rep(c(-0.2,0.1,1,-0.5,0.2,-0.5,0.1,-0.2),c(137,87,17,49,29,52,87,42))
plot(genomdat)
chrom <- rep(1:2,c(290,210))
maploc <- c(1:290,1:210)
test1 <- segment(CNA(genomdat, chrom, maploc))
segments.p(test1)
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|