nullDistDoublyTestedEdges(ppiStats)
nullDistDoublyTestedEdges()所属R语言包:ppiStats
Null distribution of number of reciprocated,
空分布投桃报李数,
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Calculate the null distribution of the number of reciprocated, unreciprocated and missing edges in a stochastic model where each edge is tested twice.
空分布计算,其中每条边两次测试中的随机模型投桃报李,unreciprocated的和丢失的边缘。
用法----------Usage----------
nullDistDoublyTestedEdges(deltaMax, n, pFP, pFN)
参数----------Arguments----------
参数:deltaMax
Integer. Distributions will be calculated for model parameter delta=0, 1, 2, ..., deltaMax.
整数。将计算分布模型参数Delta=0, 1, 2, ..., deltaMax。
参数:n
Integer. The parameter n of the model.
整数。该模型的参数n。
参数:pFP
Numeric. The parameter pFP of the model.
数字。该模型的参数亲民党。
参数:pFN
Numeric. The parameter pFN of the model.
数字。该模型的参数PFN。
Details
详情----------Details----------
The model is described in the vignette Stochastic and systematic errors in PPI data, by looking at unreciprocated in- or out-edges by W. Huber, T. Chiang and R. Gentleman.
模型中描述的小品文随机和PPI数据的系统误差,通过寻找在unreciprocated,或由T.江和河君子·胡贝尔,出边。
This function can be quite slow, its runtime grows quickly with deltaMax (and is roughly independent of n, pFP, pFN). The example below should take only a few seconds on a reasonable computer, though.
此功能可以是相当缓慢,其运行快速增长deltaMax(大约是独立的n,pFP,pFN)。下面的例子中应采取合理的计算机上,虽然只有几秒钟。
值----------Value----------
3d array with dimensions nMax+1 x nMax+1 x deltaMax+1 whose element p[nr+1, nu+1, delta+1] is the corresponding joint probability. nMax+1 is calculated (probably too conservatively) by the function to make sure that no probability leaks out of the array.
同尺寸的3D阵列nMax+1XnMax+1XdeltaMax+1元素p[nr+1, nu+1, delta+1]是相应的联合概率。 nMax+1计算功能(可能过于保守),以确保无的概率泄漏的阵列。
作者(S)----------Author(s)----------
Wolfgang Huber <a href="http://www.ebi.ac.uk/huber">http://www.ebi.ac.uk/huber</a>
举例----------Examples----------
p = nullDistDoublyTestedEdges(32, 1000, pFP=0.001, pFN=0.15)
if(interactive() && require("RColorBrewer"))
for(k in 1:dim(p)[3]) {
image(sqrt(p[,,k]), xlab=expression(N[rec]), ylab=expression(N[unrec]),
main = expression(P(N[rec], N[unrec]~";"~ delta^"*", n, p[FP], p[FN])),
x = 1:dim(p)[1], y = 1:dim(p)[2],
col = colorRampPalette(brewer.pal(9, "GnBu"))(256))
text(35, 35, paste("delta", k, sep="="))
}
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|