## ----------------------------------------------------------------------------- suppressPackageStartupMessages({库(ExperimentHub)图书馆(SingleCellExperiment)图书馆(DuoClustering2018)库(plyr ) }) ## ----------------------------------------------------------------------------- 嗯< - ExperimentHub()查询(呃,“DuoClustering2018 ") ## ----------------------------------------------------------------------------- 嗯”[[EH1500 "]] ## ----------------------------------------------------------------------------- sce_filteredExpr10_Koh () ## ----------------------------------------------------------------------------- res < - plyr:: rbind。填充(clustering_summary_filteredExpr10_Koh_v2 (), clustering_summary_filteredHVG10_Koh_v2 (), clustering_summary_filteredExpr10_Zhengmix4eq_v2 (), clustering_summary_filteredHVG10_Zhengmix4eq_v2())的(res ) ## ----------------------------------------------------------------------------- 头(res ) ## ----------------------------------------------------------------------------- method_colors < - c (CIDR =“332288 #”,FlowSOM = " # 6699 cc " PCAHC = " # 88 ccee PCAKmeans = " # 44 aa99 ",pcaReduce = " # 117733 ", RtsneKmeans = " # 999933”,修= " # DDCC77 " SC3svm =“661100 #”,SC3 = " # CC6677 " TSCAN =“grey34”=“橙色”,提升安全=“黑人”,单片眼镜=“红色”,RaceID2 = "蓝色 ") ## ----------------------------------------------------------------------------- 性能< - plot_performance (res,method_colors = method_colors)名称(穿孔)穿孔median_ari_vs_k perf median_ari_heatmap_truek美元 ## ----------------------------------------------------------------------------- 刺< - plot_stability (res method_colors = method_colors)名称(刺)刺stability_vs_k刺stability_heatmap_truek # #美元——fig.height = 8 ---------------------------------------------------------- entr < - plot_entropy (res,method_colors = method_colors)名称(entr) entr $ entropy_vs_k entr normentropy # #——fig.height = 8美元 ---------------------------------------------------------- 时间< - plot_timing (res method_colors = method_colors scaleMethod =“RtsneKmeans”)的名字(时机)时机time_normalized_by_ref # #——fig.height = 8美元 ---------------------------------------------------------- kdiff < - plot_k_diff (res,method_colors = method_colors)名称(kdiff) kdiff diff_kest_ktrue美元 ## ----------------------------------------------------------------------------- sessionInfo ()