我将df拆分为包含500个组的列表:
c1=cut(SNP_Allele_Frequency$College_SE,500)
splitc1=split(SNP_Allele_Frequency,c1,drop=FALSE)
然后我得到了列表中所有组件的均值:
sink("acb1.txt")
lapply(splitc1, function(x) mean(x$ACB))
sink()
问题是,我需要将输出保存为Excel可读或R可用格式,但输出包含所有信息,我必须手动复制每个案例才能使用它。我的意思是我只需要平均值值(例如0.3879208)而不是输出的其余部分(例如$ (-0.001,0.01]
[1])我想将它保存为列/向量,以便我可以直接使用它。
$`(-0.001,0.01]`
[1] 0.3879208
$`(0.01,0.02]`
[1] 0.3855017
$`(0.02,0.03]`
[1] 0.379335
$`(0.03,0.04]`
[1] 0.3769094
答案 0 :(得分:0)
我已经解决了这个问题:ibs1=as.character(unlist(lapply(splitc1, function(x) mean(x$IBS))))
write.table(ibs1,file="ibs1.csv",row.names=FALSE)