我有具有数百个节点和长标签的hclust层次化集群对象。例如,一个家族中多个基因的同义词。见下文。
我想把hclust切成更小的子树,然后用灵活的风格将它们可视化。在http://gastonsanchez.com/blog/how-to/2012/10/03/Dendrograms.html之后,我看到了如何切割树状图和漂亮的类人猿进化树.
我看不出有什么方法可以把切割的树状图转换成phylo对象。
> as.phylo(as.dendrogram(hc))
Error in UseMethod("as.phylo") :
no applicable method for 'as.phylo' applied to an object of class "dendrogram"我对任何呈现圆形或垂直方向子树的方法都是开放的。
--事实上,我的目标是直观地检测基因同义词中的模式,这样我就可以为它们编写类似于 八字胡 模板的东西,所以我甚至对不涉及树状图的解决方案持开放态度。有几篇关于纯文本的多序列比对的帖子,但它们有点超出了我的头脑。
> receptor.synonyms
synonym
1 alpha1B-adrenergic receptor
2 B1AR
3 adrenergic receptor, alpha 2a
4 beta 3-AR
5 alpha-2AAR
6 alpha2-C4
7 Adrb-1
8 Badm
9 beta 1-AR
10 Adrenergic, alpha2C-, receptor class I
11 alpha-1D adrenoceptor
12 beta 2-AR
13 adrenergic receptor
14 alpha-2A-adrenergic receptor
15 Adrenergic, alpha2B-, receptor class III
16 adrenergic, alpha 1B, receptor
17 α<sub>2</sub>-C2
18 adrenergic, alpha-1A-, receptor
19 ADRARL1
20 alpha-1B adrenoceptor
--- snip ---

https://stackoverflow.com/questions/31972832
复制相似问题