我有一个这样的数据集:
Row B1 B2 B3 G1 G2 G3
1 6.19 0.69 0.30 0 0 0
2 100 100 100 94.87 98.73 90.56
3 100 90.24 100 40 78.94 87.5
4 93.90 96.15 95.29 98.82 100 100
5 100 98.38 100 88.88 95.87 83.56
6 93.90 96.15 95.29 98.82 100 100
依此类推,对于 25000 行。对于每一行,我需要计算 B 列和 G 列之间的差异,并得到置信区间。这是我到目前为止所拥有的:
data=read.table("1.txt", header=T, sep="\t")
boys=data[,grep('B',names(data))]
girls=data[,grep('G',names(data))]
ci= sapply(1:nrow(data),function(i){wilcox.test(as.numeric(boys[i,]),as.numeric(girls[i,]),exact=FALSE,correct=FALSE,conf.int=TRUE,conf.level=0.95)$conf.int})
我的问题是,每次我这样做时,我都会得到置信区间以显示区间的两个数字,一个在另一个之上。例如,像这样:5 15(例如,对于第 1 行)
如何从每一行获取 CI 以显示在一行中,如下所示:5,15
将不胜感激在这方面的任何帮助!提前致谢!