偷偷问一下，关注了吗？
内容获取
1、购买打包合集（《KS科研分享与服务》付费内容打包集合），价格感人，可以加入微信VIP群（答疑交流群，甚至有小伙伴觉得群比代码更好），可以获取建号以来所有内容，群成员专享视频教程，提前更新，其他更多福利！
2、《KS科研分享与服务》公众号有QQ群，进入门槛是20元（完全是为了防止白嫖党，请理解），请考虑清楚。群里有免费推文的注释代码和示例数据（终身拥有），没有付费内容，群成员福利是购买单个付费内容半价！
需要者详情请联系作者（非需要者勿扰，处理太费时间）：
DEGs <- FindMarkers(human_data, ident.1 = "BM",ident.2 = "GM",logfc.threshold = 0.25,min.pct = 0.25, group.by = 'group')
exp_m <- GetAssayData(human_data, layer = 'counts')gene.all <- rownames(exp_m)mt.gene = grep('^MT-', gene.all, value = TRUE)#线粒体基因ribosome.gene = grep('^RPL|^RPS|^MRPS|^MRPL', gene.all, value = TRUE)#核糖体基因features = setdiff(gene.all, c(mt.gene, ribosome.gene))#去除不需要的基因
DEGs <- FindMarkers(human_data, ident.1 = "BM",ident.2 = "GM",logfc.threshold = 0.25,min.pct = 0.25, group.by = 'group',features = features)
觉得我们分享有些用的，点个赞再走呗！
关注我们获取精彩内容：
关注不迷路：扫描下面二维码关注公众号！
B站视频号链接：https://space.bilibili.com/471040659?spm_id_from=333.1007.0.0
关注 KS科研分享与服务，
认清正版优质内容和服务！
优质内容持续输出，物超所值！
合作联系：ks_account@163.com
