Bowtie2 比對 參考基因集路徑FA=/refdir/database/server/reference/genome/mm10/mm10....
# ---- 路徑 ----setwd("F:/data/解靜")#---- 包 ----packages <- c("readxl", "dp...
之前的callpeak按200bp read寫的,發(fā)現看錯數據了。。。 從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...
環(huán)境還是m6a的conda activate m6acd /home/data/t210424/ac4c call了peak結果不是很理想暫時沒...
環(huán)境還是m6a的 conda activate m6acd /home/data/t210424/ac4c rRNA測了比例竟然是0%所以沒去 ...