目录 1. 加载表达矩阵 2. 创建样本分类表 3. 加载limma包,构建如下矩阵 4. 规定哪一组数据与哪一组数据比较 5. 获取差异表达数据 6. 导出差异表达数据 1. 加载表达矩阵 2. 创建样本分类表 获得类似 ...
.概述 矩阵 函数 表达矩阵 lmFit 分组矩阵 eBayes 差异比较矩阵 topTable .读取表达矩阵: 得到表达矩阵exprSet,它的列是各个样本名称,行是各个探针ID,一个纯粹的表达矩阵,必须是数字型的 可以简单地做一下该表达矩阵的QC检测: 如果这些boxplot发现各个芯片直接数据还算整齐,则可以进行差异比较。 .制作分组矩阵 design就是分组矩阵,需要根据我们下载的芯片 ...
2020-03-09 14:45 0 697 推荐指数:
目录 1. 加载表达矩阵 2. 创建样本分类表 3. 加载limma包,构建如下矩阵 4. 规定哪一组数据与哪一组数据比较 5. 获取差异表达数据 6. 导出差异表达数据 1. 加载表达矩阵 2. 创建样本分类表 获得类似 ...
test.txt:文件格式: gene TCGA-3M-AB46-01A-11R-A414-31 TCGA-3M-AB47-01A-2 ...
简单使用limma做差异分析 Posted: 五月 12, 2017 Under: Transcriptomics By Kai no Comments 首先需要说明的是,limma是一个非常全面的用于分析芯片以及RNA-Seq的差异分析,按照其文章所说: limma ...
仍然是两年前的笔记 1. prepare-reference 如果用RSEM对比对后的bam进行转录本定量,则在比对过程中要确保比对用到的索引是由rsem-prepare-reference产生的。 可以看到,单纯用bowtie2建的索引和rsem调用bowtie2建的索引 ...
差异表达分析之FDR 随着测序成本的不断降低,转录组测序分析已逐渐成为一种很常用的分析手段。但对于转录组分析当中的一些概念,很多人还不是很清楚。今天,小编就来谈谈在转录组分析中,经常会遇到的一个概念FDR,那什么是FDR?为什么要用FDR呢?一起来学习吧! 什么是FDR ...
接前一篇: 用R和BioConductor进行基因芯片数据分析(五):芯片间归一化 经过一系列的预处理,包括缺失值填充,中位数计算以及归一化,我们的数据终于可以用啦。 下面我们就来分析一下new population和old population的个体是否有差异表达基因。 判断一个基因是否 ...
前言 本文主要演示GEO数据库的一些工具,使用的数据是2015年在Nature Communications上发表的文章Regulation of autophagy and the ubiquit ...
...