R: матричное умножение вектора
У меня следующая проблема:
myvec <- c(1:3)
mymat <- as.matrix(cbind(a = 6:15, b = 16:25, c= 26:35))
mymat
a b c
[1,] 6 16 26
[2,] 7 17 27
[3,] 8 18 28
[4,] 9 19 29
[5,] 10 20 30
[6,] 11 21 31
[7,] 12 22 32
[8,] 13 23 33
[9,] 14 24 34
[10,] 15 25 35
Я хочу умножить mymat на myvec и построить новый вектор так, чтобы
sum(6*1, 16*2, 26*3)
sum(7*1, 17*2, 27*3)
....................
sum(15*1, 25*2, 35*3)
Извините, это простой вопрос, который я не знаю...
Редактировать: опечатка исправлена
4 ответа
Решение
%*%
Оператор в R выполняет умножение матриц:
> mymat %*% myvec
[,1]
[1,] 116
[2,] 122
...
[10,] 170
Если вы хотите получить результат вектора вместо результата матрицы, вы также можете использовать..
apply(mymat*myvec,1,sum)
Немного больше обходного решения, чем у Оуэна, но оно работает.
Альтернативный, но более длинный путь может быть таким:
rowSums(t(apply(mymat, 1, function(x) myvec*x)),na.rm=T)
Это единственный способ, который я нашел, который может игнорировать NA внутри матрицы.
Матрицы являются векторами в главном порядке столбцов:
colSums( t(mymat) * myvec )
(Отредактировано после правильного прочтения вопроса на этот раз.)