생물정보학 분석에서 fpkm을 계산하기 위해 awk를 사용하는 방법은 무엇입니까? [폐쇄]

생물정보학 분석에서 fpkm을 계산하기 위해 awk를 사용하는 방법은 무엇입니까? [폐쇄]

이것은 내 코드입니다.

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

오류는 정확히 무엇이 잘못되었는지 알려줍니다. 즉, 수학적으로 불법인 0으로 나누려고 합니다. 이는 파일의 2번째 라인에서 발생합니다 BMXV02.counts. 코드의 유일한 구분은 다음 1000000*1000*$NF/($(NF-1)*counts)이므로 $(NF-1)*counts분명히 반환됩니다 0.

아마도 $(NF-1)2행의 두 번째 필드( )는 a 입니다 0.

관련 정보