RSEM-Ebseq-差异表达分析-无参

1. RSEM的安装和使用：

$ tar -xzf RSEM-1.3.0.tar.gz

$ cd RSEM-1.3.0

$ make

$ make ebseq

$ echo 'PATH=$PATH:/.../' >> ~/.bashrc

$ source ~/.bashrc

$ extract-transcript-to-gene-map-from-trinity trinity.fa gene_map （可选，转录本和基因都做）

$ rsem-prepare-reference --bowtie2 ref_trinity.fasta ref_trinity (I. Preparing Reference Sequences)

$ rsem-prepare-reference --transcript-to-gene-map ene_map --bowtie2 ref_trinity.fasta ref_trinity（可选，转录本和基因都做）

$ rsem-calculate-expression -p 48 --paired-end --bowtie2 read1.fq read2.fq ref_trinity sample.name (II. Calculating Expression Values)

$ rsem-plot-model sample.name sample.name.diagnostic.pdf （可选，看一些统计）

2. 差异表达：

$ rsem-generate-data-matrix sampleA.results sampleB.results ... > samples.matrix

$ rsem-run-ebseq samples.matrix 3,3 dge.result （两个组，每个组3个重复，如2,,3,3为3个组分别有2/3/3个重复）

$ rsem-control-fdr dge.result 0.05 dge.filter.result (controlling the false discovery rate (FDR) at level 0.05) (PostFC is recommended over the RealFC)

转载于:https://www.cnblogs.com/shawn2018/p/9186087.html

RSEM-Ebseq-差异表达分析-无参相关推荐

rsem比对_无参转录组分析：使用 Trinity 进行转录本拼接（参考脚本）
1,参考脚本: nohup /home/zxd/software/trinityrnaseq-Trinity-v2.4.0/Trinity --seqType fq --max_memory 4G - ...
分析方法升级三代测序辅助，优化无参转录组测序策略
分析方法升级&三代测序辅助,优化无参转录组测序策略无参转录组拼接升级 Corset 让"基因"概念更准确在无参转录组项目中,利用主流软件 Trinity 进行 De n ...
宏基因组有参分析和无参分析差异
宏基因组有参分析和无参分析差异分析流程解决问题结果差异宏基因组流程综述本文参考宏基因组教程Metagenomics Tutorial (HUMAnN2) 分析流程有参流程:质控–物种组成 ...
三天实现独立分析宏基因组数据(有参、无参和分箱等)
在广大粉丝的期待下,<生信宝典>联合<宏基因组>在2019年11月1-3日,北京鼓楼推出<宏基因组分析>专题培训第六期,为大家提供一条走进生信大门的捷径.为同行提供 ...
edger多组差异性分析_使用edgeR进行无重复差异表达分析
写这篇文章一部分原因是填2年前的一个坑转录组入门(7):差异表达分析. 另一部分原因是GQ最近又在搞一波无重复的差异表达分析, 所以专门去学了edgeR 我个人是不太推荐没有重复的差异表达分析,毕竟 ...
基因差异表达分析——基于RSEM对比，DESeq2操作实例
以下操作全部基于R中的DESeq2函数包使用DESeq2的两点要求: DEseq2要求输入数据是由整数组成的矩阵. DESeq2要求矩阵是没有标准化的. 下面是基本流程: 一.准备工作: 确定工作路 ...
基因表达分析（上）- 差异表达分析
基因表达什么是基因表达,如下是来自于维基百科的解释: Gene expression is the process by which information from a gene is used ...
差异表达分析之FDR
差异表达分析之FDR 随着测序成本的不断降低,转录组测序分析已逐渐成为一种很常用的分析手段.但对于转录组分析当中的一些概念,很多人还不是很清楚.今天,小编就来谈谈在转录组分析中,经常会遇到的一个概念F ...
Java类（无参、有参构造方法、静态变量）
构造方法作用: 初始化成员变量的注意: 1.构造方法没有返回值甚至不写返回值的类型 2.构造方法的方法名与类名相同书写: 关键字构造方法名(){ } 举例: 创建cat类: publi ...

RSEM-Ebseq-差异表达分析-无参

RSEM-Ebseq-差异表达分析-无参相关推荐

最新文章

热门文章