如何在MATLAB（當前為R）中編寫此行代碼？

Question

如何在MATLAB（當前為R）中編寫此行代碼？

vcov_beta_hat <- c(sigma2_hat) * solve(t(X) %*% X)

我的嘗試是

vcov_beta_hat = [sigma2_hat.*((X'*X))];

但是我在r代碼中'c'在做什么？

Answer 1

上面的答案指出解決方案是您的matlab代碼中缺少的內容，而解決方案可能意味着R中有許多不同的事物，

如果方程式中沒有逗號，那么它什么也沒求解，實際上是在求逆，

A的反函數，MATLAB：inv（A）R：solve（A）

因此， vcov_beta_hat = [sigma2_hat.*inv((X'*X))];

Answer 2

c(a,b,c)表示R中的向量。在Matlab中，您將編寫

vec = [a b c];

另外，您需要找到與R- solve()函數等效的函數。 到目前為止，您的matlab代碼僅使X'與X互變，並且不能求解方程組。

linsolve應該是一個很好的起點。