不懂矩阵代数（Matrix algebra）到底能不能做统计? 附R代码!

文章预览

许多高阶的统计模型往往都需要理解和掌握矩阵代数，这篇文章介绍几个用于运算矩阵的 常 见 R代码，也许对有些感兴趣的读者会有帮助。  首先，用于后续的运算和演示，创建一个简单的矩阵，命名为A，代码如下：  A 5, 10 ), ncol = 1 ) A 上面的 matrix() 函数创建了一个2*1的矩阵A，如上图所示。 下一步，创建一个对角矩 阵（Diagonal matrix）， 代码如下： B 5, 10 )) B B为一个对角矩阵，即数值分布在对角线上，且其他的格子为0。 下一步，创建一 个单位矩阵(Identity matrix)，代 码如下： I 1, 1 , 1 )) I I为一个单位矩阵，对角线上的值为1，其它格子为0。 可以将A矩阵进行转置，代码如下： t(A) 通过对比原始的A以及转置后的A，可以了解到 t() 函数将A的行和列进行了对调。 那如何进行矩阵的相加呢？代码如下： A + A + A 同理，也可以进行相减：  A - A 如何将A ………………………………