filename what "~/sasuser.v94/private/hw2p2_discussion_data.txt";

data look;
infile what;
time=_n_;
input z1 z2 z3 z4;
run;

title "Series z2 -- first third";
proc arima data=look(firstobs=1 obs=133);
identify var=z2;
run;

title "Series z2 -- last third";
proc arima data=look(firstobs=267 obs=400);
identify var=z2;
run;

title "Series z1 -- first half";
proc arima data=look(firstobs=1 obs=200);
identify var=z1;
run;

title "Series z1 -- second half";
proc arima data=look(firstobs=201 obs=400);
identify var=z1;
run;

title "Series z1 -- histogram of first half";
proc sgplot data=look(firstobs=1 obs=200);
histogram z1 ;
run ;

title "Series z1 -- histogram of second half";
proc sgplot data=look(firstobs=201 obs=400);
histogram z1 ;
run ;