Abstract
Abstract. Single-case experimental designs (SCEDs) are interrupted time-series designs that have recently gained recognition as being able to provide a strong basis for establishing intervention effect. Typically, SCED data are short time series and autocorrelated, which renders maximum likelihood and parametric analyses inadequate for data analysis, respectively. Although Bayesian methods overcome these challenges, most practitioners do not use Bayesian estimation because of: (a) its steep learning curve, (b) lack of Bayesian training, and (c) lack of knowledge of Bayesian software solutions. This study demonstrates two Bayesian interrupted time-series models using freeware programs R and JAGS. Practitioners could modify these codes and run them for their own data by changing the values in the codes where indicated. Providing practitioners with such tools to facilitate their analysis is one way to improve methodological rigor in applied research.
References
2004). Evidence-based practice in communication disorders: An introduction. (Technical report). Retrieved from http://shar.es/11yOzJ or http://www.asha.org/policy/TR2004-00001/
. (2003). Outcomes of an intervention to improve hospital antibiotic prescribing: Interrupted time series with segmented regression analysis. Journal of Antimicrobial Chemotherapy, 52, 842–848. https://doi.org/10.1093/jac/dkg459
(2000). Bayesian factor analysis for multilevel binary observations. Psychometrika, 65, 475–497. https://doi.org/10.1007/BF02296339
(2000). A hierarchical Bayesian approach for modeling heterogeneity in structural equation models. Marketing Science, 19, 328–347. https://doi.org/10.1287/mksc.19.4.328.11789
(1995). Impact of legislation requiring moped and motorbike riders to wear helmets. Evaluation and Program Planning, 18, 311–320.
(1970). Time series analysis: Forecasting and control. San Francisco, CA: Holden-Day.
(1998). Some issues in monitoring convergence of iterative solutions. Journal of Computational and Graphical Statistics, 7, 434–455.
(1968). The Connecticut crackdown on speeding: Time-series data in quasiexperimental analysis. Law & Society Review, 3, 33–54.
(2014). CEC’s standards for classifying the evidence base of practices in special education. Remedial and Special Education, 39, 305–318. https://doi.org/10.1177/0741932514557271
(1979). Quasi-experimentation: Design and analysis issues for field settings. Boston, MA: Houghton Mifflin.
(2016). runjags: An R package providing interface utilities, parallel computing methods and additional distributions for MCMC models in JAGS. Journal of Statistical Software, 71, 1–25. https://doi.org/10.18637/jss.v071.i09
(2010). Bayesian item response modeling: Theory and applications. New York, NY: Springer.
(1997). Designing single case treatment studies for aphasic patients. Neuropsychological Rehabilitation, 7, 401–418.
(2011). N-of-1 trials in the medical literature: A systematic review. Medical Care, 49, 761–768. https://doi.org/10.1097/MLR.0b013e318215d90d
(2013). Bayesian data analysis (3rd ed.). London, UK: Chapman & Hall.
(2014). Bayesian methods: A social and behavioral sciences approach (3rd ed.). Boca Raton, FL: Chapman & Hall.
(2012). A standardized mean difference effect size for single case designs. Research Synthesis Methods, 3, 224–239. https://doi.org/10.1002/jrsm.1052
(2013). A standardized mean difference effect size for multiple baseline designs across studies. Research Synthesis Methods, 4, 324–341. https://doi.org/10.1002/jrsm.1086
(1983). Simulation run length control in the presence of an initial transient. Operations Research, 31, 1109–1144. https://doi.org/10.1287/opre.31.6.1109
(1985). Autocorrelation in applied behavior analysis: A myth. Behavioral Assessment, 7, 107–118.
(1998). Irrelevant autocorrelation in least-squares intervention models. Psychological Methods, 3, 104–116.
(2000). A simple and powerful test for autocorrelated errors in OLS intervention models. Psychological Reports, 87, 3–20. https://doi.org/10.2466/PR0.87.5.3-20
(2012). Recovery of graded response model parameters: A comparison of marginal maximum likelihood and Markov chain Monte Carlo estimation. Applied Psychological Measurement, 36, 399–419.
(2010). Single-case designs technical documentation. Retrieved from http://ies.ed.gov/ncee/wwc/pdf/wwc_scd.pdf
(2013). Bayesian estimation supersedes the t-test. Journal of Experimental Psychology: General, 142, 573–603. https://doi.org/10.1037/a0029146
(2015). Doing Bayesian data analysis: A tutorial with R, JAGS, and Stan. (2nd ed.). Burlington, MA: Academic Press/Elsevier.
(2006). Effects of response cards on disruptive behavior and academic responding during math lessons by fourth-grade urban students. Journal of Positive Behavior Interventions, 8, 88–99. https://doi.org/10.1177/10983007060080020701
(2009). The BUGS project: Evolution, critique, and future directions. Statistics in Medicine, 28, 3049–3067. https://doi.org/10.1002/sim.3680
(2007). Introduction to applied Bayesian statistics and estimation for social scientists. New York, NY: Springer.
(2015). Further developments in summarising and meta-analysing single-case data: An illustration with neurobehavioural interventions in acquired brain injury. Neuropsychological Rehabilitation, 25, 637–662. https://doi.org/10.1080/09602011.2015.1064452
(2014). Three-level analysis of single-case experimental data: Empirical validation. The Journal of Experimental Education, 82, 1–21. https://doi.org/10.1080/00220973.2012.745470
(1992). Alcohol availability and consumption: Iowa sales data revisited. Journal of Studies on Alcohol, 53, 487–494. https://doi.org/10.15288/jsa.1992.53.487
(2015). Comparing interval estimates for small sample ordinal CFA models. Frontiers in Psychology: Quantitative Psychology and Measurement, 6, 1599.
(2019). Conceptualizing statistical simulation studies as multiple methods inquiry: A primer. International Journal of Multiple Research Approaches, 11, 77–89.
(2019 April). Adopting a meta-generative way of thinking in the field of education via the use of Bayesian methods. Paper presented at the Annual Meeting of the American Educational Research Association, Toronto, Canada.
(2017). Bayesian unknown change-point models to investigate immediacy in single case designs. Psychological Methods, 22, 743–759. https://doi.org/10.1037/met0000134
(2019 April). Accurate model vs. accurate estimates: A study of Bayesian single-case experimental designs. Paper presented at the Annual Meeting of the American Educational Research Association, Toronto, Canada.
(2016). Bayesian prior choice in IRT estimation using MCMC and Variational Bayes. Frontiers in Psychology, 7, 1–11. https://doi.org/10.3389/fpsyg.2016.01422
(2003).
(JAGS: A program for analysis of bayesian graphical models using gibbs sampling . In K. HornikF. LeischA. ZeileisEds., Proceedings of the 3rd International Workshop on Distributed Statistical Computing. Vienna, Austria: Technical University.2014). R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. Retrieved from http://www.R-project.org/
. (2014). Nonlinear Bayesian analysis for single case designs. Journal of School Psychology, 52, 179–189. https://doi.org/10.1016/j.jsp.2013.12.003
(2011). Effect of pay for performance on the management and outcomes of hypertension in the United Kingdom: Interrupted time series study. British Medical Journal, 342, 108–114. https://doi.org/10.1136/bmj.d108
(2002). Experimental and quasi-experimental designs for generalized causal inference. Boston, MA: Houghton Mifflin.
(2015). The role of between-case effect size in conducting, interpreting, and summarizing single-case research (NCER 2015-002). Washington, DC: National Center for Education Research, Institute of Education Sciences, US Department of Education. Retrieved from http://ies.ed.gov/
(2013). Bayesian estimates of autocorrelations in single-case designs. Behavior Research Methods, 45, 813–821. https://doi.org/10.3758/s13428-012-0282-1
(2011). Characteristics of single-case designs used to assess treatment effects in 2008. Behavior Research Methods, 43, 971–980. https://doi.org/10.3758/s13428-011-0111-y
(2014). Using generalized additive (mixed) models to analyze single case designs. Journal of School Psychology, 52, 149–178. https://doi.org/10.1016/j.jsp.2013.11.004
(2014). Assisting students with autism to actively perform collaborative walking activity with their peers using dance pads combined with preferred environmental stimulation. Research in Autism Spectrum Disorders, 8, 1591–1596. https://doi.org/10.1016/j.rasd.2014.08.011
(2011). Stan version 2.19. Retrieved from https://mc-stan.org/about/
. (2014). Examining the internal validity and statistical precision of the comparative interrupted time series design by comparison with a randomized experiment. American Journal of Evaluation, 35, 311–327.
(2014). An effect size measure and Bayesian analysis of single-case designs. Journal of School Psychology, 52, 213–230. https://doi.org/10.1016/j.jsp.2013.12.002
(2013). Opioid overdose rates and implementation of overdose education and nasal naloxone distribution in Massachusetts: Interrupted time series analysis. British Medical Journal, 346, 174–185. https://doi.org/10.1136/bmj.f174
(2015). Adding design elements to improve time series designs: No child left behind as an example of causal pattern-matching. Journal of Research on Educational Effectiveness, 8, 245–279. https://doi.org/10.1080/19345747.2013.878011
(