GAUSS:Whiteの標準誤差

YをＸに対してＯＬＳしたとき，
err=Y-X*bとすると，
whiteの標準誤差を計算するプロシージャー．

proc whitesde(Y,X,err);
local n, asyV , q , OMEGA,S0;
n=rows(X);OMEGA=zeros(n,n);
for q (1,n,1);
OMEGA[q,q]=(err[q]^2);
endfor;
asyV=INVPD(X'X)*X'*OMEGA*X*INVPD(X'X);
asyV=sqrt(diag(asyV));
retp(asyV);
endp;

真ん中の部分は，もっと簡単に書けるコマンドが，
あるかも，しれません．