R 中"Solve a System of Equations"的简单解释



我是一名政治学学生,刚刚开始与R合作。除其他外,我使用此程序进行贝叶斯统计。就这一点而言,我有一个问题:

有人能解释一下(有点简单)命令solve(t(x)%*%x)的作用吗?

给定示例:

s<-seq(1,5)
x<-cbind(1,s)

求解输出(t(x)%*%x):

          s
   1.1 -0.3
s -0.3  0.1

我知道t(x)%*%x的作用,但我不知道solve命令的作用。这可能很简单,但我的数学和编程知识非常有限。

提前感谢!

根据?解决:这个通用函数求解x的方程a%*%x=b,其中b可以是向量或矩阵并且:"如果缺失,b被取为单位矩阵,求解将返回a的逆。"

在你的情况下,b不见了。

A <- t(x)%*%x

所以A是

      s
   5 15
s 15 55

现在

B <- solve(A)

给你A的倒数。试试看:

A %*% B给出

    s
  1 0
s 0 1

这就是身份。

最新更新