我想从44.400个基因的宇宙中创建1000个随机列表中的1652个基因。 我决定更换。我使用下面的指令来创建随机列表:更换的随机列表
randomMatrix<-replicate(1000, sample(gene_list, 1652, replace = T))
重点是在每个列表中复制一个基因。对于我的研究,基因可以在列表之间复制,但不在每个列表中复制。我如何强制不复制每个单子中的基因?
在此先感谢
我想从44.400个基因的宇宙中创建1000个随机列表中的1652个基因。 我决定更换。我使用下面的指令来创建随机列表:更换的随机列表
randomMatrix<-replicate(1000, sample(gene_list, 1652, replace = T))
重点是在每个列表中复制一个基因。对于我的研究,基因可以在列表之间复制,但不在每个列表中复制。我如何强制不复制每个单子中的基因?
在此先感谢
一个reproducible example将是很好的说明你的问题,因为你没有给我们这样的例子,我只是假设一个List
并提出了一些复制
List <- list(c(2,1,3,4,5,6), c(1,4,5,7,0,6), c(2,4,7,9,3,1))
set.seed(001)
replicate(3, lapply(List, sample, 7, replace=TRUE), simplify = FALSE)
产生
[[1]]
[[1]][[1]]
[1] 1 3 4 6 1 6 6
[[1]][[2]]
[1] 7 7 1 4 4 0 5
[[1]][[3]]
[1] 3 7 3 1 7 3 1
[[2]]
[[2]][[1]]
[1] 1 4 2 1 3 2 3
[[2]][[2]]
[1] 6 5 5 7 5 4 0
[[2]][[3]]
[1] 3 3 2 3 7 3 9
[[3]]
[[3]][[1]]
[1] 5 4 4 5 2 3 5
[[3]][[2]]
[1] 0 5 6 5 4 1 1
[[3]][[3]]
[1] 4 9 9 7 1 4 7
请注意,此方法将为您的原始列表中的每个元素提供重新采样数据(包含替换)这就是为什么每个复制都是一个由三个元素组成的列表。
如果你写sapply
而不是lapply
内replicate(...)
输出结果会更好。
set.seed(001)
replicate(3, sapply(List, sample, 7, replace=TRUE), simplify = FALSE)
[[1]]
[,1] [,2] [,3]
[1,] 1 7 3
[2,] 3 7 7
[3,] 4 1 3
[4,] 6 4 1
[5,] 1 4 7
[6,] 6 0 3
[7,] 6 5 1
[[2]]
[,1] [,2] [,3]
[1,] 1 6 3
[2,] 4 5 3
[3,] 2 5 2
[4,] 1 7 3
[5,] 3 5 7
[6,] 2 4 3
[7,] 3 0 9
[[3]]
[,1] [,2] [,3]
[1,] 5 0 4
[2,] 4 5 9
[3,] 4 6 9
[4,] 5 5 7
[5,] 2 4 1
[6,] 3 1 4
[7,] 5 1 7
它应与replace = FALSE
工作:
randomMatrix<-replicate(1000, sample(gene_list, 1652, replace = FALSE))
这当然需要gene_list
至少1652个唯一值。
通过强制我的老板,我必须更换。但我需要的是取代每个列表中没有基因必须复制的条件。当然,列表之间的基因被复制,因为替换= TRUE – Elb
@Elb样本命令被复制1000次。在每次运行中,如果'replace = FALSE',您将获得唯一的值。但这对其他运行没有影响。因此,在列表之间你实际上取代了“基因”。 –
我们喜欢[可重现的例子](http://stackoverflow.com/q/5963269)这里。 – rinni