我是一名政治学学生,刚刚开始与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
这就是身份。