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上查看代码片派生到我的代码片
</pre><pre code_snippet_id="1943999" snippet_file_name="blog_20161023_4_6616005" name="code" class="plain">#求矩阵的各列平均值
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
至于其他函数,大家可以自己去试试练习一下。如果有什么错误,还请大家留言批评指出
数据分析咨询请扫描二维码