here's the code:
/****************************************************************
* NOW CODE A MACRO TO REGRESS CO2 FOR 60 SECOND LINEAR FIT WITH *
* A MOVING STARTING TIME EVERY 15 SECONDS *
****************************************************************/
/* DIVERT LOG TO A FILE FOR THE BATCH PROCESSING- ELSE IT HALTS THE PROGRAMME WHEN IT GETS FULL */
proc printto new log="&OUTFOLDER.\CRICHTONLINEARlog.LOG" print="&OUTFOLDER.\CRICHTONLINEAR.LOG";
run;
%MACRO WINDOWS;
%DO h= 01 %TO 21;
%DO i= 20 %TO 170 %BY 15;
DATA START&h._&i.;
SET DATASET.ONE;
IF ETIM < &i. THEN DELETE;
IF ETIM > &i. + 60 THEN DELETE;
IF LICORFILECODE NE &h. THEN DELETE;
RUN;
PROC REG DATA= START&h._&i. OUTEST= MODEL_CO2&h._&i. ;
MODEL CO2_uMOL= ETIM;
ODS OUTPUT PARAMETERESTIMATES= LINFLUX_PARAM&h._&i. FITSTATISTICS= LINFLUX_FIT&h._&i.;
BY LICORFILE OBS;
RUN;
DATA FLUX_CO2&h._&i.;
SET MODEL_CO2&h._&i.;
LINFLUX= ETIM/ 0.12;
DEADBAND= &i.;
PROC PRINT;
RUN;
DATA LINFLUX_FIT2_&h._&i.;
SET LINFLUX_FIT&h._&i. (RENAME=(cValue2= LINFLUX_R2));
IF Label1 ^= 'Dependent Mean' THEN DELETE;
DEADBAND= &i.;
RUN;
DATA LINFLUX_CV&h._&i.;
SET LINFLUX_FIT&h._&i. (RENAME=(cValue1= LINFLUX_CV));
IF MISSING(Label2)= 0 THEN DELETE;
DEADBAND= &i.;
RUN;
DATA LINFLUX_PARAMS&h._&i.;
SET LINFLUX_PARAM&h._&i. (RENAME=(tValue= LINFLUX_T Probt= LINFLUX_P));
IF Variable = 'Intercept' THEN DELETE;
DEADBAND= &i.;
RUN;
DATA LINEARFIT&h._&i.;
MERGE FLUX_CO2&h._&i. LINFLUX_FIT2_&h._&i. LINFLUX_CV&h._&i. LINFLUX_PARAMS&h._&i.;
DROP _MODEL_
_TYPE_
_DEPVAR_
_RMSE_
INTERCEPT
ETIM
CO2_uMOL
MODEL
DEPENDENT
LABEL1
CVALUE1
CVALUE2
NVALUE1
LABEL2
NVALUE2
VARIABLE
DF
ESTIMATE
STDERR;
BY LICORFILE OBS;
RUN;
DATA DATASET.LINEARFIT&h._&i.;
SET LINEARFIT&h._&i.;
RUN;
%END;
%END;
%MEND WINDOWS;
%WINDOWS;
RUN;
DATA ALL_LINEARFIT;
SET DATASET.LINEARFIT:;
BY LICORFILE OBS;
RUN;
DATA DATASET.ALL_LINEARFIT;
SET ALL_LINEARFIT;
RUN;
