【发布时间】:2020-07-31 06:15:53
【问题描述】:
我想通过 chr1-chr22 标记相同的 22 个数据集。这是一个文件的样子。 als.sumstats.meta.chr10.txt.gz。我想解压缩这些文件并从列中提取所有数据 $2,$1,$3,$9 并在新列中合计 $14 和 $15。完成后我想压缩这些文件。我想返回一个附加表,其中包含来自上述每一列的所有数据以进行下游分析。这就是我在压缩/解压缩部分遇到的问题。
#!/bin/bash
FILES=/ALSsummaryGWAS/Summary_Statistics_GWAS_2016/als.sumstats.meta.chr{1..22}.txt.gz
for f in $FILES;
do
echo "$FILES"
echo "extracting columns 2,1,3,9"
awk '{print $2,$1,$3,$9, $14+$15}' > ALSGWAS.txt
done
【问题讨论】: