Explorar el Código

extract_multiqc

master
LUYAO REN hace 4 años
padre
commit
dd056011c3
Se han modificado 2 ficheros con 2 adiciones y 1 borrados
  1. +1
    -0
      codescripts/extract_multiqc.py
  2. +1
    -1
      tasks/D5_D6.wdl

+ 1
- 0
codescripts/extract_multiqc.py Ver fichero

dat = dat.loc[:, dat.columns.str.endswith('ALL')] dat = dat.loc[:, dat.columns.str.endswith('ALL')]
dat_transposed = dat.T dat_transposed = dat.T
benchmark = dat_transposed.loc[:,['sample_id','METRIC.Precision','METRIC.Recall']] benchmark = dat_transposed.loc[:,['sample_id','METRIC.Precision','METRIC.Recall']]
benchmark['sample_id'] = benchmark.index
benchmark.columns = ['Sample','Precision','Recall'] benchmark.columns = ['Sample','Precision','Recall']


#output #output

+ 1
- 1
tasks/D5_D6.wdl Ver fichero

sample_j=$(echo $j | cut -f7 -d_) sample_j=$(echo $j | cut -f7 -d_)
rep_i=$(echo $i | cut -f8 -d_) rep_i=$(echo $i | cut -f8 -d_)
rep_j=$(echo $j | cut -f8 -d_) rep_j=$(echo $j | cut -f8 -d_)
if [[ $sample_i == "LCL5" ] && [ $sample_j == "LCL6" ] && [ $rep_i == $rep_j ]];then
if [ $sample_i == "LCL5" ] && [ $sample_j == "LCL6" ] && [ $rep_i == $rep_j ];then
cat $i | grep -v '##' | cut -f1,2,4,5,10 | cut -d ':' -f1 > LCL5.txt cat $i | grep -v '##' | cut -f1,2,4,5,10 | cut -d ':' -f1 > LCL5.txt
cat $j | grep -v '##' | cut -f10 | cut -d ':' -f1 > LCL6.txt cat $j | grep -v '##' | cut -f10 | cut -d ':' -f1 > LCL6.txt
paste LCL5.txt LCL6.txt > sister.txt paste LCL5.txt LCL6.txt > sister.txt

Cargando…
Cancelar
Guardar