11.1常用的矩陣計算
> m
C1 C2 C3 C4 C5
R1 1 2 3 4 5
R2 6 7 8 9 10
R3 11 12 13 14 15
R4 16 17 18 19 20
> m+1 #各個元素都加1
C1 C2 C3 C4 C5
R1 2 3 4 5 6
R2 7 8 9 10 11
R3 12 13 14 15 16
R4 17 18 19 20 21
> m*2 #各個元素都乘2
C1 C2 C3 C4 C5
R1 2 4 6 8 10
R2 12 14 16 18 20
R3 22 24 26 28 30
R4 32 34 36 38 40
> m+m #對應(yīng)元素相加
C1 C2 C3 C4 C5
R1 2 4 6 8 10
R2 12 14 16 18 20
R3 22 24 26 28 30
R4 32 34 36 38 40
> n <- matrix(1:20,5,4)
> n
[,1] [,2] [,3] [,4]
[1,] 1 6 11 16
[2,] 2 7 12 17
[3,] 3 8 13 18
[4,] 4 9 14 19
[5,] 5 10 15 20
> m+n #m和n的行和列數(shù)不同所以沒法進(jìn)行運(yùn)算
Error in m + n : 非整合陳列
矩陣之間的運(yùn)算需要相同的行或列习柠，如果想進(jìn)行運(yùn)算可以取出單獨(dú)的行或列進(jìn)行運(yùn)算
> t <- m[1,]
> sum(t)
[1] 15
11.2矩陣計算的函數(shù)
| 序號 |
函數(shù) |
作用 |
| 01 |
colSums(m) |
求m矩陣中每列的和 |
| 02 |
rowSums(m) |
求m矩陣中每行的和 |
| 03 |
colMeans(m) |
求m矩陣中每列的平均值 |
| 04 |
rowMeans(m) |
求m矩陣中每行的平均值 |
> colSums(m)
C1 C2 C3 C4 C5
34 38 42 46 50
> rowSums(m)
R1 R2 R3 R4
15 40 65 90
> colMeans(m)
C1 C2 C3 C4 C5
8.5 9.5 10.5 11.5 12.5
> rowMeans(m)
R1 R2 R3 R4
3 8 13 18
11.3矩陣的乘積（內(nèi)積&外積）
矩陣的內(nèi)積
> n <- matrix(1:9,3,3)
> t <- matrix(2:10,3,3)
> n
[,1] [,2] [,3]
[1,] 1 4 7
[2,] 2 5 8
[3,] 3 6 9
> t
[,1] [,2] [,3]
[1,] 2 5 8
[2,] 3 6 9
[3,] 4 7 10
> n*t
[,1] [,2] [,3]
[1,] 2 20 56
[2,] 6 30 72
[3,] 12 42 90
矩陣的外積
| 序號 |
函數(shù) |
作用 |
| 01 |
diag(m) |
求m階對角矩陣的對角線启涯，返回對角線位置的值 |
| 02 |
t(m) |
將m矩陣進(jìn)行轉(zhuǎn)置扭粱，將行和列進(jìn)行互換日丹，只換數(shù)字不換名稱（row=F) |
> n %*% t
[,1] [,2] [,3]
[1,] 42 78 114
[2,] 51 96 141
[3,] 60 114 168
> diag(n) #返回對角線位置的值
[1] 1 5 9
> n
[,1] [,2] [,3]
[1,] 1 4 7
[2,] 2 5 8
[3,] 3 6 9
> m
C1 C2 C3 C4 C5
R1 1 2 3 4 5
R2 6 7 8 9 10
R3 11 12 13 14 15
R4 16 17 18 19 20
> t(m) # 將m矩陣進(jìn)行轉(zhuǎn)置剿另，將行和列進(jìn)行互換
R1 R2 R3 R4
C1 1 6 11 16
C2 2 7 12 17
C3 3 8 13 18
C4 4 9 14 19
C5 5 10 15 20
