f.pvalue(sva)
f.pvalue()所属R语言包:sva
Quickly calculate F-test p-values
快速计算F检验的p值
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Calculate p-values from a parametric F-test comparing the models mod and mod0 for each row of the data set dat.
计算从一个参数比较模型国防部和每一行的数据集DAT MOD0 F检验的p值。
用法----------Usage----------
f.pvalue(dat,mod,mod0)
参数----------Arguments----------
参数:dat
A m genes by n arrays matrix of expression data
由n阵列基因表达数据矩阵A M基因
参数:mod
A n by k model matrix corresponding to the primary model fit (see model.matrix)
K型矩阵对应的主模型拟合A N(见model.matrix)的
参数:mod0
A n by k0 model matrix corresponding to the null model to be compared to mod.
一个n到K0模型矩阵对应的空模型进行比较,以国防部。
Details
详情----------Details----------
The data for test i should be in the ith row of dat, if mod0 is null, the first column of mod is used as the null model.
我应该是在第i行的DAT,如果MOD0是空的测试数据,国防部第一列作为空模型。
值----------Value----------
参数:p
A vector of p-values one for each row of dat.
一个向量p值每个DAT行之一。
作者(S)----------Author(s)----------
Jeffrey T. Leek <a href="mailto:jleek@jhsph.edu">jleek@jhsph.edu</a>, John Storey <a href="mailto:jstorey@princeton.edu">jstorey@princeton.edu</a>
举例----------Examples----------
## Not run: [#无法运行:]
## Load data[#将数据]
library(bladderbatch)
data(bladderdata)
## Obtain phenotypic data[#获取表型数据。]
pheno = pData(bladderEset)
edata = exprs(bladderEset)
batch = pheno$batch
mod = model.matrix(~as.factor(cancer), data=pheno)
mod0 = model.matrix(~1, data=pheno)
#Calculate the p-values[计算p值]
p <- f.pvalue(edat,mod,mod0)
hist(p)
## End(Not run)[#结束(不运行)]
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|