我有这样的数据集:
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})
我的问题是,每次我这样做时,我都会获得置信区间,并显示间隔的2个数字。例如,像这样: 五 15 (对于第1行,比方说)
如何让每行的CI显示在一行中,如下所示:5,15
非常感谢这方面的任何帮助! 提前谢谢!