连夜更新---别说两组了,这个cellchat多组比较气泡图函数10组也能做了

学术   2024-08-16 10:01   重庆  

偷偷问一下,关注了吗

内容获取


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


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


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


事情是这样的:前面我们发布了这样一款函数((原创函数)-Cellchat受配体多组比较气泡图函数)。但是限制是最多只能进行三组的比较。小伙伴也强烈要求更新多组分析。其实在第一版函数完成后,我们就考虑到了这个问题,也有这个意向,并在睡觉前冥想了一下思路,所以比较顺畅,连夜进行了函数更新!
看看更新函数:此次更新不仅可以进行多组比较了,还能够自选受配体!

具体的使用和更新请看视频吧,可能比写来的更清楚!(注意:购买过此函数的小伙伴记得微信联系作者免费获取更新哦!)

视频解说:点击即可播放!

(背景音乐:长寿六宝;视频解说:TS的美梦)
效果如下:
setwd('D:\\KS项目\\公众号文章\\Cellchat受配体多组比较气泡图函数\\更新-cellchat多组受配体可视化函数(不限组数)')# group_name <- c("HD","MDA")# analysis_cells <- "Tcell"# pathway=c("CCL", "MIF","MHC-II")#两组object_list1 <- list(HD.cellchat, MDA.cellchat)ks_ccGroup_bubble(object_list=object_list1,                  analysis_cells = "Tcell",                  pathway=c("CCL", "MIF","MHC-II"),                  group_name = c("HD","MDA"),                  tag_pos = c(0.5, 0.085))
#三组-假装有三组object_list2 <- list(HD.cellchat, MDA.cellchat,HD.cellchat)ks_ccGroup_bubble(object_list=object_list2, analysis_cells = "Tcell", pathway=c("CCL", "MIF","MHC-II"), group_name = c("HD","MDA","EE"), tag_pos = c(0.5, 0.085))

#4组-假装有4object_list3 <- list(HD.cellchat, MDA.cellchat,HD.cellchat,MDA.cellchat)ks_ccGroup_bubble(object_list=object_list3, analysis_cells = "Tcell", pathway=c("CCL", "MIF","MHC-II"), group_name = c("HD","MDA","EE","AA"), tag_pos = c(0.5, 0.085))


#选定LRobject_list2 <- list(HD.cellchat, MDA.cellchat,HD.cellchat)LR <- read.csv('hh.csv', header = F)
ks_ccGroup_bubble(object_list=object_list2, analysis_cells = "Tcell", select_LR=LR$V1, group_name = c("HD","MDA","EE"), tag_pos = c(0.55, 0.085))



END

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

关注我们获取精彩内容:


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




关注 KS科研分享与服务,

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

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

合作联系:ks_account@163.com

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

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