偷偷问一下,关注了吗?
内容获取
1、购买打包合集(《KS科研分享与服务》付费内容打包集合),价格感人,可以加入微信VIP群(答疑交流群,甚至有小伙伴觉得群比代码更好),可以获取建号以来所有内容,群成员专享视频教程,提前更新,其他更多福利!
2、《KS科研分享与服务》公众号有QQ群,进入门槛是20元(完全是为了防止白嫖党,请理解),请考虑清楚。群里有免费推文的注释代码和示例数据(终身拥有),没有付费内容,群成员福利是购买单个付费内容半价!
需要者详情请联系作者(非需要者勿扰,处理太费时间):
Monocle3往期
Monocle3(1):单细胞拟时分析---分析(详细注释版)
##run
cds_data <- ks_run_monocle3(object=cytotrace2_sce,
idents="cluster",
use_partition=F,
learn_graph_control=list(minimal_branch_len=9.5,euclidean_distance_ratio=1),
define_root=T,
know_root=T,
root_state="YSMP")
#plot
label_cell_groups = F,
color_cells_by = "pseudotime",
label_branch_points = F,
label_roots =F,
label_leaves =F,
graph_label_size = 0,
cell_size=2,
trajectory_graph_color='black',
trajectory_graph_segment_size = 2)
#在monocle3分析中,我们使用graph_test函数分析过拟时差异基因,这个过程比较慢
modulated_genes <- graph_test(cds_data, neighbor_graph = "principal_graph", cores = 5)
plot_genes <- row.names(subset(modulated_genes, q_value < 0.01 & morans_I > 0.25))
#plot
ht1 = ks_monocle3_heatmap(cds = cds_data,
graph_gene = plot_genes,
celltype_color = c("#E69253", "#EDB931", "#E4502E", "#4f372d"),
cluster_color =c("#8f657d", "#42819F", "#86AA7D", "#CBB396"),
num_clusters=4)
#label gene
genes <- c("C1QB","C1QC","C1QA","MRC1","LGMN","MS4A7","MAF","FOLR2",
"HLA-DPA1","CLEC10A","IL10RA","CD163","KCTD12","CLEC7A","MS4A6A","CD14",
"ITM2A","CYTL1","MDK","SELP","CD24",
"S100A8","S100A9","S100A12")
ht2 = ks_monocle3_heatmap(cds = cds_data,
graph_gene = plot_genes,
celltype_color = c("#E69253", "#EDB931", "#E4502E", "#4f372d"),
cluster_color =c("#8f657d", "#42819F", "#86AA7D", "#CBB396"),
num_clusters=4,
labels = T,
label_genes = genes)
觉得我们分享有些用的,点个赞再走呗!
关注我们获取精彩内容:
关注不迷路:扫描下面二维码关注公众号!
B站视频号链接:https://space.bilibili.com/471040659?spm_id_from=333.1007.0.0
关注 KS科研分享与服务,
认清正版优质内容和服务!
优质内容持续输出,物超所值!
合作联系:ks_account@163.com