*
* Example BLOOD.DAT from pp 504-508
*
open data blood.dat
data(format=free,org=columns) 1 187 week blood event
*
* Run BJIDENT over the pre-intervention data
*
@bjident(diffs=1) blood 1 134
*
boxjenk(const) blood * 134
@regcorrs
*
boxjenk(const,inputs=1) blood
# event
@regcorrs
*
prj fitted
graph(footer="Figure 12-7 Intervention Model for Blood Demand") 2
# blood
# fitted