矩阵转置
1、使用t完成标准的矩阵转置
> array(h,dim=c(2,5))->mya
> mya
[,1] [,2] [,3] [,4] [,5] > t(mya)
[,1] [,2]
[1,] 1 2
[2,] 3 4
[3,] 5 6
[4,] 7 8
[5,] 9 10
2、使用aperm函数实现矩阵转置
aperm有2个常用的参数
第一个参数是需要转置的矩阵,第二个参数perm指示新矩阵相对于第一个参数矩阵的维度的下标,比如说,将行转换为列,将列转换为行,将行列次序更换,将第一维的元素与第二维的元素互换,perm设为c(2,1),perm中是维度下标,不是矩阵下标。
> array(h,dim=c(2,5))->mya
> mya
[,1] [,2] [,3] [,4] [,5][5,] 9 10
> aperm(mya,perm=c(2,1))->myb
> myb
[,1] [,2]
[1,] 1 2
[2,] 3 4
[3,] 5 6
[4,] 7 8
[5,] 9 10
>
如果将perm设为c(1,2)表示不交换原矩阵的维度,即不做操作
> mya
[,1] [,2] [,3] [,4] [,5]> aperm(mya,perm=c(1,2))->myb
> myb
[,1] [,2] [,3] [,4] [,5]
[1,] 1 3 5 7 9
[2,] 2 4 6 8 10
>
我们再来看一个3维数组
> array(mya,c(2,2,5))->mya1
> mya1
, , 1
[,1] [,2]
[1,] 1 3
[2,] 2 4
, , 2
[,1] [,2]
[1,] 5 7
[2,] 6 8
, , 3
[,1] [,2]
[1,] 9 1
[2,] 10 2
, , 4
[,1] [,2]
[1,] 3 5
[2,] 4 6
, , 5
[,1] [,2]
[1,] 7 9
[2,] 8 10
> aperm(mya1,perm=c(2,1,3))->myb1
> myb1
, , 1
[,1] [,2]
[1,] 1 2
[2,] 3 4
, , 2
[,1] [,2]
[1,] 5 6
[2,] 7 8
, , 3
[,1] [,2]
[1,] 9 10
[2,] 1 2
, , 4
[,1] [,2]
[1,] 3 4
[2,] 5 6
, , 5
[,1] [,2]
[1,] 7 8
[2,] 9 10
>
> aperm(mya1,perm=c(1,3,2))->myb1
> myb1
, , 1
[,1] [,2] [,3] [,4] [,5]
[1,] 1 5 9 3 7
[2,] 2 6 10 4 8
, , 2
[,1] [,2] [,3] [,4] [,5]
[1,] 3 7 1 5 9
[2,] 4 8 2 6 10
矩阵的维数与行列数
> ncol(mya)
[1] 5
> nrow(mya)
[1] 2
> dim(mya)
[1] 2 5
>
8)矩阵乘积
若A为m×n矩阵,B为n×r矩阵,则他们的乘积AB(有时记做A· B)会是一个m×r矩阵,但前提是m与n相同时才有定义。
> a
[,1] [,2] [,3] [,4] [,5]
[1,] 1 3 5 7 9
[2,] 2 4 6 8 10
> b
[,1] [,2]
[1,] 1 6
[2,] 2 7
[3,] 3 8
[4,] 4 9
[5,] 5 10
> a %*% b
[,1] [,2]
[1,] 95 220
[2,] 110 260
>
>