mt.sample.teststat(multtest)
mt.sample.teststat()所属R语言包:multtest
Permutation distribution of test statistics and raw (unadjusted) p-values
置换分布的检验统计量和原料(未经调整)的P-值
译者:生物统计家园网 机器人LoveR
描述----------Description----------
These functions provide tools to investigate the permutation distribution of test statistics, raw (unadjusted) p-values, and class labels.
这些功能提供的工具来调查测试统计,原料(未经调整)p的值和类标签的排列分布。
用法----------Usage----------
mt.sample.teststat(V,classlabel,test="t",fixed.seed.sampling="y",B=10000,na=.mt.naNUM,nonpara="n")
mt.sample.rawp(V,classlabel,test="t",side="abs",fixed.seed.sampling="y",B=10000,na=.mt.naNUM,nonpara="n")
mt.sample.label(classlabel,test="t",fixed.seed.sampling="y",B=10000)
参数----------Arguments----------
参数:V
A numeric vector containing the data for one of the variables (genes).
一个数字变量(基因)的数据向量。
参数:classlabel
A vector of integers corresponding to observation (column) class labels. For k classes, the labels must be integers between 0 and k-1. For the blockf test option, observations may be divided into n/k blocks of k observations each. The observations are ordered by block, and within each block, they are labeled using the integers 0 to k-1.
观察(列)类的标签对应的整数向量。对于k类,标签必须是0k-1之间的整数。 blockf测试选项,观测n/kk观察每块可分为。的意见,下令块,每块内,他们都使用整数0k-1标记。
参数:test
A character string specifying the statistic to be used to test the null hypothesis of no association between the variables and the class labels.<br> If test="t", the tests are based on two-sample Welch t-statistics (unequal variances). <br> If test="t.equalvar", the tests are based on two-sample t-statistics with equal variance for the two samples. The square of the t-statistic is equal to an F-statistic for k=2. <br> If test="wilcoxon", the tests are based on standardized rank sum Wilcoxon statistics.<br> If test="f", the tests are based on F-statistics.<br> If test="pairt", the tests are based on paired t-statistics. The square of the paired t-statistic is equal to a block F-statistic for k=2. <br> If test="blockf", the tests are based on F-statistics which adjust for block differences (cf. two-way analysis of variance).
如果一个字符串指定的统计,被用来测试空假设之间没有关联的变量和类的标签。参考test="t",测试是基于两样本韦尔奇t-统计(不平等的差异)。 <br>如果test="t.equalvar",测试是基于两样本t-统计量等于两个样本方差。 t-统计量的平方等于k=2F-统计。参考如果test="wilcoxon",测试标准化秩秩统计基础。参考如果test="f",测试F-统计基础。参考如果test="pairt"测试配对t-统计的基础上。配对的t-统计的平方等于块k=2F-统计。 <br>如果test="blockf",测试是根据F-统计调整块的差异(参见双向方差分析)。
参数:side
A character string specifying the type of rejection region.<br> If side="abs", two-tailed tests, the null hypothesis is rejected for large absolute values of the test statistic.<br> If side="upper", one-tailed tests, the null hypothesis is rejected for large values of the test statistic.<br> If side="lower", one-tailed tests, the null hypothesis is rejected for small values of the test statistic.
如果一个字符串指定排斥区域类型。参考如果side="abs",一个side="upper",双尾检验,检验统计量的大绝对值拒绝零假设。参考尾检验,零假设被拒绝的检验统计量的大值。如果side="lower",单尾测试参考,拒绝零假设检验统计量的小值。
参数:fixed.seed.sampling
If fixed.seed.sampling="y", a fixed seed sampling procedure is used, which may double the computing time, but will not use extra memory to store the permutations. If fixed.seed.sampling="n", permutations will be stored in memory. For the blockf test, the option n was not implemented as it requires too much memory.
fixed.seed.sampling="y"如果,一个固定的种子抽样程序使用,这可能会增加一倍的计算时间,但不会使用额外的内存来存储的排列。如果fixed.seed.sampling="n",排列将存储在内存中。 blockf测试,选项n没有实施,因为它需要太多的内存。
参数:B
The number of permutations. For a complete enumeration, B should be 0 (zero) or any number not less than the total number of permutations.
排列数。对于一个完整的枚举,B应该是0(零)或任何数量不超过总数的排列。
参数:na
Code for missing values (the default is .mt.naNUM=--93074815.62). Entries with missing values will be ignored in the computation, i.e., test statistics will be based on a smaller sample size. This feature has not yet fully implemented.
缺失值的代码(默认为.mt.naNUM=--93074815.62)。与缺失值的参赛作品将在计算中忽略,即检验统计量将基于一个小样本大小。此功能尚未完全落实。
参数:nonpara
If nonpara="y", nonparametric test statistics are computed based on ranked data. <br> If nonpara="n", the original data are used.
如果nonpara“=”Y“,非参数检验统计排名数据计算。如果nonpara=“N”的参考,使用原始数据。
值----------Value----------
For mt.sample.teststat, a vector containing B permutation test statistics. <br> <br> For mt.sample.rawp, a vector containing B permutation unadjusted p-values. <br><br> For mt.sample.label, a matrix containing B sets of permuted class labels. Each row corresponds to one permutation.
mt.sample.teststat,向量B置换检验统计量。参考参考mt.sample.rawp,B置换未经调整的p值向量。参考参考mt.sample.label,矩阵包含B套置换类标签。每一行对应一个排列。
作者(S)----------Author(s)----------
Yongchao Ge, <a href="mailto:yongchao.ge@mssm.edu">yongchao.ge@mssm.edu</a>, <br>
Sandrine Dudoit, <a href="http://www.stat.berkeley.edu/~sandrine">http://www.stat.berkeley.edu/~sandrine</a>.
参见----------See Also----------
mt.maxT, mt.minP, golub.
mt.maxT,mt.minP,golub。
举例----------Examples----------
# Gene expression data from Golub et al. (1999)[Golub等基因表达数据。 (1999)]
data(golub)
mt.sample.label(golub.cl,B=10)
permt<-mt.sample.teststat(golub[1,],golub.cl,B=1000)
qqnorm(permt)
qqline(permt)
permt<-mt.sample.teststat(golub[50,],golub.cl,B=1000)
qqnorm(permt)
qqline(permt)
permp<-mt.sample.rawp(golub[1,],golub.cl,B=1000)
hist(permp)
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|