Counts到差异表达基因(二)- 有重复数据 - DiffAnalysis.r

写在前面

在《原始数据到Counts数据》中,已经讲明白了,如何得到Counts表格。
总之,运行结束,你应该会得到以下文件

ls -1 *.csv
gene_count_matrix.csv
transcript_count_matrix.csv

其中一个是gene水平的Counts数据,一个是转录本水平的。除非有特殊要求,一般我们只使用基因水平的Counts。
那么接下来就是进行差异表达分析

准备输入文件

  1. Counts文件
    CSV格式,在上述我们已经准备好了,大体如下



    2.样本信息文件,sample.info
    注意,header必须是Sample type condition



    样本信息,应该是将会用到的所有样本,比如(编一个)

3.两两比较的信息,comp.info
注意,header必须是Sample1 Sample2



这个文件比较重要,可以一次设置无限个比较方案,比如



这个输入,就会进行7个两两比较的差异表达分析

进行差异表达分析

进行差异表达分析,比较简单...(如果报错,找Junting Feng

Rscript /tools/XiaLabRNAseqPipe/DiffAnalysis.r gene_count_matrix.csv sample.info comp.info

等着,等到最后就完成了


输出的文件比较多,包括差异表达分析结果和PCA等


写在最后

看不懂结果的,找JuntingFeng

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容

  • Spring Cloud为开发人员提供了快速构建分布式系统中一些常见模式的工具(例如配置管理,服务发现,断路器,智...
    卡卡罗2017阅读 135,027评论 19 139
  • 原文:http://blog.sina.com.cn/s/blog_751bd9440102v72b.html 一...
    简单点lili阅读 4,160评论 0 8
  • 用两张图告诉你,为什么你的 App 会卡顿? - Android - 掘金 Cover 有什么料? 从这篇文章中你...
    hw1212阅读 13,045评论 2 59
  • 我的母亲 世界上有一种最动听的声音,那便是母亲的呼唤 。小时的我,很贪玩。一放学,就喜欢到佛岭街玩个够,不到妈妈...
    深圳彭秀清阅读 562评论 0 3
  • 2007年,国内电子商务最热门的关键词应该是“PPG”了。PPG在国内电子商务企业中几乎是一个神话,好像一夜之间,...
    包云东阅读 618评论 0 4