post-7-plot-run-C-mcmc-a.R

## Item discrimination 
check.sampler.graph(
    a.disc,
    all.eap[ get.2pl.params(NULL,1:I.items,NULL,NULL)],
    desc="Item discrimination", ylab="EAP Estimates", col="blue" )