作业要求: 我们统一选择p<0.05而且abs(logFC)大于一个与众的基因为显著差异表达基因集,对这个基因集用R包做KEGG/GO超几何分布检验分析。 然后把表达矩阵和分组信息分别作出cls和gct文件,导入到GSEA软件分析。 基本任务是完成这个分析。 【1】环境 ...
作业要求: 使用R语言,载入表达矩阵,然后设置好分组信息,统一用DEseq 进行差异分析,当然也可以走走edgeR或者limma的voom流程。 基本任务是得到差异分析结果,进阶任务是比较多个差异分析结果的异同点。 安装DESeq DESeq 对于输入数据的要求: .DEseq 要求输入数据是由整数组成的矩阵。 .DESeq 要求矩阵是没有标准化的。 DESeq 进行差异表达分析 DESeq 分析 ...
2018-07-03 21:47 0 5049 推荐指数:
作业要求: 我们统一选择p<0.05而且abs(logFC)大于一个与众的基因为显著差异表达基因集,对这个基因集用R包做KEGG/GO超几何分布检验分析。 然后把表达矩阵和分组信息分别作出cls和gct文件,导入到GSEA软件分析。 基本任务是完成这个分析。 【1】环境 ...
引入clusterProfiler与注释数据 GO(gene ontology)分析 GO,Gene Ontology,是基因功能国际标准分类体系。它旨在建立一个适用于各种物种的,对基因和蛋白质功能进行限定和描述的,并能随着研究不断深入而更新的语言词汇标准。GO分为分子功能 ...
edgeR的介绍 背景 RNA-seq表达谱与生物复制的差异表达分析。 实现一系列基于负二项分布的统计方法,包括经验贝叶斯估计,精确检验,广义线性模型和准似然检验。 与RNA-seq一样,它可用于产生计数的其他类型基因组数据的差异信号分析,包括ChIP-seq,Bisulfite-seq ...
任务列表 1.在UCSC下载hg19参考基因组; 2.从gencode数据库下载基因注释文件,并且用IGV去查看感兴趣的基因的结构,比如TP53,KRAS,EGFR等等。 3.截图几个基因的IGV可视化结构 4.下载ENSEMBL,NCBI的gtf,也导入IGV看看,截图 ...
转录组差异表达分析小实战(二) Posted: 八月 14, 2017 Under: Transcriptomics By Kai no Comments 差异基因表达分析 我按照前面的流程转录组差异表达分析小实战(一),将小鼠的4个样本又重新跑了一遍,从而获得一个新 ...
Ballgown是分析转录组差异表达的R包。 软件安装: 运行R, source(“http://bioconductor.org/biocLite.R”) biocLite(“ballgown”) R会自动安装Ballgown,及相应的依赖包。 Ballgown的输入文件 ...
转录组差异表达分析小实战(一) Posted: 七月 28, 2017 Under: Transcriptomics By Kai no Comments 读文献获取数据 文献名称:AKAP95 regulates splicing through ...
作业要求: 在UCSC下载hg19参考基因组,我博客有详细说明,从gencode数据库下载基因注释文件,并且用IGV去查看你感兴趣的基因的结构,比如TP53,KRAS,EGFR等等。 作业,截图几个基因的IGV可视化结构!还可以下载ENSEMBL,NCBI的gtf,也导入IGV看看,截图基因结构 ...