r умножение каждой строки матрицы на элементы из столбца - PullRequest
Новая
       17

r умножение каждой строки матрицы на элементы из столбца

1 голос
/ 18 апреля 2020

Если это первые 6 строк моей матрицы nx6 

  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

Цель состоит в том, чтобы создать новую nx5 матрицу на основе этой логики c. Каждая строка из 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

Любые указатели для выполнения этого оценены. Спасибо.

1 Ответ

1 голос
/ 18 апреля 2020

Мы могли бы l oop по последовательности строк с 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")))
...