我正在尝试在两个数据集之间进行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值可能会给我带来错误消息
答案 0 :(得分:0)
如果dataset1
和dataset2
的列号相等,您可以编写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
,然后执行测试。