Question

我是政治学的学生，刚刚开始与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的作用，但我无法弄清楚解决命令的作用。这可能很简单，但我对数学和编程的知识非常有限。

提前致谢！

Answer 1

从？解决的描述： “这个通用函数解决了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

这是身份。