优化阅读性和性能 if $length -eq 1 ; then grep -v ## $vcf awk -v patient=$patient BEGINOFS=tNR==1fori=1;i=NF;i++a$i=iNR=1print $a#CHROM$aPOS$aID$aREF$aALT$aQUAL$aFILTER$aINFO$aFORMAT$apatient $dirdatevcf_no
如果$length等于1,将从$vcf中提取出指定病人的信息,并将结果存储在$dirdate/vcf_noheader_patient文件中。如果$length大于1,则先将$vcf文件中的数据提取出来,并将结果存储在$dirdate/vcf_noheader文件中。然后,对于每个病人,从$vcf文件中提取出该病人的信息,并将结果存储在$dirdate/vcf_patient_"$i"文件中。最后,将所有文件合并到$dirdate/vcf_noheader_patient文件中。在处理数据时,使用grep和awk命令来过滤和提取所需的信息。为了提高代码的可读性和性能,使用了变量和循环来简化代码,同时使用了管道和重定向来将数据传递和存储在不同的文件中。
原文地址: https://www.cveoy.top/t/topic/cFDg 著作权归作者所有。请勿转载和采集!