someMTP-package(someMTP)
someMTP-package()所属R语言包:someMTP
Some Multiple Testing Procedures
一些多次测试程序
译者:生物统计家园网 机器人LoveR
描述----------Description----------
It is a collection of functions for Multiplicty Correction and Multiple Testing.
这是一个集合多重数校正和多种测试功能。
Details
详细信息----------Details----------
</table>
</ TABLE>
(作者)----------Author(s)----------
livio finos
Maintainer: <livio@stat.unipd.it>
参考文献----------References----------
For weighted methods:
Benjamini, Hochberg (1997). Multiple hypotheses testing with weights. Scand. J. Statist. 24, 407-418.
Finos, Salmaso (2007). FDR- and FWE-controlling methods using data-driven weights. Journal of Statistical Planning and Inference, 137,12, 3859-3870.
For LSD test:
J. Lauter, E. Glimm and S. Kropf (1998). Multivariate test based on Left-Spherically Distributed Linear Scores. The Annals of Statistics, Vol. 26, No. 5, 1972-1988
L. Finos (2011). A note on Left-Spherically Distributed Test with covariates, Statistics and Probabilty Letters, Volume 81, Issue 6, June 2011, Pages 639-641
实例----------Examples----------
set.seed(13)
y <- matrix(rnorm(5000),5,1000) #create toy data[创建玩具数据]
y[,1:100] <- y[,1:100]+3 #create toy data[创建玩具数据]
p <- apply(y,2,function(y) t.test(y)$p.value) #compute p-values[计算p-值]
M2 <- apply(y^2,2,mean) #compute ordering criterion[计算排序标准]
fdr <- p.adjust(p,method="BH") #(unweighted) procedure, fdr control[(未加权)程序,FDR控制]
sum(fdr<.05)
fdr.w <- p.adjust.w(p,method="BH",w=M2) #weighted procedure, weighted fdr control[加权程序,加权FDR控制]
sum(fdr.w<.05)
fwer <- p.adjust(p,method="holm") #(unweighted) procedure, fwer control[(未加权)程序,fwer控制]
sum(fwer<.05)
fwer.w <- p.adjust.w(p,method="BHfwe",w=M2) #weighted procedure, weighted fwer (=fwer) control[加权程序,的加权fwer(= fwer)控制]
sum(fwer.w<.05)
plot(M2,-log10(p))
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|