發(fā)簡(jiǎn)信
IP屬地：廣東
Bowtie2 比對(duì) 參考基因集路徑FA=/refdir/database/server/reference/genome/mm10/mm10....
# ---- 路徑 ----setwd("F:/data/解靜")#---- 包 ----packages <- c("readxl", "dp...
之前的callpeak按200bp read寫(xiě)的，發(fā)現(xiàn)看錯(cuò)數(shù)據(jù)了。。。 從peakcall開(kāi)始的腳本重新記一下 peaksingle2.sh #...
分組的圖 # ==============================# 分組 metagene plot 腳本（Sham vs SNI）#...
library(tidyverse)# 1. 讀入counts_ip <- read.table("counts_ip.txt", header...
DESeq2分析出來(lái)的結(jié)果很爛，換了limma先用IP/INPUT在做組間比較 library(tidyverse) library(limm...
一開(kāi)始感覺(jué)要直接做三組平均，寫(xiě)了個(gè)peakaverage.sh 暫時(shí)沒(méi)用上 #!/bin/bash# Define input and outp...
環(huán)境還是m6a的conda activate m6acd /home/data/t210424/ac4c call了peak結(jié)果不是很理想暫時(shí)沒(méi)...
環(huán)境還是m6a的 conda activate m6acd /home/data/t210424/ac4c rRNA測(cè)了比例竟然是0%所以沒(méi)去 ...