Как использовать R, чтобы сделать этот матричный расчет?

Question

У меня есть матричное вычисление, и мне нужно использовать R. Для расчета.

Answer 1

Попробуйте:

> z <- matrix(c(1,2,3,4,5,6,7,8,9), 3, byrow=TRUE)
> z
     [,1] [,2] [,3]
[1,]    1    2    3
[2,]    4    5    6
[3,]    7    8    9

> z[1,]
[1] 1 2 3

> matrix(z[1,])
     [,1]
[1,]    1
[2,]    2
[3,]    3

Матричное умножение (вектор-столбец * Вектор-строка):

> matrix(z[1,]) %*% z[1,] + matrix(z[2,]) %*% z[2,] + matrix(z[3,]) %*% z[3,]
     [,1] [,2] [,3]
[1,]   66   78   90
[2,]   78   93  108
[3,]   90  108  126

См .: http://www.philender.com/courses/multivariate/notes/matr.html

Нетот же результат при вычислении Row-Vector * Column-Vector:

> z[1,] %*% matrix(z[1,]) + z[2,] %*% matrix(z[2,]) + z[3,] %*% matrix(z[3,])
     [,1]
[1,]  285

Answer 2

Попробуйте это (с умножением матрицы):

library(optimbase)
Z <- matrix(1:9, byrow = TRUE, nrow = 3)
A <- transpose(Z[1, ]) %*% Z[1, ] + transpose(Z[2, ]) %*% Z[2, ] + transpose(Z[3, ]) %*% Z[3, ]
A # output
     [,1] [,2] [,3]
[1,]   66   78   90
[2,]   78   93  108
[3,]   90  108  126

Answer 3

z <- matrix(1:9,nrow = 3,byrow = TRUE)
z
     [,1] [,2] [,3]
[1,]    1    2    3
[2,]    4    5    6
[3,]    7    8    9

t(z)%*%z
     [,1] [,2] [,3]
[1,]   66   78   90
[2,]   78   93  108
[3,]   90  108  126

Answer 4

Попробуйте это:

z<-matrix(seq(1:9),ncol=3,nrow=3,byrow=TRUE)
A=z[1,]*t(z[1,])+z[2,]*t(z[2,])+z[3,]*t(z[3,])
A
#       [,1] [,2] [,3]
#  [1,]   66   93  126

Edit1: следующий код работает для любого измерения, если это квадратная матрица:

colSums(t(apply(z,1,function(x)(x[1:nrow(z)]*t(x[1:nrow(z)])))))
#[1]  66  93 126

Edit2:

t(z)%*%z #the following 3 lines will all give you the same thing!
crossprod(z) 
matrix(colSums(t(apply(z,1,
function(x)(matrix(x[1:nrow(z)])%*% x[1:nrow(z)])))),ncol=nrow(z),byrow = TRUE)
  #     [,1] [,2] [,3]
 #[1,]   66   78   90
 #[2,]   78   93  108
 #[3,]   90  108  126