R语言之矩阵操作
R语言作为一门统计语言,肯定缺少不了矩阵处理。下面是我参考《R语言与数据挖掘》总结出来以下关于矩阵运算的函数,其中包括矩阵的求和,转置,等。
接下来我们通过一些例子来了解一下这些函数。
[plain] view plain copy 在CODE上查看代码片派生到我的代码片
#定义两个都是2*2的矩阵
a<-matrix(c(1:4),ncol=2)
b<-matrix(c(5:8),ncol=2)
其中a,b分别为
[plain] view plain copy 在CODE上查看代码片派生到我的代码片
> a
[,1] [,2]
[1,] 1 3
[2,] 2 4
> b
[,1] [,2]
[1,] 5 7
[2,] 6 8
[plain] view plain copy 在CODE上查看代码片派生到我的代码片
#求矩阵的各列平均值 colMeans(a) [1] 1.5 3.5 [plain] view plain copy 在CODE上查看代码片派生到我的代码片 #求矩阵a的转置矩阵 t(a) [,1] [,2] [1,] 1 2 [2,] 3 4 [plain] view plain copy 在CODE上查看代码片派生到我的代码片 #矩阵的四则运算 3*a-b+b/a [,1] [,2] [1,] 3 4.333333 [2,] 3 6.000000 [plain] view plain copy 在CODE上查看代码片派生到我的代码片 #方阵求解行列式 det(a) [1] -2 [plain] view plain copy 在CODE上查看代码片派生到我的代码片 #矩阵内积 crossprod(a,b) [,1] [,2] [1,] 17 23 [2,] 39 53 #注意;crossprod(a,b)等价于t(a)%*%b 至于其他函数,大家可以自己去试试练习一下。如果有什么错误,还请大家留言批评指出
数据分析咨询请扫描二维码