'Open data
wfopen .\card_data

'Figure-1
freeze(mode=overwrite,figure1) cards.line
figure1.axis dual
figure1.axis(l) units(b)
figure1.axis(r) units(b)

'Figure-2
series change = dlog(cards)*100
freeze(mode=overwrite,figure2) change.line
figure2.axis dual

'Generate seasonal dummies (excluding January)
group seasons.add @expand(@month,@dropfirst)

'Run the models
cards.locallinear(log,save,output) 'stochastic trend
cards.locallinear(log,smooth,save,output) 'smooth stochastic trend
cards.locallinear(xser=seasons,log,smooth,save,output) 'smooth stochastic trend with seasonality

'Generate week-on-week percentage changes
series change_sm = dlog(cards_smlevel)*100
series change_sm01 = dlog(cards_smlevel01)*100
series change_sm02 = dlog(cards_smlevel02)*100

'Figure-4
group compare.add change change_sm change_sm01 change_sm02
freeze(mode=overwrite,figure4) compare.line
figure4.axis dual
figure4.axis(l) range(-4,4)
figure4.setelem(1) linecolor(@rgb(194,210,235))
figure4.setelem(2) linepattern(SOLID)
figure4.setelem(3) linepattern(SOLID)
figure4.setelem(4) linepattern(SOLID)
figure4.legend columns(4)
figure4.setelem(1) legend(Actual)
figure4.setelem(2) legend(Stochastic)
figure4.setelem(3) legend(Smooth Stochastic)
figure4.setelem(4) legend(Smooth Stochastic (w. seasonality))
figure4.legend columns(2) position(BOTCENTER)