gauss怎么导出回归残差？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

具体程序如下：
y = {
2,


3,


1,


7,


5 };


  

x = {
1 3 2,


2 3 1,


7 1 7,


5 3 1,


3 5 5 };

cls;

output file = ols.out reset;
call ols(0,y,x);
output off;
想把残差单独导出用re表示，以便于以后的应用，谢谢。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：GAUSS USS Aus output outpu 残差 GAUSS

output file = 文件夹/ols.out reset;

版主好，可能我没有说清楚。我想得到e=y-b*x这样一个残差序列。不知道怎么实现？

3# tianjianqiang

{ vnam,m,b,stb,vc,stderr,sigma,cx,rsq,resid,dwstat }
= ols(dataset,depvar,indvars);

其中的resid

iooo好，但怎么赋值给e,使之成为一个序列呢？

5# tianjianqiang


//1, 直接
xxi=inv(x'*x);
b=xxi*(x'*y);
e=y-x*b;

//2,或者
{ vnam,m,b,stb,vc,stderr,sigma,cx,rsq,resid,dwstat }
= ols(dataset,depvar,indvars);
//之后
e=reisid;

//3,或者直接
{ vnam,m,b,stb,vc,stderr,sigma,cx,rsq,e,dwstat }
= ols(dataset,depvar,indvars);

谢谢iooo。我试了一下，方法一不包含常数项。虽然方法2和3的resid中的值为零（不知道为什么），但我利用b中的值解决了问题。刚开始学习，以后请多指教。

不需要截距的b=y/x;e=y-x*b;
需要截距的，xx=x~ones(n,1);b=y/xx;e=y-xx*b;

qingqing 的方法也不错，看来还是要多讨论啊。

load data[100,5] = exer.txt; @导入数据@
t = rows(data); @数据的行@
k = 5; @解释变量个数@
y = data[.,1];  x2 = data[.,2];x3 = data[.,3];x4 = data[.,4];x5 = data[.,5];
yy = y; vny = {"y"};
xx = ones(t,1)~x2~x3~x4~x5;  @解释变量@
vnx = {"cons", "x2", "x3", "x4", "x5" };@输出符号变量@

b = invpd(xx'xx)*(xx'yy); @回归系数估计@
e = yy - xx*b ; @回归残差@
s2 = (e'e)/(t-k); @ e'e是残差平方和，t-k是自由度，s2是回归误差的方差@
vb = s2*invpd(xx'xx); @系数的协方差矩阵@
se = sqrt(diag(vb)); @回归估计量的标准差@
tst = b./se; @ 回归估计量显著性的t统计量@
p = 2*cdftc(abs(tst),t-k); @ t统计量p值@
econ = b~se~tst~p;
econ = vnx~econ;

tss = yy'yy - t*meanc(yy)^2; @总平方和@
rss = e'e; @残差平方和@
ess = tss-rss; @回归平方和@
r2 = 1 - rss/tss; @可决系数@
ar2 = 1-(rss/(t-k))/(tss/(t-1)); @调整的可决系数@
f=(ess/(k-1))/(rss/(t-k)); @方程整体显著性F检验@
d = (e[1:(t-1),1] - e[2:t,1])'(e[1:(t-1),1] - e[2:t,1])/(e'e); @DW统计量@