我有一个方差协方差矩阵 S:
我正在尝试找到它的倒数。
我拥有的代码是:
然而,如果我使用solve(),我会得到这个:
为什么invS是不正确的?我应该做哪些修改来纠正它?
> S
[,1] [,2]
[1,] 4 -3
[2,] -3 9
我正在尝试找到它的倒数。
我拥有的代码是:
>invS <- (1/((S[1,1]*S[2,2])-(S[1,2]*S[2,1])))*S
[,1] [,2]
[1,] 0.1481481 -0.1111111
[2,] -0.1111111 0.3333333
然而,如果我使用solve(),我会得到这个:
>invSalt <- solve(S)
[,1] [,2]
[1,] 0.3333333 0.1111111
[2,] 0.1111111 0.1481481
为什么invS是不正确的?我应该做哪些修改来纠正它?