如何在MATLAB(当前为R)中编写此代码行?
vcov_beta_hat <- c(sigma2_hat) * solve(t(X) %*% X)
我的尝试是
vcov_beta_hat = [sigma2_hat.*((X'*X))];
但是我正在为r代码中的'c'做些什么?
答案 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应该是一个很好的起点。