# --quantMode TranscriptomeSAM will output alignments translated into tran-script coordinates,为了使用RSEM 进行定量分析做准备;# STAR align 复杂版本 样本R2 STAR --runThreadN 12 --genomeDir arab_STAR_genome \ --readFilesCommand zcat \ --readFilesIn 02clean_data/sample2_1_paired_clean_R1....
1.在Rstudio中安装DESeq2 if(!requireNamespace("BiocManager",quietly=TRUE))install.packages("BiocManager")library(BiocManager)install.packages("DESeq2") 2.分析 ##批量处理library(DESeq2)library(stringr)library(dbplyr)file_name<-c("L7_24h-vs-L7_CK","L7_24h-vs-L7_Re24h","L7_24h-vs-Le_...
STAR+RSEM 先比对再进行定量 # build index 建立索引;STAR--runThreadN10--runMode genomeGenerate \--genomeDir arab_STAR_genome \--genomeFastaFiles00ref/TAIR10_Chr.all.fasta \--sjdbGTFfile00ref/Araport11_GFF3_genes_transposons.201606.gtf \--genomeSAindexNbases12\--sjdbOverhang149# 一般为reads...