☰

gawk: creación dinámica de un archivo de texto utilizando campos de archivo de entrada

tag-icon • tag-icon text-processing gawk
gawk: creación dinámica de un archivo de texto utilizando campos de archivo de entrada

Tengo un archivo de texto que se parece a este:

UICEX_0001  UICEX_0001_T1.bam   UICEX_0001_C2.bam   chr1:16946335   chr19:9064309   chr8:10480278                                                                               
UICEX_0003  UICEX_0003_T1.bam   UICEX_0003_C2.bam   chr1:16974893

Estoy intentando utilizar esta información para crear algo como lo siguiente, combinando cadenas e ingresando información dinámicamente desde el archivo de texto. Yo quiero:

  • recorrer cada línea
  • asignar las primeras tres columnas a variables
  • imprimir algo de texto con esas variables
  • luego recorra la columna del cuarto extremo
  • imprimir algo específico.

Aquí hay un ejemplo de mi archivo de salida:

Archivo de salida

load UICEX_0001_T1.bam
load UICEX_0001_C2.bam
goto chr1:16946335
collapse
snapshot UICEX_0001_chr1:16946335.png
goto chr19:9064309
collapse
snapshot UICEX_0001_chr19:9064309.png
goto chr8:10480278
collapse
snapshot UICEX_0001_chr8:10480278.png

load UICEX_0003_T1.bam
load UICEX_0003_C2.bam
collapse
snapshot UICEX_0003_chr1:16974893.png

lo que he probado

Pensé que podría hacer esto anidando comandos gawk. Aquí hay uno que he probado:

SAMPLEFILE = "2016-10-13_mutation_table.txt"
gawk -F";" 'BEGIN{
gawk -F";" -v SAMPLE=$(cat $SAMPLEFILE | cut -d"\t" -f1) -v BAMT=$(cat $SAMPLEFILE | cut -d"\t" -f2) -v BAMN=$(cat $SAMPLEFILE | cut -d"\t" -f3);}{print "new \nload " $BAMN;}
{print "new \nload " $BAMT;}
{awk "{for(i=4; i<=NF-1; i++){ print "goto " $i ; print "collapse\nsnapshot " $SAMPLE"_"$i".png";} }" 2016-10-13_mutation_table.txt;

}END{print "exit \n"}'

Pero al ejecutar esto se obtienen errores como este:

gawk: cmd. line:2: gawk -F";" -v SAMPLE=$(cat $SAMPLEFILE | cut -d"\t" -f1) ...
gawk: cmd. line:2:                     ^ syntax error

Respuesta1

try this

$ awk '{printf("load %s\nload %s\n",$2,$3);for(i=4;i<=NF;i++){printf("goto %s\ncollapse\nsnapshot %s_%s.png\n",$i,$1,$i)}}' test.t
    load UICEX_0001_T1.bam
    load UICEX_0001_C2.bam
    goto chr1:16946335
    collapse
    snapshot UICEX_0001_chr1:16946335.png
    goto chr19:9064309
    collapse
    snapshot UICEX_0001_chr19:9064309.png
    goto chr8:10480278
    collapse
    snapshot UICEX_0001_chr8:10480278.png
    load UICEX_0003_T1.bam
    load UICEX_0003_C2.bam
    goto chr1:16974893
    collapse
    snapshot UICEX_0003_chr1:16974893.png

Respuesta2

Lo que tuparecerpedir se puede hacer como

awk '
  NR > 1 {print ""}
  {
    printf("load %s\nload %s\n", $2, $3);
    for (i=4; i<=NF; i++) {
      if (NF > 4) printf("goto %s\n", $i);
      printf("collapse\nsnapshot %s_%s.png\n", $1, $i);
    }
  }
' samplefile

información relacionada