你有我有全都有呀!-连续出击cellphonedb v5受配体多组比较气泡图(原创函数)

学术   2024-08-17 10:00   重庆  

偷偷问一下,关注了吗

内容获取


1、购买打包合集(《KS科研分享与服务》付费内容打包集合),价格感人,可以加入微信VIP群(答疑交流群,甚至有小伙伴觉得群比代码更好),可以获取建号以来所有内容,群成员专享视频教程,提前更新,其他更多福利!


2、《KS科研分享与服务》公众号有QQ群,进入门槛是20元(完全是为了防止白嫖党,请理解),请考虑清楚。群里有免费推文的注释代码和示例数据(终身拥有),没有付费内容,群成员福利是购买单个付费内容半价!


需要者详情请联系作者(非需要者勿扰,处理太费时间):


前面我们发布了关于cellchat的函数(连夜更新---别说两组了,这个cellchat多组比较气泡图函数10组也能做了)。因为cellchat比较好入手,所以先开刀。很多小伙伴说有没有cpdb的,其实在写函数之初,我们就考虑到了,只不过先从cellchat好入手,本来以为套用可能大差不差,结果cpdb在数据上有很大出入,所以这次费了点时间。However,最终效果刚刚的!
微信VIP已提前发布,自行下载,需要购买函数可微信联系作者!
函数B站解说视频(一定要看使用方法哦!):
https://www.bilibili.com/video/BV1EreueiE7Q/?spm_id_from=333.999.0.0&vd_source=05b5479545ba945a8f5d7b2e7160ea34

函数主体:也是支持多组,支持自选受配体,自选pathway,自定义分类!

看看演示:load data
library(ggplot2)library(tidyr)
#load dataGO_pvals <- read.delim("./GO_cpdb/statistical_analysis_pvalues_08_15_2024_132104.txt", check.names = FALSE)GO_means <- read.delim("./GO_cpdb/statistical_analysis_means_08_15_2024_132104.txt", check.names = FALSE)

WT_pvals <- read.delim("./WT_cpdb/statistical_analysis_pvalues_08_15_2024_132617.txt", check.names = FALSE)WT_means <- read.delim("./WT_cpdb/statistical_analysis_means_08_15_2024_132617.txt", check.names = FALSE)

data = list(list(pval=GO_pvals, means=GO_means), list(pval=WT_pvals, means=WT_means))
测试1:选定通路
#测试1:cpdb_anno没有pathway,用cpdb默认的,用通路选择cpdb_interLR <- read.csv(file="cpdb_interLR",header = T)#选定pathway,注释文件中没有pathwayks_cpdb_Group_bubble(cpdb_data = data,                     group_names = c("GO","WT"),                     analysis_cells = "Endothelial",                     pathway = c("Signaling by Transforming growth factor","Signaling by Semaphorin"),                     cpdb_anno = cpdb_interLR,                     tag_pos = c(0.5,0.12),                     sig = F)

#随机换种celltype试试ks_cpdb_Group_bubble(cpdb_data = data, group_names = c("GO","WT"), analysis_cells = "Macrophages", pathway = c("Adhesion by Laminin","Signaling by Integrin"), cpdb_anno = cpdb_interLR, tag_pos = c(0.4,0.2), sig = F)
#只显示显著的,sig=Tks_cpdb_Group_bubble(cpdb_data = data, group_names = c("GO","WT"), analysis_cells = "Endothelial", pathway = c("Signaling by Transforming growth factor","Signaling by Semaphorin"), cpdb_anno = cpdb_interLR, tag_pos = c(0.5,0.12), sig = T)

测试2:自选受配体,自定义分类!
#测试2#自选受配体对,注释文件带pathway注释cpdb_interLR_anno <- read.csv(file = 'cpdb_interLR_anno.csv', header = T, row.names = 1)select_LR <- read.csv('plot_pairs.csv', header = F)

ks_cpdb_Group_bubble(cpdb_data = data, group_names = c("GO","WT"), analysis_cells = "Endothelial", select_LR = select_LR$V1, cpdb_anno = cpdb_interLR_anno, tag_pos = c(0.4,0.2), sig = F)

没毛病,非常完美!希望对你有所帮助!



END

觉得我们分享有些用的,点个赞再走呗!

关注我们获取精彩内容:


关注不迷路:扫描下面二维码关注公众号!
B站视频号链接https://space.bilibili.com/471040659?spm_id_from=333.1007.0.0




关注 KS科研分享与服务,

认清正版优质内容和服务!

优质内容持续输出,物超所值!

合作联系:ks_account@163.com

新的板块-重要通知-双向选择

KS科研分享与服务
科研学习交流于分享,生信学习笔记,科研经历和生活!
 最新文章