测试数据

#定义两个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,但计算速度更快

计算过程



R语言CRAN r语言crossprod_r语言

R语言CRAN r语言crossprod_测试数据_02

R语言CRAN r语言crossprod_测试数据_03

R语言CRAN r语言crossprod_测试数据_04


tcrossprod()

#矩阵内积
> tcrossprod(a,b)
     [,1] [,2]
[1,]   26   30
[2,]   38   44
#注意;tcrossprod等价于x %*% t(y) ,但计算速度更快

计算过程



R语言CRAN r语言crossprod_r语言_05

R语言CRAN r语言crossprod_测试数据_06

R语言CRAN r语言crossprod_测试数据_07

R语言CRAN r语言crossprod_r语言_08