R语言之矩阵基本操作

1. 创建矩阵

a <- matrix(1:16,nrow = 4,ncol = 4) #创建一个4行4列的矩阵
b <- matrix(1:12,nrow = 4,ncol = 3) #创建一个4行3列的矩阵
d <- diag(a) #矩阵的对角线元素
c <- diag((diag(d))) #利用对角线元素创建对角矩阵

c=a
c+a #矩阵相加
c-a #矩阵相减
c%*%a #矩阵相乘

2.矩阵信息

dim(a) #矩阵的维度
ncol(a) #矩阵的列数
nrow(a) #矩阵的行数

3. 基本操作

rowSums(a) #矩阵按行求和
rowMeans(a) #矩阵按行求均值
colSums(a) #矩阵按列求求和
colMeans(a) #矩阵按列求均值

4. apply（X, MARGIN, FUN ,…）函数.X:矩阵； MARGIN: 行列 FUN:函数

apply(a, 1, sum) #矩阵按行求和

apply(a, 2, sum) #矩阵按列求和

apply(a, 1,mean) #矩阵按行求均值
b <- matrix(rnorm(100),20,5)

apply(b, 2,var) #矩阵按列求方差
c <- matrix(1:12,3,4)

apply(c, 2, function(x,a)x*2,a=2) #每列乘2