Gretl - 如何计算矩阵

Question

我有一个线性回归模型：

yi = α + βxi + ui

我想计算：

(\sigma_u)^2(X'X)^(-1)

我可以在 gretl 中做到这一点吗？怎么做？ 如果不是，如何从 gretl 中取出 X 矩阵？

非常感谢您的回复！！！

Answer 1

这是一种方法。如果您想了解更多信息，我强烈建议您阅读gretl guide 的第 13 章。

open galton
list xlist = const parent child # or list xlist 0 1 2 
matrix X = {xlist}
ols const child parent --quiet
scalar sigma_u = $sigma
matrix res = sigma_u^2 * inv(X'X)
res

##  res
## res (3 x 3)

##    0.0010799  -1.2533e-05  -3.2771e-06 
##  -1.2533e-05   2.7165e-07  -8.8464e-08 
##  -3.2771e-06  -8.8464e-08   1.3688e-07 

Answer 2

dickoa：我想你的回答是：

open galton
list xlist = const parent # or list xlist 0 1 2 
matrix X = {xlist}
ols child const parent --quiet
scalar sigma_u = $sigma
matrix res = sigma_u^2 * inv(X'X)
res

无论如何，这是一种更简单的方法：

open galton
list xlist = const parent # or list xlist 0 1 2 
matrix X = {xlist}
ols child const parent --quiet
matrix res = $vcv
res