%include data(preacher); goptions vsize=6in; title 'Model PCB: Preacher, Curran, Bower 2007'; proc glm data=preacher; model math = anti | hyper age grade female minority / ss2; run; title 'PCB model: Anti * Hyper'; %rug(data=preacher, var=anti, on=X, at=0, ht=1, color=blue, jitter=0.5); axis1 order=(30 to 45 by 5) label=(a=90) minor=(n=1); %effplot(data=preacher, y=math, ylab=Math ability, model=anti | hyper age grade female minority, class=female minority, grid=%str(anti=0 to 6, hyper=0 to 6 by 2, age=60 to 140 by 10, grade=0 to 6, female=0 1, minority=0 1), xvar=anti, cvar=hyper, vaxis=axis1, print=Y, anno=_rug_, proc=glm);