y原创:hxj7
关键词:awk
awk是生信人必须要掌握的命令行工具。为啥?因为它太强大了。我们举一个例子来说明。
假设我们有一个1000万行的文件,大概长这样:

怎么求第四列的平均数呢?
Python版本
我们可以用几行Python代码解决,比如这样:

其耗时:

R版本
用R来做计算也是很适合的,比如像这样:

其耗时:

awk版本
awk用一行代码就可以解决问题,像这样(注意耗时):

至此,我们可以看出,awk代码简单,但是性能却不差!在同样的机器上处理同样的文件,awk的运行时间是Python的一半左右,是R的大概十分之一。可以说,awk已经非常快了!
C版本
都说C快,让我们看看到底有多快。代码如下:
其耗时:

可以看出,C的版本也仅比awk的稍快一点点。但是,C的代码复杂多了!由此,我们可以粗略比较出awk是一个非常完美的文本处理工具!
如果有任何问题,欢迎交流!
(公众号:生信了)

网友评论