plotFDR(BGmix)
plotFDR()所属R语言包:BGmix
Plot estimated FDR etc. for BGmix output.
图,估计BGmix输出FDR等。
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Given a threshold on the posterior probabilities, genes are declared as null or differentially expressed. For any given threshold, the FDR (false discovery rate) and FNR (false non-discovery rate) can be estimated using the posterior probabilities. This function plots these quantities twice, once versus the threshold probabilities, and once versus the number of declared positives.
鉴于上的后验概率的阈值,基因被宣布为空或差异表达。 FDR(错误发现率)和FNR(假非发现率)对于任何给定的阈值,可以使用后验概率估计。此功能图数量的两倍,一旦与阈值的概率,而一旦与声明阳性。
用法----------Usage----------
plotFDR(res, ylim = NULL, q.plotfnr = F, q.plotpcut = T, q.plotnpos = T, ...)
参数----------Arguments----------
参数:res
list object output from 'calcFDR'
列表对象输出从calcFDR
参数:ylim
optional argument specifying limit for y-axis
可选的参数指定限额为y轴
参数:q.plotfnr
Include FNR in plots?
包括在图的FNR?
参数:q.plotpcut
Include the plot of error rates v. threshold on posterior probabilites?
包括错误率诉后probabilites阈值的图吗?
参数:q.plotnpos
Include the plot of error rates v. no. positives.
包括错误率诉没有的图。阳性。
参数:...
arguments passed to 'plot'
参数传递的图“
值----------Value----------
No value is returned to R. Results from BGmix model are output to files.
没有返回值从BGmix模型R.的结果输出到文件。
作者(S)----------Author(s)----------
Alex Lewin
举例----------Examples----------
## Note this is a very short MCMC run![#注意:这是一个很短的MCMC运行!]
## For good analysis need proper burn-in period.[#为了获得良好的分析,需要适当的老化期。]
data(ybar,ss)
outdir <- BGmix(ybar, ss, c(8,8), nburn=0, niter=100, nthin=1)
params <- ccParams(outdir)
fdr <- calcFDR(params)
par(mfrow=c(1,2))
plotFDR(fdr)
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|