我有一个带有一组基因组坐标的数据框。我希望使用nearest.gene()
一次一个打印结果来找到这些坐标周围的基因。我一直在努力在循环中运行该函数:
apply(gene_lst, 1, function (x) nearest.gene(chr=gene_lst$Chr, pos=gene_lst$Pos))
1 2 3 4 5 6 7 8 9 10
"ACBD3" "ACBD3" "ACBD3" "ACBD3" "ACBD3" "ACBD3" "ACBD3" "ACBD3" "ACBD3" "ACBD3"
它正在覆盖接下来的九个坐标的第一个输出。有没有更好的方法来运行这个函数?
我想您可以尝试使用以下代码,您应该在其中传递x
给您的函数nearest.gene()
apply(gene_lst, 1, function (x) nearest.gene(chr=x["Chr"], pos=x["Pos"]))
本文收集自互联网,转载请注明来源。
如有侵权,请联系 [email protected] 删除。
我来说两句