【科研猫·出品】TCGA超大批量生存分析教程
上一期的教程【科研猫】生存分析的正确姿势(1)视频+R代码分享,给大家讲解了批量对TCGA中的基因进行生存分析的第一步。教大家如何通过cBioportal去下载TCGA中的数据和患者的表型数据,然后通过R语言做了一个生存分析的KM图。
最近很多同学都在催我们啦,何时出第二季,何时教大家进行批量的生存分析。被催到脸红了呀,实在是羞愧,最近团队给接了太多数据分析的项目。【小编插话:每天跟这些晚期强迫症技术要稿子,我容易嘛我。。。】

大部分科研的套路都是从一个表型出发,筛出一个有意义的基因,然后对其进行深入的挖掘和功能实验。但是如何才能筛选到真正重要的目标基因呢?这个问题困扰着很多人。作为一名临床科研工作者,如果一个分子能够跟患者预后扯上关系,那么文章的档次立刻就上来了。如果我们能分分钟完成数百上千个基因的生存分析,那么以后的“筛基因”“找靶点”是不是就有了指路明灯了呢?有木有一点小激动呢~
经过我们科研猫团队生信分析团队的数日辛苦劳作,我们完成了全网首个批量进行TCGA生存分析的程序,只要输入你想分析的基因列表,下载相关肿瘤的数据文件,运行代码即可。我把这个程序命名为Survival_Analysis_Terminator.R,没错又是“Terminator系列”,为什么命名为“终结者”呢?因为有了这个代码,终结所有相关问题,无需求助其他软件。
为了演示,我在这里提交了500个基因,总运行时长仅3分多钟!是不是很心潮澎湃嗯?先来看个操作演示视频,这次配了音效,还挺燃的呢~
【科研猫】TCGA批量生存分析-代码-教学视频
代码运行后,会生成如下的三个结果文件:
1) 1.LogRank.Pvalues.csv,包括500个基因的生存分析Log-rank test的P值。
2) 2.Kaplan-Meier_plots.pdf,500个基因的KM图。你没看错,是一个500张高清矢量大图的pdf文件。
3) 3.KM_plots_with_significant_Pvalues文件夹,这个文件夹中包括数据所有log-rank test之后P值显著的所有基因的KM图,是tiff格式,方便大家写文章使用,可直接放到paper里面。

这三个文件和文件夹的截图如上图所示。所有分析都是自动化完成,无需修改任何代码。

为了完成这个代码,【科研猫】的整个生信团队付出了诸多心血,要知道这样的一个代码,就可以作为某些小型生物公司数据分析服务的技术核心。包括之前的所有代码和教程(GEO数据库挖掘、功能富集分析气泡图等),如果外包给公司做成case,每次分析费用至少2000+。懂行的小伙伴纷纷夸赞咱们公众号良心!

可是,让人气愤的是,不少同学反映网络上出现了多个盗版科研猫的教程。由于教程和代码都是免费获取的,竟然有人随便换个抬头就堂而皇之地转手高价出售(淘宝和微信已各警告过一家,请侵权者好自为之;再若盗版,必追究法律责任)。还有小白同学竟然来问,咱们不收费是不是因为是盗版!
盗版也就算了,盗版发了之后没有答疑,不会解答大家的疑问,导致很多同学在学习过程中屡屡受挫乃至怀疑起咱们的代码来。着实坏了我们“科研猫“的招牌,实在让人心累。因此,和我们的生信团队反复讨论后,我们决定从技术部拨出来部分人力提供网上教程答疑服务。
另外,仰仗各位同仁转发朋友圈一直是咱们公众号的主要推广方式,但由于交际环境的复杂,渐渐考虑到全部可见的诸多不便,经反复考量和商议,综合此次代码的高含金量:
此次采用付费服务的方式,初步仅定价298元,具体服务打包如下:
市价5000元的批量生存分析代码和教学视频;为期3个月VIP答疑服务
生信分析师1对1技术指导及支持;100元生信分析优惠券(有效期2年)
提供试剂/服务费发票(专票/普票均可,15日内寄出)
科研猫的推文怎么会没有福利呢?看到这篇半个广告贴,会不会有点失望,
不会!不能!绝不要!
福利通道
1)推文投稿:一经审核发表,该代码免费赠送;但凡对诸位科研同道有所裨益之论,皆可投稿,具体事宜,联系客服小姐姐(下方二维码)。
2)凡与科研猫所属 上海腾云生物科技有限公司 ,进行过合作、委托过数据分析及论文合作者,免费赠送。(需提供发票或合作合同)
3)凡转发过【科研猫】公众号往期推文,保留至今未删除者,可享受8折。(需大于15日)
4)朋友圈打卡集赞,全部可见,集满30赞,可享受7折。
欲购买该套分析套餐的同学,赶紧联系咱们折耳猫小姐姐吧。
欲购买该套分析套餐的同学,赶紧联系咱们折耳猫小姐姐吧。
更多科研新鲜资讯、文献精度和生物信息技能,请关注科研猫公众号
下方点好看,更多好看。
往期热文:
生存分析
挖掘GEO速成SCI文章系列教程(代码免费分享)
重磅:GEO数据库挖掘教程(4)一体化分析代码(带视频+R代码分享)
GO/KEGG功能富集系列(代码免费分享)
TCGA数据挖掘系列
隔壁实验室的“秃鹫”师兄又发SCI啦--TCGA数据挖掘实战
生物信息入门系列(书籍、R包免费分享)