*libname edisto 'z:\Projects\Gratis'; *run;
libname edisto 'e:\STAT 540'; run;
data mon; set edisto.mon;
/* Monthly data collected from the Edisto River */
/* for over 50 years  */
/* The variable cycle represents La Nina/El Nino/normal */
/* climate conditions, but is not part of the data   */
/*x1W= (cycle=1)&(month=12 or month=1 or month=2);
x1Sp= (cycle=1)&(month=3 or month=4 or month=5);
x1Su= (cycle=1)&(month=6 or month=7 or month=8);
x1F= (cycle=1)&(month=9 or month=10 or month=11);
x2W= (cycle=2)&(month=12 or month=1 or month=2);
x2Sp= (cycle=2)&(month=3 or month=4 or month=5);
x2Su= (cycle=2)&(month=6 or month=7 or month=8);
x2F= (cycle=2)&(month=9 or month=10 or month=11);
x1=(cycle=1);
x2= (cycle=2); */
LogQ=log(flow); run;
/* The data was stored so that all January's appeared */
/* first, then all Februaries, etc.  This caused      */
/* problems in the subsequent ARIMA analysis until    */
/* it was noticed then fixed with PROC SORT           */

proc sort data=mon out=month; by year month; run;