很多人是在线性代数课学的矩阵,当时什么问题都没有,除了一个问题:学习矩阵到底有什么用呢?矩阵是一个集合,它里面可以存放很多对象,比如一个行就是一个对象(或者说记录),每一个对象又有很多(属性)列。如果把一组对象~属性表示成矩阵,我们就能很容易取出每个对象对应的某个属性了,并且根据线性代数方法考察两个对象之间的联系(相似性)。矩阵的行列数我们一般称作维数。
对于矩阵而言,我们当然想实现以下操作与功能:
矩阵的加减、乘除运算
矩阵的行列切片
最值的快速获取
线性代数运算
好在R语言中的矩阵可以很轻易帮我们实现这些功能,有了这一神兵利器,我们就可以游刃有余地操作应该算是数据分析的基本单位——矩阵了。
创建矩阵
R中直接调用函数matrix()可以快速自定义矩阵,下面一行命令可以快速创建一个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
这里相当于先创建一个向量,再将其转化为一个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
矩阵行、列、元素的选取(切片)
取第一行第二列元素
a[1,2]
[1] 2
取第一行元素,这与Matlab很相像