....
ods listing close;
proc reg data = p010;
model y = x4;
ods output anova=temp;
run;
quit;
ods listing;
data temp;
set temp;
if source = 'Error' then call symput ('sse', ss);
if source = 'Model' then call symput ('p', df );
run;
%put &sse &p;
data resid2;
set resid;
id = _N_;
newid + 1;
d = .;
d = r/sqrt(&sse);
hadi = h/(1-h) + (&p + 1)*d**2/((1-h)*(1-d**2));
keep hadi d id dffits cookd h;
run;
主要是最后的sse和p找不到,会警告说:
WARNING: 没有解析符号引用 SSE。
WARNING: 没有解析符号引用 P。急求高人帮忙啊