Início / unix / Perguntas / 421316
Accepted
user977828
Asked: 2018-02-02 15:04:27 +0800 CST

Concatenação de strings não funciona

  • 772
DASqv -v -H$H -c$cov $db $i | grep Recommend - | sed "s|Recommend ||g" - | sed "s|'||g" -` by itself produces `DAStrim -g20 -b25

Meu objetivo é combinar os resultados anteriores awk '{print $1 " " $2 " "$3 " $db $i"}'e canalizar todo o comando para um arquivo de saída > $(basename $i .las).DAStrim.

Infelizmente, só obtenho como resultado bananaDB ./bananaDB.100.lase não DAStrim -g20 -b25 bananaDB ./bananaDB.100.lascom o seguinte código:

#!/bin/bash

db=bananaDB
H=6973
cov=38

for i in $(find . -type f -name "*.*.las");
do
  #cat <<EOF
  qsub <<EOF

#!/bin/bash -l

#PBS -N DASqv
#PBS -l walltime=48:00:00
#PBS -j oe
#PBS -l mem=1G
#PBS -l ncpus=1
#PBS -M [email protected]
##PBS -m bea

cd \$PBS_O_WORKDIR

source activate thegenemyers


DASqv -v -H$H -c$cov $db $i | grep Recommend - | sed "s|Recommend ||g" - | sed "s|'||g" - | awk '{print $1 " " $2 " "$3 " $db $i"}' > $(basename $i .las).DAStrim

EOF

done

ATUALIZAR 

DASqv -v -H$H -c$cov $db $i

produzido:

DASqv -c38 bananaDB ./bananaDB.100.las

Input:   16,450reads,   210,758,575 bases (another 9,934 were < H-length)

Histogram of q-values (average 10 best)

                 Input                 QV

    50:    1494189    0.2%       380302   18.0%

    49:     364713    0.0%          484    0.0%
    48:     545846    0.1%          423    0.1%
    47:     650479    0.2%          466    0.1%
    46:     835282    0.3%          548    0.1%
    45:    1054589    0.4%          648    0.1%
    44:    1299423    0.5%          775    0.2%
    43:    1644281    0.7%          895    0.2%
    42:    2036915    0.9%         1193    0.3%
    41:    2571126    1.2%         1334    0.4%
    40:    3518594    1.5%         1647    0.5%
    39:    3641660    1.9%         2046    0.6%
    38:    5026473    2.4%         2291    0.7%
    37:    6243982    3.1%         2708    0.9%
    36:    7600704    3.9%         3301    1.1%
    35:    9313754    4.9%         4002    1.3%
    34:   11257936    6.0%         4676    1.6%
    33:   13508338    7.5%         5544    1.9%
    32:   15981847    9.1%         6552    2.3%
    31:   18648809   11.1%         7771    2.7%
    30:   22290239   13.4%         9124    3.3%
    29:   25083448   16.0%        10624    3.9%
    28:   29566164   19.1%        12874    4.6%
    27:   33339712   22.6%        15482    5.5%
    26:   37891335   26.6%        18869    6.6%
    25:   44146531   31.2%        23307    7.9%
    24:   44948068   35.9%        28142    9.5%
    23:   50951224   41.3%        33590   11.5%
    22:   55009718   47.1%        42157   13.9%
    21:   57456151   53.1%        52181   16.9%
    20:   60635065   59.4%        63207   20.6%
    19:   58423422   65.6%        76426   25.0%
    18:   58472922   71.7%        91565   30.2%
    17:   55127848   77.5%       107289   36.4%
    16:   50395382   82.7%       123758   43.6%
    15:   43893354   87.3%       136465   51.4%
    14:   36509552   91.2%       145632   59.8%
    13:   28654550   94.2%       145540   68.2%
    12:   21245809   96.4%       138232   76.2%
    11:   14560980   97.9%       121403   83.2%
    10:    9345155   98.9%        98071   88.8%
     9:    5395169   99.5%        73996   93.1%
     8:    2894210   99.8%        52246   96.1%
     7:    1335673   99.9%        33845   98.0%
     6:     581470  100.0%        19476   99.2%
     5:     201756  100.0%         9367   99.7%
     4:      76322  100.0%         3760   99.9%
     3:      18979  100.0%         1082  100.0%
     2:       4751  100.0%          264  100.0%
     1:        456  100.0%           41  100.0%
     0:       2686  100.0%           38  100.0%

  Recommend 'DAStrim -g20 -b25'

O que eu perdi?

Agradeço antecipadamente.

bash awk

1 respostas

  1. Best Answer

    Você está tornando as coisas mais difíceis do que precisam ser e encontrando espaços em branco e problemas de citação. Tente algo como o seguinte:

    Etapa 1: crie um script autônomo que faça o que você deseja com um ou mais de seus arquivos de dados, dados os argumentos e nomes de arquivo apropriados na linha de comando.

    #!/bin/sh

# use the first 3 arguments for the values to pass to DASqv
db="$1"
H="$2"
cov="$3"

# use shift to get rid of them once we have them in variables, ...
shift 3

# ... so we can loop over the remaining filenames (1 or more) on the command line
for filename in "$@" ; do
  outfile="$(basename "$filename" .las).DAStrim"
  qsub <<EOF
#!/bin/bash -l

#PBS -N DASqv
#PBS -l walltime=48:00:00
#PBS -j oe
#PBS -l mem=30G
#PBS -l ncpus=1
#PBS -M [email protected]
##PBS -m bea

cd "\$PBS_O_WORKDIR"

source activate thegenemyers

DASqv -v -H"$H" -c"$cov" "$db" "$filename" | 
  sed -n -e '/Recommend/ {
               s/Recommend //;
               s/\x27//g;
               s:$: "$db" "$filename":;
               p
             }' > "$outfile"

EOF

done

    (o sedscript no meio disso pode estar todo em uma linha, mas as alimentações de linha extras e o recuo tornam-no mais legível sem alterar o que ele faz / como é executado. Além disso, observe o uso de \x27para remover todos os single- caracteres de aspas. 0x27é a notação hexadecimal para o caractere de aspas simples ASCII)

    salve-o como, por exemplo, submit-jobs.she torne-o executável com chmod +x submit-job.sh.

    Passo 2: Teste

    Teste se o script faz o que você deseja usando-o manualmente para enviar tarefas. por exemplo, executar:

    /path/to/submit-jobs.sh bananaDB 6973 38 /path/to/somefile.las

    Modifique o script, se necessário, até que ele faça exatamente o que você deseja.

    Etapa 3: agora use findpara enviar vários trabalhos usando o script:

    find . -type f -name '*.las' -exec /path/to/submit-jobs.sh bananaDB 6973 38 {} +

    Etapa 4: (opcional) transforme a etapa 3 em um script que você pode executar com diferentes argumentos para evitar que você tenha que digitar o find ...comando toda vez que quiser outra execução com valores ligeiramente diferentes. por exemplo

    #!/bin/sh
find . -type f -name '*.las' -exec /path/to/submit-jobs.sh "$1" "$2" "$3"

    Se você salvou isso como find-and-submit.she o tornou executável com chmod +x, você o executaria como:

    find-and-submit.sh bananaDB 6973 38

    Esse script da etapa 4 pode até ter um loop for para as variáveis ​​para que, por exemplo, ele envie trabalhos para $covvalores de 35 a 45, em vez de exigir $covque seja um dos argumentos.

    • 2

relate perguntas