- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我想找到这样的模式:当前行中的第 2 列是“C”,下一行中的第 2 列是“G”。文件的第 4 列是“CG”。我想比较第 1 与第 2、第 3 与第 4、第 5 与第 6,依此类推。然后打印几行当前行和下一行。“C”可以出现在偶数行和奇数行中。
输入如下:
chr1 C 10467 CHH CT 0.0 0 1
chr1 C 10469 CG CG 0.0 0 1
chr1 G 10470 CG CG 0.0 0 8
chr1 C 10471 CG CG 0.0 0 1
chr1 G 10472 CG CG 1.0 8 8
预期输出是,用制表符分隔符分隔:
chr1 C 10469 CG CG 0.0 0 1
chr1 G 10470 CG CG 0.0 0 8
chr1 C 10471 CG CG 0.0 0 1
chr1 G 10472 CG CG 1.0 8 8
我的代码是:
awk '{a=$2; c=$4; d=$0; e=NR; getline; f=$2; g=$4} {if (a == "C" && f == "G" && c == "CG" && g == "CG") {print d,e,"\n",$0,NR}}' input_file
我使用 getline 并检查下一行是否有“G”。问题是,如果我这样做,awk 将直接转到第三行,并且会错过一些行。例如,输入的第 2 列是:
Line 1: G
Line 2: C
Line 3: G
Line 4: C
预期的输出是第2行和第3行。但是,awk直接从第一行转到第三行,而不是逐行。因此,输出为无。
亲切的问候!
最佳答案
编辑(要将每一行与其下一行进行比较,请使用此行):现在使用 OP 的新示例添加此解决方案。
awk '
FNR>1{
if(secCol=="C" && $2=="G" && fourthCol=="CG" && $4=="CG"){
print prevLine ORS $0
}
}
{
secCol=$2
fourthCol=$4
prevLine=$0
}
' Input_file
说明:为上述内容添加详细说明。
awk '
##Starting awk program from here.
FNR>1{
##Checking condition if current line number is more than 1 then do following.
if(secCol=="C" && $2=="G" && fourthCol=="CG" && $4=="CG"){
##Checking condition if secCol is C AND 2nd column is G AND fourthCol is CG and 4th column is CG then do following.
print prevLine ORS $0
##Printing prevLine ORS and current line.
}
}
{
secCol=$2
##Creating secCol with 2nd column of current line.
fourthCol=$4
##Creating fourthCol with 4th column of current line.
prevLine=$0
##Setting prevLine to current line value.
}
' Input_file ##Mentioning Input_file name here.
初始解决方案(比较每个奇数行和偶数行):(OP 的示例在编辑后变得更加清晰,但也将这个解决方案保留在这里,以供将来的读者使用,以防它有帮助)您能否尝试以下操作,仅按照所示示例编写。这会检查上一行是否有第四列(fourthCol)也是 CG
,以防万一您不需要它,然后从以下内容中删除 && foruthCol=="CG"
.
awk '
FNR%2==0{
if(secCol=="C" && $2=="G" && fourthCol=="CG" && $4=="CG"){
print prevLine ORS $0
}
prevLine=secCol=fourthCol=""
next
}
{
secCol=$2
fourthCol=$4
prevLine=$0
}
' Input_file
输出如下。
chr1 C 10469 CG CG 0.0 0 1
chr1 G 10470 CG CG 0.0 0 8
chr1 C 10471 CG CG 0.0 0 1
chr1 G 10472 CG CG 1.0 8 8
说明:为上述内容添加详细说明。
awk ' ##Starting awk program from here.
FNR%2==0{ ##Checking condition if line number is divided by 2 or not.
if(secCol=="C" && $2=="G" && fourthCol=="CG" && $4=="CG"){
##Checking condition if secCol is C AND 2nd column is G AND fourthCol is CG and 4th column is CG then do following.
print prevLine ORS $0 ##Printing prevLine ORS and current line.
}
prevLine=secCol=fourthCol="" ##Nullifying prevLone, secCol, fourthCol here.
next ##next will skip all further statements from here.
}
{
secCol=$2 ##Creating secCol with 2nd column of current line.
fourthCol=$4 ##Creating fourthCol with 4th column of current line.
prevLine=$0 ##Setting prevLine to current line value.
}
' Input_file ##Mentioning Input_file name here.
关于awk - 比较 awk 中的当前行和下一行,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/65360525/
我试图通过这段代码读取未知数量的整数: while (1) { int c = getchar (); if (c == EOF) break;
我正试图找到一个类似于谷歌分析日期选择器的日期选择器: 知道 jQuery 是否提供了类似的东西吗? 最佳答案 这个 Twitter Bootstrap 风格的日期范围选择器非常接近。 https:/
我正在使用 javascript。如何获取当前 URL 的路径并将其分配给我的代码?这是我的代码: $(document).ready(function() { $(".share").hides
如何获得今天的Julian day number (JDN)相等的?或任何日期? 我看了又看,但只发现了一些产生“year-dayOfYear”的函数,而不是:2457854。 最佳答案 在 bash
我有相当简单的 UDP 服务器写在 c 上。 有时我需要知道在套接字中排队的所有 udp 数据包(字节)的当前长度。 据我了解,getsockopt 没有得到这样的信息。 欢迎使用 Linux 和 F
我一直在寻找几个小时来找到一个可以在图像中添加诸如“填充:5px”之类的东西的插件。每个人都通过纯 html 做到这一点吗?我们的客户需要一种方法来简单地使用按钮或右键单击上下文菜单来添加它。有什么建
是否有可能获得当前正在执行的 TCL 脚本的完整路径? 在 PHP 中,它将是:__FILE__ 最佳答案 根据“当前正在执行的 TCL 脚本”的含义,您实际上可能会寻找 info script ,甚
我最近从直接使用 ISession 转向了包装的 ISession,即工作单元类型模式。 我曾经使用 SQL Lite(内存中)对此进行测试。我有一个简单的帮助器类,它配置我的 SessionFact
我按照步骤操作 here在 WebStorm 中配置代码完成和其他内容,但我仍然收到以下语法错误。 我该如何解决这个问题? 最佳答案 通过相应地将“JavaScript 语言版本”(Settings/
我可以为我团队的 TFS 当前 Sprint 任务板添加书签吗?我们有两周的冲刺,因此 URL 每两周更改一次。 默认 URL 的形式为: http://[Server]/tfs/[Project]/
是否有 Subversion 命令可以显示当前版本号? 在svn checkout之后,我想启动一个脚本并需要变量中的修订号。如果有像 svn info get_revision_number 这样的
我正在编写表单的一个组件 首次安装组件时,sources={{}} ,一本空字典。由于该组件包装了现有的 Javascript 库,因此我正在实现一个自定义比较函数。为了让这个 diffing 函数
无论系统时间设置为多少以及机器所在的时区,我都需要正确的 UTC 时间。 (即使我必须打电话到互联网才能同步......) 是否有一些库或其他方法可以优雅地做到这一点? 最佳答案 如果您想获得准确可靠
我一边编码,一边拿出一些我和 friend 建立的旧网站来重新开始工作。我已经有一段时间没有做过任何 AJAX 了,当我试图找出我的代码失败的地方时,我发现没有显示很多资源。我猜这是因为我使用的是旧方
由于对性能的巨大影响,我从不怀疑我现在的桌面CPU是否有分支预测。当然可以。但各种 ARM 产品又如何呢? iPhone或Android手机有分支预测吗?较旧的任天堂 DS?基于 PowerPC 的
我有一个具有以下有效负载的 JWT: { "id": "394a71988caa6cc30601e43f5b6569d52cd7f6df", "jti": "394a71988caa6cc30
从其他一些帖子中,我能够通过以下方式获取当前 URI: 但是以下方法不起作用: 我很好奇为什么上面的方法不起作用,以及如何将当前 URI 分配给字符串。 最佳答案 每the javadocs ,g
我在表格 View 中有几个单元格。现在在任何给定的时间点,我想计算 View 中单元格的当前高度,即如果它是 View 的 3/4,它应该返回 (cellheight)*3/4 高度。 我通过以下方
这是网站的身份验证脚本。这安全吗?是最近的节目吗?它已经过时了吗?是否有“更好更安全的方法”我很新,但我没有看到太多地方使用 header 授权。 如有任何帮助,我们将不胜感激!这是我制作的第一个登录
我已经在其他 stackoverflow 线程上检查过这个错误,但在我的代码中没有发现任何错误。也许我累了,但我觉得还好。 网站.urls.py: from django.conf.urls impo
我是一名优秀的程序员,十分优秀!