R語言學(xué)習(xí)之矩陣

很多人是在線性代數(shù)課學(xué)的矩陣，當時什么問題都沒有，除了一個問題：學(xué)習(xí)矩陣到底有什么用呢？矩陣是一個集合，它里面可以存放很多對象，比如一個行就是一個對象（或者說記錄），每一個對象又有很多（屬性）列。如果把一組對象~屬性表示成矩陣，我們就能很容易取出每個對象對應(yīng)的某個屬性了，并且根據(jù)線性代數(shù)方法考察兩個對象之間的聯(lián)系（相似性）。矩陣的行列數(shù)我們一般稱作維數(shù)。

對于矩陣而言，我們當然想實現(xiàn)以下操作與功能：

  • 矩陣的加減、乘除運算
  • 矩陣的行列切片
  • 最值的快速獲取
  • 線性代數(shù)運算

好在R語言中的矩陣可以很輕易幫我們實現(xiàn)這些功能，有了這一神兵利器，我們就可以游刃有余地操作應(yīng)該算是數(shù)據(jù)分析的基本單位——矩陣了。

創(chuàng)建矩陣

R中直接調(diào)用函數(shù)matrix()可以快速自定義矩陣，下面一行命令可以快速創(chuàng)建一個4行3列的矩陣：

>a<-matrix(c(1:12),nrow=4,ncol=3,byrow=TRUE)
> a
     [,1] [,2] [,3]
[1,]    1    2    3
[2,]    4    5    6
[3,]    7    8    9
[4,]   10   11   12

這里相當于先創(chuàng)建一個向量，再將其轉(zhuǎn)化為一個4x3的矩陣，bynow=TRUE表示會按照把第一行排滿，接著排第二行。

還有一些其它小技巧，比如有時候我們需要初始化一個矩陣，以便于后面對其進行賦值：

> a1<-matrix(0,3,4)
> a1
     [,1] [,2] [,3] [,4]
[1,]    0    0    0    0
[2,]    0    0    0    0
[3,]    0    0    0    0

矩陣行、列、元素的選?。ㄇ衅?/strong>

取第一行第二列元素

> a[1,2]
[1] 2

取第一行元素，這與Matlab很相像

> a[1,]
[1] 1 2 3

取第一行除了第二個元素之外的元素

> a[1,-2]
[1] 1 3

取第一列元素

> a[,1]
[1]  1  4  7 10

取第一列除了第二個元素之外的元素

> a[-2,1]
[1]  1  7 10

矩陣全部元素

> a[,]
     [,1] [,2] [,3]
[1,]    1    2    3
[2,]    4    5    6
[3,]    7    8    9
[4,]   10   11   12

矩陣的基本運算

這里創(chuàng)建一個新的矩陣b、c1，b與a的維數(shù)相同，c1的列、行數(shù)與a的行、列數(shù)分別相等，便于做實驗。

> b<-matrix(c(13:24),nrow=4,ncol=3,byrow = TRUE)
> b
     [,1] [,2] [,3]
[1,]   13   14   15
[2,]   16   17   18
[3,]   19   20   21
[4,]   22   23   24
> c1<-matrix(c(13:24),nrow=3,ncol=4,byrow = TRUE)
> c1
     [,1] [,2] [,3] [,4]
[1,]   13   14   15   16
[2,]   17   18   19   20
[3,]   21   22   23   24

獲取矩陣維數(shù)

> dim(a)
[1] 4 3

加減法運算

矩陣的加減法運算表示兩個矩陣對應(yīng)元素分別進行加減法運算，返回兩個矩陣對應(yīng)元素分別進行加減法運算的矩陣。當然了，矩陣加減法運算前提是兩個矩陣的維數(shù)必須一樣，否則會報錯。

> a+b
     [,1] [,2] [,3]
[1,]   14   16   18
[2,]   20   22   24
[3,]   26   28   30
[4,]   32   34   36

乘除法運算

> a*b
     [,1] [,2] [,3]
[1,]   13   28   45
[2,]   64   85  108
[3,]  133  160  189
[4,]  220  253  288

矩陣的乘除法運算表示兩個矩陣對應(yīng)元素分別進行乘除法運算，返回兩個矩陣對應(yīng)元素分別進行乘除法運算的矩陣。當然了，矩陣乘除法運算前提是兩個矩陣的維數(shù)必須一樣，否則會報錯。

還有就是矩陣的乘法，要求是前面矩陣的列數(shù)等于后面矩陣的列數(shù)，返回一個左邊矩陣行數(shù)x右邊矩陣列數(shù)的矩陣。

> a%*%c1
     [,1] [,2] [,3] [,4]
[1,]  110  116  122  128
[2,]  263  278  293  308
[3,]  416  440  464  488
[4,]  569  602  635  668

線性代數(shù)運算

R語言提供了很多用于線性代數(shù)運算的函數(shù)，常用的列出如下：

eigen() #求特征值和特征向量

solve() #求逆矩陣

chol() #Choleski分解

svd() #奇異值分解

qr() #QR分解

det() #求行列式

dim() #給出行列數(shù)

t() #矩陣轉(zhuǎn)置

矩陣的拼接

R語言矩陣的拼接主要用到兩個函數(shù)，rbind()、cbind()

按行拼接要求兩個矩陣列數(shù)要相同rbind()

> rbind(a,b)
     [,1] [,2] [,3]
[1,]    1    2    3
[2,]    4    5    6
[3,]    7    8    9
[4,]   10   11   12
[5,]   13   14   15
[6,]   16   17   18
[7,]   19   20   21
[8,]   22   23   24

按列拼接要求兩個矩陣行數(shù)要相同cbind()

> cbind(c1,matrix(c(1:6),nrow = 3,byrow = TRUE))
     [,1] [,2] [,3] [,4] [,5] [,6]
[1,]   13   14   15   16    1    2
[2,]   17   18   19   20    3    4
[3,]   21   22   23   24    5    6

其它函數(shù)的靈活結(jié)合

矩陣相關(guān)計算求法還可以靈活應(yīng)用其它函數(shù)，比如求和函數(shù)sum()，平均值函數(shù)mean()，最值函數(shù)max()等。

> a
     [,1] [,2] [,3]
[1,]    1    2    3
[2,]    4    5    6
[3,]    7    8    9
[4,]   10   11   12
> max(a)
[1] 12

第一行最大值

>max(a[1,])
[1] 3
> max(a[,1])
[1] 10

對第一行求和

> sum(a[1,])
[1] 6
> mean(a[,1])
[1] 5.5

還有就是結(jié)合apply()函數(shù)，后面會講到。

用法，舉個例子。apply(Matrix,1,FUN=mean)，這里，F(xiàn)UN=mean計算矩陣Matrix每一行的平均值，以向量的形式返回,中間的參數(shù)‘1’表示求每一行均值，如果是‘2’，表示求每一列均值

比如：

求a每一行平均值

> apply(a,1,mean)
[1]  2  5  8 11

對a每一列分別求和

> apply(a,2,sum)
[1] 22 26 30

好了，關(guān)于矩陣就講到這里，希望對你們有用。下次會講數(shù)組哦，盡情期待！





作者：柯廣的網(wǎng)絡(luò)日志

微信公眾號：Java大數(shù)據(jù)與數(shù)據(jù)倉庫