![バイオインフォマティクス分析でfpkmを計算するためにawkを使用する方法は？ [閉鎖]](https://linux33.com/image/207992/%E3%83%90%E3%82%A4%E3%82%AA%E3%82%A4%E3%83%B3%E3%83%95%E3%82%A9%E3%83%9E%E3%83%86%E3%82%A3%E3%82%AF%E3%82%B9%E5%88%86%E6%9E%90%E3%81%A7fpkm%E3%82%92%E8%A8%88%E7%AE%97%E3%81%99%E3%82%8B%E3%81%9F%E3%82%81%E3%81%ABawk%E3%82%92%E4%BD%BF%E7%94%A8%E3%81%99%E3%82%8B%E6%96%B9%E6%B3%95%E3%81%AF%EF%BC%9F%20%5B%E9%96%89%E9%8E%96%5D.png)
これは私のコードです。
awk -v "counts=170712" '{if(NR>1){print $1 "\t" 1000000*1000*$NF/($(NF-1)*counts)} else{print $1"\t"$NF}}' a.counts > a.fpkm.txt
ただし、次のエラーが発生します。
awk: cmd. line:1: (FILENAME=BMXV02.counts FNR=2) fatal: division by zero attempted
何を変えるべきですか？
答え1
エラーは正確に何が間違っているかを伝えます。つまり、数学的に違法なゼロに分けようとします。これはファイルの2行目で発生しますBMXV02.counts。コードの唯一の区別は次1000000*1000*$NF/($(NF-1)*counts)のため$(NF-1)*counts、明らかに返されます0。
おそらく$(NF-1)2行目の2番目のフィールド（）はaです0。