刚入手mata , 有下面的算法
mata
: st_view(x=. , . , ("weight","foreign","cons"))
: st_view(y=. , . , "mpg")
: b=invsym(x'*x)*x' *y
: e=y-x*b
: k=cols(x)
: s2=(e:^2)/[rows(x)-k]
: n=rows(x)
: s2=(e' * e)/(n-k)
这里有个问题,为什么我用s2=(e:^2)/[rows(x)-k] 或者是 : s2=(e'*e)/[rows(x)-k] 就要报错
如果想计算robust方差 输入
xe=x:*e
Vr=invsym(x'*x)*(xe'*xe)*invsym(x'*x)rows(x)/[rows(x)-k]
诚心 请教错误原因