1.6.1按树结构分组

gs_groups <- gs_split_by_tree(Gslist) - gs1, gs2, gs3, gs4, gs5

##按门控树分组…

长度(gs_groups)

[1] 4

这就把所有GatingSetS分成不同的组，每个组共享相同的树结构。这里有4##检查是否可以通过删除叶节点来解决差异

res <- try(gs_check_redundant_nodes(gs_groups)， silent = TRUE) print(res[[1]])

[1] "错误在(函数(thisNodeSet, thisObj):不能删除非终端节点:singlets "

显然，非叶节点(单线态)检查失败，由用户决定是否隐藏该节点或保持该组分离，以免进一步合并。在这里，我们试图隐藏它。

对于(gp in gs_groups) plot(gp[[1]])

基于每个组的树形结构(通常没有像GatingSet对象本身)，我们将隐藏单线态为组2而且组4．

for(i in c(2,4)) for(gs in gs_groups[[i]]) invisible(gs_pop_set_visibility(gs， "singlets"， FALSE))

现在再检查一下.gs_check_redundant_nodes

toRm <- gs_check_redundant_nodes(gs_groups) toRm

[[1]] [1] " ccr7 + 45ra + " " ccr7 + 45ra - "

[[2]] [1] " dnt " " dpt "

[[3]]字符(0)

[[4]]字符(0)

基于此，可以通过删除*来合并这些组CCR7 + 45 ra +而且CCR7 + 45 ra -从组1．＊二硝基酚而且DPT从组2．

有时，通过简单地绘制整个门通树或查看这个简单的平面节点列表(特别是当某些组中缺少整个子树时)，可能很难检查和区分树差异。使用helper函数只可视化和突出显示树层次结构的差异是有帮助的gs_plot_diff_tree

gs_plot_diff_tree (gs_groups)

为了继续删除这些节点，.gs_remove_redundant_nodes可以用来代替手工做吗

gs_remove_redundant_nodes (gs_groups暴雨)

##去除CCR7+ 45RA+

##去除CCR7+ 45RA-

##删除DNT

##删除DPT

现在它们可以合并为一个GatingSetList．

GatingSetList (gslist)

##警告:“GatingSetList”已弃用。##使用'merge_list_to_gs'代替。##查看帮助(“已弃用”)

一个包含5个唯一样本的GatingSetList。