小编Lut*_*ett的帖子

如何在Matlab中从矩阵中随机选择x个索引

我正在尝试在零矩阵中高效地生成随机分散但数量有限的1。

假设我有一个10x10的零矩阵(zeros(10)),我想随机放置10个1,所以它看起来像:

     0     0     0     0     0     0     0     0     0     1
     0     0     0     0     0     1     0     0     0     0
     0     0     1     0     0     0     0     0     0     0
     0     0     0     0     0     0     0     0     1     0
     1     0     0     0     0     1     0     0     0     0
     0     0     0     0     0     0     0     0     0     0
     0     0     0     1     0     0     0     0     0     0
     0     1     0     0     0     0 …
Run Code Online (Sandbox Code Playgroud)

random matlab matrix indices

2
推荐指数
1
解决办法
55
查看次数

如何在ggtree中按组为树的尖端着色?

如何根据样本所属的组简单地为树的尖端着色?

nwk <- system.file("extdata", "sample.nwk", package="treeio")
tree <- read.tree(nwk)

group_1 <- rep("1.1.1", 5)
group_2 <- rep("1.1.2", 3)
group_3 <- rep("1.2", 2)
group_4 <- rep("1.2.1", 2)
group_5 <- "1.2"

meta_data <- data.frame(ID = LETTERS[1:13], 
                        group = c(group_1, group_2, group_3, group_4, group_5))

Run Code Online (Sandbox Code Playgroud)

我该如何处理元数据才能使提示显示为彩色圆圈或正方形(或其他形状)并根据这些颜色创建图例?

我在 ggtree 文档中看到了类似的树,但是它们的代码非常复杂(例如,请参见https://guangchuangyu.github.io/ggtree-book/chapter-ggtree.html中的第 4.3.7.3 节)。

另外,我无法重现这些示例,因为它们依赖于没有明确说明可用的元数据:

treefile <- "RAxML_bestTree.Aln_All_H3.nwk"
tipseqfile <- "Aln_All_H3_filted.fas"
Run Code Online (Sandbox Code Playgroud)

例如。

必须有一种简单的方法来使用此元数据为提示着色。

r ggtree

1
推荐指数
1
解决办法
1007
查看次数

标签 统计

ggtree ×1

indices ×1

matlab ×1

matrix ×1

r ×1

random ×1