有时,我们需要获得矩阵的逆矩阵
solve()可用来获得逆矩阵,在其文档中是如此介绍的:
solve(a, b, ...)是用来求解a %*% x = b的,但是当参数b省略时,b会被设为单位矩阵,此时solve()返回a的逆
如果a是不可逆的,那么solve()将会报错,提示线性系统正好是奇异的
当a不可逆时,可以使用MASS包中的ginv()获得广义逆,当a可逆时,ginv()返回的逆矩阵与solve()的结果是一致的
比如下面的例子:
有时,我们需要获得矩阵的逆矩阵
solve()可用来获得逆矩阵,在其文档中是如此介绍的:
solve(a, b, ...)是用来求解a %*% x = b的,但是当参数b省略时,b会被设为单位矩阵,此时solve()返回a的逆
如果a是不可逆的,那么solve()将会报错,提示线性系统正好是奇异的
当a不可逆时,可以使用MASS包中的ginv()获得广义逆,当a可逆时,ginv()返回的逆矩阵与solve()的结果是一致的
比如下面的例子: