2012-10-05 69 views
-1

我想从44.400个基因的宇宙中创建1000个随机列表中的1652个基因。 我决定更换。我使用下面的指令来创建随机列表:更换的随机列表

randomMatrix<-replicate(1000, sample(gene_list, 1652, replace = T)) 

重点是在每个列表中复制一个基因。对于我的研究,基因可以在列表之间复制,但不在每个列表中复制。我如何强制不复制每个单子中的基因?

在此先感谢

+1

我们喜欢[可重现的例子](http://stackoverflow.com/q/5963269)这里。 – rinni

回答

0

一个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而不是lapplyreplicate(...)输出结果会更好。

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 
+0

嗨吉贝特。非常感谢你的帮助,但可能我不明白你的建议。我想避免的是为前例。在[[1]]中重复[,1]中的数字6,在[[2]中重复7和4等等。 – Elb

+0

@Elb如此设置'replace = FALSE'。它不会像你期望的那样工作,然后给我们一些数据来处理并向我们展示预期的输出,以便我们可以完全理解你的问题。 –

+0

mmmmm ...我不能....我的老板告诉我根据我们的数据进行替换样品.. – Elb

2

它应与replace = FALSE工作:

randomMatrix<-replicate(1000, sample(gene_list, 1652, replace = FALSE)) 

这当然需要gene_list至少1652个唯一值。

+0

通过强制我的老板,我必须更换。但我需要的是取代每个列表中没有基因必须复制的条件。当然,列表之间的基因被复制,因为替换= TRUE – Elb

+0

@Elb样本命令被复制1000次。在每次运行中,如果'replace = FALSE',您将获得唯一的值。但这对其他运行没有影响。因此,在列表之间你实际上取代了“基因”。 –