/*把原始的txt数据文件全部放入 d:\sean\ 中,不要改动,复制运行下面代码即可,结果文件也是在d:\sean\中
http://www.pinggu.org/bbs/thread-705581-1-1.html*/
%let path=d:\sean\;
%let path1=d:\;
x "dir &path.*.txt /b > &path1.fname.txt";
proc import datafile="d:\fname.txt" out=tmp dbms=tab replace;
getnames=no;delimiter=' ';run;
data fname;set tmp;
fname=compress('%a('||var1||')','.txt');keep fname;
file "d:\ex.txt";put fname;run;
%macro a(f);data &f;
infile "d:\sean\&f..txt" firstobs=2;
input date:yymmdd10. lighttrapdata modelprediction;
format date yymmdd10.;run;
ODS HTML file="d:\sean\&f..HTML";
title 'Light-trap-data_my_gnuplot_Prediction vs Lighttrap'&f;
PROC PLOT;
PLOT ModelPrediction*LighttrapData='*'; /* 作原始数据散点图,以了解是否适合作直线相关和回归分析 */
RUN;
PROC CORR;
VAR LighttrapData ModelPrediction; /* 调用CORR过程作相关分析 */
RUN;
PROC REG;
MODEL ModelPrediction=LighttrapData;
run;
PROC REG;
MODEL ModelPrediction=LighttrapData / noint R CLI CLM;
PLOT (U95. L95. P.)*LighttrapData='-' ModelPrediction*LighttrapData='o' / OVERLAY;
OUTPUT OUT=my1980 P=yhat R=yresid;
RUN;
PROC PLOT DATA=my1980;
PLOT yresid*yhat='R'; /* 调用PLOT过程绘残差图,纵轴为残差,横轴为Y的预测值 */
RUN;
ODS HTML close;
%mend a;%include "d:\\ex.txt";run;