晚上好,我正在做一个类项目,我尝试做多个未配对的2个样本t-测试,然后存储它们的p-值,以便以后只处理下面的p值,这是我一直在尝试的代码:
pVals_1Beta <-vector("numeric", length = nrow(group1_Y_Beta))
for (i in 1:nrow(group1_Y_Beta)) {
pVals_1Beta[i] <- t.test(x = group1_Y_Beta$values[i,],
y = group1_N_Beta$values[i,],
paired = FALSE,
var.equal =FALSE,
conf.level = 0.95)$p.value
}其中group1_Y_Beta和group1_N_Beta有两列(values和ind)和大约110312行,我想要运行未配对的t检验,比较两个组的值,并存储所有110312个p值。当我尝试运行这个时,我得到:
group1_Y_Beta$valuesi中的错误:不正确的维数
任何帮助,如何调整我的代码,使其工作,将不胜感激。
THanks,LIz
发布于 2018-11-23 03:51:50
由于group1_N_Beta和group1_Y_Beta是2D对象,因此需要(1)行和(2)列标识符才能获得特定单元格的值。但是,由于您已经使用$符号指定了列的名称,所以只需要提供一个数字(或数字向量)就可以完成查询。将[i,] ("ith行,所有列“)替换为[i]
https://stackoverflow.com/questions/53439897
复制相似问题