R中的T.test ...获得多个p值

时间:2013-10-06 01:39:03

标签: r p-value

我正在尝试在两个数据集之间进行T.test ..

这是'数据集1'

TIME        5    10      15      20      25
Specimen no.                        
1          15.2 30.5     41     12.5         16.2
2          13.1 16.2    12.5    Na       13.2
3          16.11 45.7   11.4    18.9    11.7
4          11.2  Na    9.11    20.7      19

另一个'数据集2'

TIME        5   10        15      20    25

Specimen no.                        
1       11.8    34.8    14.2    19.9    23.4
2         NA    6.4     29.2    32.7    17.1
3       10.0    35.5    38.5    28.3    27.3
4       18.7    NA      11.5    14.6    18.9

我只想使用t.test比较从dataset1到dataset2的每个5秒间隔。我如何获得p值列表

到目前为止我只有这个

t.test(dataset1[[2]],dataset2[[2]])$p.value...

显然这是错误的......因为我不知道如何从dataset1和dataset2访问某些列...我也有NA值可能会给我带来错误消息

1 个答案:

答案 0 :(得分:0)

如果dataset1dataset2的列号相等,您可以编写for循环:

pval = rep(0, ncol(dataset1))
for(i in 1:ncol(dataset1)){
    pval[i] = t.test(dataset1[, i], dataset2[, i], na.action = na.omit)$p.value
}

此处,dataset1[, i]表示i列,na.omit将删除NA,然后执行测试。