1.矩阵生成
matrix()函数
as.matrix() 将对象强制转换为矩阵
is.matrix()判断对象是否为矩阵
matrix(data = NA, nrow = 1, ncol = 1, byrow = FALSE,dimnames = NULL)
nrow/nr为矩阵的行数，ncol/nc为矩阵的列数，默认值为1,
byrow,TRUE表示将数据按行放置,默认为FALSE，按列放置
dimnames,行和列的名称
> mdat <- matrix(c(1,2,3, 11,12,13), nrow = 2, ncol = 3, byrow = TRUE,
dimnames = list(c("row1", "row2"),
c("C.1", "C.2", "C.3")))
> mdat
C.1 C.2 C.3
row1 1 2 3
row2 11 12 13
> matrix(1:10,2,5)
[,1] [,2] [,3] [,4] [,5]
[1,] 1 3 5 7 9
[2,] 2 4 6 8 10
#生成空矩阵
> A<-matrix(nr=2,nc=5) #2行5列
> #赋值
> A[1,1]<-1
> A[1,5]<-5
> A[2,1]<-6
> A[2,5]<-10
> A
[,1] [,2] [,3] [,4] [,5]
[1,] 1 NA NA NA 5
[2,] 6 NA NA NA 10
#dim获取对象的维度
> dim(A)
[1] 2 5
#nrow 获取行数;ncol获取列数
> nrow(A)
[1] 2
> ncol(A)
[1] 5
2.矩阵运算
(1)矩阵的合并
rbind()按行合并（上下拼接）,需要相同的列数
> A<-as.matrix(c(1:5))
> B<-as.matrix(c(100:104))
> A
[,1]
[1,] 1
[2,] 2
[3,] 3
[4,] 4
[5,] 5
> B
[,1]
[1,] 100
[2,] 101
[3,] 102
[4,] 103
[5,] 104
> rbind(A,B)
[,1]
[1,] 1
[2,] 2
[3,] 3
[4,] 4
[5,] 5
[6,] 100
[7,] 101
[8,] 102
[9,] 103
[10,] 104
cbind() 按列合并（左右拼接）,需要相同的行数
> cbind(A,B)
[,1] [,2]
[1,] 1 100
[2,] 2 101
[3,] 3 102
[4,] 4 103
[5,] 5 104
(2)矩阵的拉直，as.vector()，将矩阵按列拉直
> dat<-cbind(A,B)
> as.vector(dat)
[1] 1 2 3 4 5 100 101 102 103 104
(3)矩阵行列的命名
dim()获得矩阵维度
dimnames()获得矩阵行列的名称
rownames() 获得矩阵行名，或命名行名
colnames()获得矩阵列名，或命名列名
> dat
[,1] [,2]
[1,] 1 100
[2,] 2 101
[3,] 3 102
[4,] 4 103
[5,] 5 104
> dim(dat)
[1] 5 2
> dimnames(dat) #此时无行名和列名
NULL
> rownames(dat)<-paste0("row.",1:5) ##命名行
> colnames(dat)<-paste0("col.",1:2) ##命名列
> dat
col.1 col.2
row.1 1 100
row.2 2 101
row.3 3 102
row.4 4 103
row.5 5 104
> rownames(dat)
[1] "row.1" "row.2" "row.3" "row.4" "row.5"
> colnames(dat)
[1] "col.1" "col.2"
> dimnames(dat)
[[1]]
[1] "row.1" "row.2" "row.3" "row.4" "row.5"
[[2]]
[1] "col.1" "col.2"
（4）矩阵下标
要访问矩阵某个元素或为该元素赋值的时候，只要写出矩阵名和方括号，中间有逗号分开两个下标
> dat
col.1 col.2
row.1 1 100
row.2 2 101
row.3 3 102
row.4 4 103
row.5 5 104
> dat[3,2] #第3行，第2列的值
[1] 102
> dat[3,2]<-500 #第3行，第2列的值改为500
> dat
col.1 col.2
row.1 1 100
row.2 2 101
row.3 3 500
row.4 4 103
row.5 5 104
> dat[1:3,1] ##第1:3行，第1列
row.1 row.2 row.3
1 2 3
若是想要去掉矩阵对应行和/或列，用负号即可
> dat2<-dat[-1,] ##去除矩阵的第一行
> dat2
col.1 col.2
row.2 2 101
row.3 3 500
row.4 4 103
row.5 5 104
