处理hic的pca

paste PCA.WT.50kb.chr19.bedgraph PCA.TU.50kb.chr19.bedgraph | awk 'BEGIN{FS="\t";OFS="\t"}{if ($4>0 && $8>0) print $1,$2,$3,"A2A";else if ($4<0 && $8<0) print $1,$2,$3,"B2B";else if ($4>0 && $8<0) print $1,$2,$3,"A2B"; else if($4<0 && $8>0) print $1,$2,$3,"B2A"; else print $1,$2,$3,"."}' > PCA.WT2TU.50kb.chr19.bedgraph

前提是A/B已经正确矫正！

查看A/B转换的情况：

awk '{print $4}' PCA.WT2TU.50kb.chr19.bedgraph | sort | uniq -c

    109 A2A

    498 A2B

    457 B2A

    96 B2B

输出A2B和B2A：

grep 'A2B' PCA.WT2TU.50kb.chr19.bedgraph > A2B.bed

grep 'B2A' PCA.WT2TU.50kb.chr19.bedgraph > B2A.bed