![생물정보학 분석에서 fpkm을 계산하기 위해 awk를 사용하는 방법은 무엇입니까? [폐쇄]](https://linux55.com/image/207992/%EC%83%9D%EB%AC%BC%EC%A0%95%EB%B3%B4%ED%95%99%20%EB%B6%84%EC%84%9D%EC%97%90%EC%84%9C%20fpkm%EC%9D%84%20%EA%B3%84%EC%82%B0%ED%95%98%EA%B8%B0%20%EC%9C%84%ED%95%B4%20awk%EB%A5%BC%20%EC%82%AC%EC%9A%A9%ED%95%98%EB%8A%94%20%EB%B0%A9%EB%B2%95%EC%9D%80%20%EB%AC%B4%EC%97%87%EC%9E%85%EB%8B%88%EA%B9%8C%3F%20%5B%ED%8F%90%EC%87%84%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
오류는 정확히 무엇이 잘못되었는지 알려줍니다. 즉, 수학적으로 불법인 0으로 나누려고 합니다. 이는 파일의 2번째 라인에서 발생합니다 BMXV02.counts
. 코드의 유일한 구분은 다음 1000000*1000*$NF/($(NF-1)*counts)
이므로 $(NF-1)*counts
분명히 반환됩니다 0
.
아마도 $(NF-1)
2행의 두 번째 필드( )는 a 입니다 0
.