我正在尝试对包含数百到数千列的数据表进行某些统计,例如 t 检验。数据的格式是我比较的两组值在同一列中。
所以,基本上我的第一次尝试是像下面这样剪切和粘贴;
NN <-read.delim("E:/output.txt")
View(NN)
attach(NN)
#output p-values of 100 t-tests
sink(file="E:/ttest.txt", append=TRUE, split=FALSE)
t.test(Tree1[1:13],Tree1[14:34])$p.value
t.test(Tree2[1:13],Tree2[14:34])$p.value
t.test(Tree3[1:13],Tree3[14:34])$p.value
……………………
随着数据的增长,这变得越来越不切实际。有没有办法通过每一列依次循环这些 t 检验并将输出保存到文件中?
提前致谢。