R语言-matrix生成矩阵

主要介绍一下利用matrix函数和rep生成矩阵

在R语言中可以使用matrix()函数来创建矩阵，其语法格式如下：

matrix(data=NA, nrow = 1, ncol = 1, byrow = FALSE, dimnames = NULL)

参数含义如下：

data：矩阵的元素，默认为NA，即未给出元素值的话，各项为NA

nrow：矩阵的行数，默认为1；

ncol：矩阵的列数，默认为1；

byrow：元素是否按行填充，默认按列；

dimnames：以字符型向量表示的行名及列名。

 

rep函数是R语言中重复操作的函数

rep(x,times,each,length.out)

x：代表的是你要进行复制的对象，可以是一个向量或者是一个因子。

times：代表的是复制的次数，只能为正数。负数以及NA值都会为错误值。复制是指的是对整个向量进行复制。

each：代表的是对向量中的每个元素进行复制的次数。

length.out：代表的是最终输出向量的长度。

 

下面我们利用这两个函数生成特定的矩阵

> matrix(rep(1:4,times = 2),nrow = 4 , ncol = 2 ,byrow =T)
     [,1] [,2]
[1,]    1    2
[2,]    3    4
[3,]    1    2
[4,]    3    4
> matrix(rep(1:4,times = 2),nrow = 4 , ncol = 2 ,byrow =F)
     [,1] [,2]
[1,]    1    1
[2,]    2    2
[3,]    3    3
[4,]    4    4

从上面两个执行结果，可以看出，rep生成的是一行向量，matrix把一行的向量按照byrow蚕食指定的填充方向填充

下面我们生成一个4x4的矩阵，要求矩阵的每个位置上的元素等于1/(i+j-1),体会rep两个参数的含义

> I <- matrix(rep(1:4,times = 4),nrow = 4 , byrow = F)
> J <- matrix(rep(1:4,each = 4),nrow = 4 , byrow = T)#实际上J是I的转置
> A <- 1/(I+J-1)
> A
          [,1]      [,2]      [,3]      [,4]
[1,] 1.0000000 1.0000000 1.0000000 1.0000000
[2,] 0.3333333 0.3333333 0.3333333 0.3333333
[3,] 0.2000000 0.2000000 0.2000000 0.2000000
[4,] 0.1428571 0.1428571 0.1428571 0.1428571