之前的callpeak按200bp read写的，发现看错数据了。。。 从peakcall开始的脚本重新记一下 peaksingle2.sh #...

分组的图 # ==============================# 分组 metagene plot 脚本（Sham vs SNI）#...

library(tidyverse)# 1. 读入counts_ip <- read.table("counts_ip.txt", header...

DESeq2分析出来的结果很烂，换了limma先用IP/INPUT在做组间比较 library(tidyverse) library(limm...

一开始感觉要直接做三组平均，写了个peakaverage.sh 暂时没用上 #!/bin/bash# Define input and outp...

环境还是m6a的conda activate m6acd /home/data/t210424/ac4c call了peak结果不是很理想暂时没...

环境还是m6a的 conda activate m6acd /home/data/t210424/ac4c rRNA测了比例竟然是0%所以没去 ...

环境还是m6a的 conda activate m6acd /home/data/t210424/test QC mkdir ./QCclean...

直接改了/home/data/t210424/m/pcnolambda.sh的output的文件夹 peakcall0423nolambda和p...