如何在MATLAB(当前为R)中编写此行代码?

时间:2019-01-16 08:15:08

标签: r matlab

如何在MATLAB(当前为R)中编写此代码行?

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

我的尝试是

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

但是我正在为r代码中的'c'做些什么?

2 个答案:

答案 0 :(得分:1)

虽然以上答案解决了matlab代码中缺少的解决方案,但解决方案可能意味着R中有许多不同的事物,

如果方程式中没有逗号,则它不求解任何东西,并且实际上取反数,

MATLAB的A的逆函数:inv(A)R:solve(A)

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

答案 1 :(得分:-1)

c(a,b,c)表示R中的向量。在Matlab中,您应该编写

vec = [a b c];

此外,您需要找到与R-solve()函数等效的函数。到目前为止,您的matlab代码仅使X'与X互变,并且无法求解方程组。

linsolve应该是一个很好的起点。