While traditional statistics courses teach students to calculate intervals and test for binomial proportions using a normal or t approximation, this method does not always work well. Agresti and Coull ("Approximate is better than "exact' for interval ...

Read more »

Tags: %sysevalf, binom.test(), Clopper-Pearson CI, coverage probabilities, mode=include, pdf function, Plus 4 estimator, PROC GPLOT, prop.test(), qt(), quantile function, Wilson estimator

Posted in SAS | Comments Off on Example 9.37: (Mis)behavior of binomial confidence intervals

A colleague recently asked "why should the average get closer to the mean when we increase the sample size?" We should interpret this question as asking why the standard error of the mean gets smaller as n increases. The central limit theorem shows t...

Read more »

Tags: a*b=c syntax, apply(), graphics, looping, PROC GPLOT, running average, t() function, type="n"

Posted in SAS | Comments Off on Example 9.19: Demonstrating the central limit theorem

A facebook friend posted the picture reproduced above-- it makes the case that President Obama has been a successful creator of jobs, and also paints GW Bush as a president who lost jobs. Another friend pointed out that to be fair, all of Bush's presi...

Read more »

Tags: abline(), axis statement, Bureau of Labor Statistics, col option, date formats, href option, job creation, mtext(), offset option, plot.ts(), PROC GPLOT, reflabel option, time series, ts()

Posted in SAS | Comments Off on Example 9.11: Employment plot

A facebook friend posted the picture reproduced above-- it makes the case that President Obama has been a successful creator of jobs, and also paints GW Bush as a president who lost jobs. Another friend pointed out that to be fair, all of Bush's presi...

Read more »

Tags: abline(), axis statement, Bureau of Labor Statistics, col option, date formats, href option, job creation, mtext(), offset option, plot.ts(), PROC GPLOT, reflabel option, time series, ts()

Posted in SAS | Comments Off on Example 9.11: Employment plot

A facebook friend posted the picture reproduced above-- it makes the case that President Obama has been a successful creator of jobs, and also paints GW Bush as a president who lost jobs. Another friend pointed out that to be fair, all of Bush's presi...

Read more »

Tags: abline(), axis statement, Bureau of Labor Statistics, col option, date formats, href option, job creation, mtext(), offset option, plot.ts(), PROC GPLOT, reflabel option, time series, ts()

Posted in SAS | Comments Off on Example 9.11: Employment plot

The book "Bayesian Computation with R" by Jim Albert is an easy to read entry level book on applied Bayesian Statistics. While the book was written for R users, it is not difficult to translate the languages between R and SAS and I believe it is a go...

Read more »

Tags: Bayesian, Graph, PROC GPLOT, PROC MEANS, PROC REG

Posted in SAS | Comments Off on Bayesian Computation with SAS (1).