测试数据
#定义两个2*2的矩阵
a<-matrix(c(1:4),ncol=2)
b<-matrix(c(5:8),ncol=2)
那么,有
> a
[,1] [,2]
[1,] 1 3
[2,] 2 4
> b
[,1] [,2]
[1,] 5 7
[2,] 6 8
进行试验
crossprod()
#矩阵内积
crossprod(a,b)
[,1] [,2]
[1,] 17 23
[2,] 39 53
#注意;crossprod(a,b)等价于t(a)%*%b,但计算速度更快
计算过程
左 | 右 |
tcrossprod()
#矩阵内积
> tcrossprod(a,b)
[,1] [,2]
[1,] 26 30
[2,] 38 44
#注意;tcrossprod等价于x %*% t(y) ,但计算速度更快
计算过程
左 | 右 |