我们可以用sapply
遍历行的序列,然后在最后一列和其余列之间进行交叉生产,而我们将从该行到最后一列的元素用于最后一列
如果这些是我的nx6矩阵的前6行
x1 x2 x3 x4 x5 a1
1 30 3 43 0 0.20
3 40 2 35 0 0.21
3 50 7 72 0 0.75
4 40 58 63 10 0.55
1 20 19 61 10 0.43
4 20 5 49 0 0.62
目标是基于此逻辑创建new nx5
矩阵。 X1-X5的每一行都应乘以a1中的标量值,例如]
a2 [1,30,3,43,0]*0.20 + [1,30,3,43,0]*0.21 + [1,30,3,43,0]*0.75 + [1,30,3,43,0]*0.55 + [1,30,3,43,0]*0.43 + [1,30,3,43,0]*0.62 [3,40,2,35,0]*0.21 + [3,40,2,35,0]*0.75 + [3,40,2,35,0]*0.55 + [3,40,2,35,0]*0.43 + [3,40,2,35,0]*0.62 [3,50,7,72,0]*0.75 + [3,50,7,72,0]*0.55 + [3,50,7,72,0]*0.43 + [3,50,7,72,0]*0.62 [4,40,58,63,10]*0.55 + [4,40,58,63,10]*0.43 + [4,40,58,63,10]*0.62 [1,20,19,61,10]*0.43 + [1,20,19,61,10]*0.62 [4,20,5,49,0]*0.62
最后得到的是下面的nx5矩阵
2.76 82.8 8.28 118.68 0 7.68 102.4 5.12 89.6 0 7.05 117.5 16.45 169.2 0 6.4 64 92.8 100.8 16 1.05 21 19.95 64.05 10.5 2.48 12.4 3.1 30.38 0
任何用于完成此操作的指标,请参见。谢谢。
如果这些是我的nx6矩阵的前6行x1 x2 x3 x4 x5 a1 1 30 3 43 0 0.20 3 40 2 35 0 0.21 3 50 7 72 0 0.75 4 40 58 63 10 0.55 1 ...]]] >
我们可以用sapply
遍历行的序列,然后在最后一列和其余列之间进行交叉生产,而我们将从该行到最后一列的元素用于最后一列
t( sapply(seq_len(nrow(m1)), function(i) { x1 <- m1[i:nrow(m1),6] x2 <- m1[i, -6] matrix(x2, length(x2), length(x1)) %*% x1 }))
数据
m1 <- structure(c(1, 3, 3, 4, 1, 4, 30, 40, 50, 40, 20, 20, 3, 2, 7,
58, 19, 5, 43, 35, 72, 63, 61, 49, 0, 0, 0, 10, 10, 0, 0.2, 0.21,
0.75, 0.55, 0.43, 0.62), .Dim = c(6L, 6L), .Dimnames = list(NULL,
c("x1", "x2", "x3", "x4", "x5", "a1")))
我们可以用sapply
遍历行的序列,然后在最后一列和其余列之间进行交叉生产,而我们将从该行到最后一列的元素用于最后一列