R中的随机元素乘法

时间:2016-06-23 14:34:37

标签: r matrix random elementwise-operations

我一直在网站上搜索我的问题的答案,我是R的新人所以我希望这甚至是可能的。我有两个大型模拟矩阵(A = 100,000 x 50和B = 10,000 x 50),我想逐行地逐行乘法。

基本上我希望A中的每一行从B中随机选择一行进行逐元素乘法。

A:

      [,1] [,2] [,3] [,4] [,5]
[1,]    1    1    1    1    1
[2,]    1    1    1    1    1
[3,]    1    1    1    1    1
[4,]    1    1    1    1    1
[5,]    1    1    1    1    1
[6,]    1    1    1    1    1
[7,]    1    1    1    1    1
[8,]    1    1    1    1    1
[9,]    1    1    1    1    1
[10,]   1    1    1    1    1

和 B:

      [,1] [,2] [,3] [,4] [,5]
[1,]    1    1    1    1    1
[2,]    2    2    2    2    2
[3,]    3    3    3    3    3
[4,]    4    4    4    4    4
[5,]    5    5    5    5    5

是否有一个运算符可以通过A的行并从B中随机选择一行来进行元素乘法运算?对于这样的结果:

 C <- A&*&B
 C
 A[1,]*B[3,]
 A[2,]*B[1,]
 A[3,]*B[2,]  
 A[4,]*B[5,] 
 A[5,]*B[3,]  
 A[6,]*B[4,] 
 A[7,]*B[1,]  
 A[8,]*B[5,] 
 A[9,]*B[2,] 
 A[10,]*B[2,]

谢谢!

1 个答案:

答案 0 :(得分:5)

试试这个:

row_id <- sample(1:nrow(B), nrow(A), replace = TRUE)
A * B[row_id, ]

我想我只需要解释sample()的作用。考虑:

sample(1:5, 10, replace = TRUE)
[1] 4 5 2 4 1 2 2 1 2 5

我没有按set.seed()设置随机种子,所以当你运行它时,你会得到不同的结果。但你需要知道的是:它是随机的。