首页
学习
活动
专区
工具
TVP
发布

生信技能树

专栏成员
1964
文章
3840901
阅读量
1150
订阅数
生存分析新技能:小洁老师带你全代码实现生存曲线与风险表的完美融合
为了让代码具有可重复性,保存图片也最好是用代码来实现,而不是用点鼠标的方式。最近有一个需求是将生存曲线和表格一起保存,尝试了经典的三段论、ggsave、图片数据类型转换、cowplot包的recordPlo函数都没能实现我的目的。恰好昨天是生信技能树的福利:每月1次的答疑时间。那就求助小洁老师来解决问题吧,小洁老师很快就找到了解决办法,实在是太厉害啦!下面是解答过程:
生信技能树
2024-07-05
970
单细胞转录组可以这样简单计算相关性吗
各种数据挖掘文章本质上都是要把目标基因集缩小,比如表达量矩阵通常是2万多个蛋白编码基因,不管是表达芯片还是RNA-seq测序的,采用何种程度的差异分析,最后都还有成百上千个目标基因。如果是临床队列,通常是会跟生存分析进行交集,或者多个数据集差异结果的交集,比如:多个数据集整合神器-RobustRankAggreg包 ,这样的基因集就是100个以内的数量了,但是仍然有缩小的空间,比如lasso等统计学算法,最后搞成10个左右的基因组成signature即可顺利发表。其实还有另外一个策略方向,有点类似于人工选择啦,通常是可以往热点靠,比如肿瘤免疫,相当于你不需要全部的两万多个基因的表达量矩阵进行后续分析,仅仅是拿着几千个免疫相关基因的表达矩阵即可。最近比较热门的有:自噬基因,铁死亡,EMT基因,核受体基因家族,代谢基因。还有一个最搞笑的是m6a基因的策略,完全是无厘头的基因集搞小,纯粹是为了搞小而搞小。目前单细胞转录组大行其道,所以很多人喜欢使用公共的单细胞转录组数据集来缩小基因范围。学员在微信交流群分享了一个2024年5月的单细胞数据挖掘文章,标题是:《Single-cell combined with transcriptome sequencing to explore the molecular mechanism of cell communication in idiopathic pulmonary fibrosis》,研究者们重新分析了 GSE122960 这个单细胞转录组数据集,主要是第一层次降维聚类分群后,提取了巨噬细胞的特异性基因,然后走了随机森林生存分析算法,得到了 five most related key genes (CD163, IFITM2, IGSF6, S100A14 and SOD3). 有了目标的5个基因就可以很方便的各种简单分析来强调他们的生物学意义。比如去跟PDCD1基因看相关性:
生信技能树
2024-07-05
1220
这配色方案让人费解啊
randomcoloR和paletteer的使用方式类似,都提供了直观的函数来生成和应用颜色方案。randomcoloR 包可以生成随机的颜色方案,非常适合当你需要快速创建一个颜色方案时使用。
生信技能树
2024-07-05
880
月更单细胞图表复现-文献1-第四和五集
但是我们的文字版推文还在第一篇文献,前面已经分享了3个:胃癌单细胞数据集GSE163558复现(二):Seurat V5标准流程,接下来是图表美化和单细胞亚群比例探讨:
生信技能树
2024-07-05
1160
单细胞水平的差异分析是不是就会合理的产出这样的火山图呢
而且很明显,第一层次降维聚类分群其实是没办法区分 28200 epithelial cells (markers: EPCAM, KRT8, and KRT18), 里面的 恶性与否的上皮细胞,这个时候需要走inferCNV等流程对上皮细胞进行恶性程度的判断。
生信技能树
2024-07-05
1010
单细胞转录组差异分析的8大痛点
文献中提出的新范式(new paradigm)是直接在原始UMI计数上应用广义线性混合模型(GLMM),这种方法可以在执行批次校正、标准化、插补或特征选择之前,就对原始数据进行差异表达分析。这种方法保留了样本特定的结构和生物信号,并可以调整任何潜在的混杂因素,如批次、年龄、性别或种族,通过将它们作为固定效应的协变量纳入模型中。这种方法可能有助于克服上述一些挑战,提高单细胞转录组差异分析的准确性和可靠性。
生信技能树
2024-07-05
1440
拆分你的百万级别单细胞数据集后做降维聚类分群
因为纳入的数据集有点多,来源于12篇文章:232 single cell transcriptome samples (normal = 31; adjacent = 54; tumor = 148) ,分别来源于:
生信技能树
2024-07-05
1230
单细胞数据的GSVA
Seurat v5 提示建议用AggregateExpression做伪bulk转录组分析,那个是用来求和的,目前查到的文献和教程都是使用平均值,这里就木有改动.
生信技能树
2024-07-05
1120
单细胞数据的二次分群
Seurat里的FindClusters函数设置的resolution数值越大,分群的数量就越多,但是当单细胞数量太多的时候,会遇到resolution再变大,分群的数量也不再增加的情况。一次分群分不开时就会需要二次分群。
生信技能树
2024-07-05
1040
monocle单样本拟时序分析
做拟时序分析是为了探索自己感兴趣的几种细胞之间的发育关系,一般不是用全部类型的细胞来做的。例如本例中选择了CD14和CD16单核细胞。
生信技能树
2024-07-05
970
monocle多样本拟时序分析
前面已经是介绍了单个样品的单细胞转录组表达量矩阵的monocle分析,接下来分享一下多样品的时候如何注意个体差异因素。
生信技能树
2024-07-05
1000
生信爱好者周刊(第 128 期):什么事情使你放弃了学术理想?
本杂志开源(GitHub: ShixiangWang/weekly[1]),欢迎提交 issue,投稿或推荐生信相关内容。
生信技能树
2024-07-05
820
单细胞组间比较分析
这篇文章介绍的是有分组的单细胞数据怎样分析,数据来自GEO的GSE231920,有3个treat,3个control样本,代码完整,可以自行下载数据跑一跑,但请注意细胞数量是6w,对计算资源要求较高,自己的电脑跑不动,需要在服务器上跑。
生信技能树
2024-07-05
1280
为什么你的R语言不能默认显示英文呢
因为学习新技能,怎么样都要遇到问题并解决问题的,英文的报错信息更容易搜索得到答案,而中文的就不怎么好搜。mac的同学就不用看了,默认就是英文的,不用设置。 而天煞的windows默认就是中文,不设置就难受死啦。
生信技能树
2024-06-21
1100
CellChat细胞通讯
细胞通讯是单细胞数据高级分析中比较常见的一个,我们习惯使用的R包是CellChat。
生信技能树
2024-06-21
2060
单细胞亚群的命名的细致程度需要有一个度
可以看到,这个时候已经是有了 M1 macrophages, natural killer T (NKT) cells, plasma cells 这样的单细胞亚群:
生信技能树
2024-06-19
1420
单细胞亚群合理命名是数据分析基石啊
其实这样的笑话在单细胞数据挖掘文章里面层出不穷,比如另外一个数据挖掘文章;《Identification of Five Hub Genes Based on Single-Cell RNA Sequencing Data and Network Pharmacology in Patients With Acute Myocardial Infarction》,链接是:https://www.frontiersin.org/articles/10.3389/fpubh.2022.894129/full ,里面的bug就更可怕了,居然是把成纤维细胞错误的命名成为了单核细胞。如下所示:
生信技能树
2024-06-08
1600
作者仅提供了fpkm格式表达量矩阵的转录组测序数据集该如何重新分析呢
研究者们在GEO数据库是有数据分享:https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE182923
生信技能树
2024-06-08
1700
差异分析的火山图为什么不喷发呢
数据集是:https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE33113
生信技能树
2024-06-08
1410
代码分享|| 见过将代码整理成wiki资源的吗?这是篇极好的学习单细胞与scATAC-Seq组学以及联合分析的文献资源!
由于最近学习多组学方向的思路,顺便随手将以前整理的资源进行了翻阅,发现了非常多好东西,~本次分享的文献是一篇极好的学习scRNA-Seq与scATAC-Seq组学以及联合分析的文献资源!囊括了上游分析的bash代码和下游每一张Figure的复现代码,是一个非常好的学习单细胞与ATAC-Seq联合分析的文献资源,代码可以说整理的非常好了,还包括各种readme说明。作者甚至还将代码整理了一个wiki版本的,这说明有各种详细的说明以及代码大纲,简直不要太棒!
生信技能树
2024-06-08
1370
点击加载更多
社区活动
【纪录片】中国数据库前世今生
穿越半个世纪,探寻中国数据库50年的发展历程
Python精品学习库
代码在线跑,知识轻松学
博客搬家 | 分享价值百万资源包
自行/邀约他人一键搬运博客,速成社区影响力并领取好礼
技术创作特训营·精选知识专栏
往期视频·千货材料·成员作品 最新动态
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档