Computes and Plots Compatibility (Confidence) Intervals, P-Values, S-Values, & Likelihood Intervals to Form Consonance, Surprisal, & Likelihood Functions

Allows one to compute compatibility (confidence) intervals for various statistical tests along with their corresponding P-values, S-values, and likelihoods. The intervals can be plotted to create consonance, surprisal, and likelihood functions allowing one to see what effect sizes are compatible with the test model at various compatibility levels rather than being limited to one interval estimate such as 95\%. Functions can also be compared to one another to see how much they overlap with one another and differ. Results can also be exported for Word, Powerpoint, and TeX documents. The package currently supports bootstrapping, linear models, generalized linear models, linear mixed-effects models, survival analysis, and meta-analysis. These methods are discussed by Poole C. (1987) , Schweder T, Hjort NL. (2002) , Singh K, Xie M, Strawderman WE. (2007) , Rothman KJ, Greenland S, Lash TL. (2008, ISBN:9781451190052), Greenland S. (2019) , Chow ZR, Greenland S. (2019) , and Greenland S, Chow ZR. (2019) .

A Single Interval Isn’t Enough

Interval estimates such as confidence compatibility/consonance intervals are now widely reported in many journals alongside the exact P-value of a statistical test and point estimate.

Source: Kristoffer Magnusson

While this is a large improvement over what constituted statistical reporting in the past two decades, it is still largely inadequate.

Take for example, the 95% compatibility interval. As many have stated before, there is nothing special about 95%, yet we rarely see intervals of any other level. Choosing to compute a 95% interval is as mindless as choosing a 5% alpha level for hypothesis testing. A single compatibility interval is only a slice of a wide range of compatibility intervals at different levels. Reporting 95% intervals only promotes cargo-cult statistics since there is not much thought behind the choice. (1)

rather than conscientious practice*.” - Stark & Saltelli, 2018

Thus, we propose that instead of only calculating one interval estimate, every interval associated with a compatibility level be calculated, along with its corresponding P-value and S-value, and plotted to form a function. (2-8)

This can be accomplished using the concurve package in R.

Install the Package from CRAN

install.packages("concurve")

Install the Developer Version

library(devtools)
install_github("zadchow/concurve")

Check out the Articles to See the Package in Action.

"Statistical software enables and promotes cargo-cult statistics. Marketing and adoption of statistical software are driven by ease of use and the range of statistical routines the software implements. Offering complex and “modern” methods provides a competitive advantage. And some disciplines have in effect standardised on particular statistical software, often proprietary software.

Statistical software does not help you know what to compute, nor how to interpret the result. It does not offer to explain the assumptions behind methods, nor does it flag delicate or dubious assumptions. It does not warn you about multiplicity or p-hacking. It does not check whether you picked the hypothesis or analysis after looking at the data, nor track the number of analyses you tried before arriving at the one you sought to publish – another form of multiplicity. The more “powerful” and “user-friendly” the software is, the more it invites cargo-cult statistics." - Stark & Saltelli, 2018

References

1. Stark PB, Saltelli A. Cargo-cult statistics and scientific crisis. Significance. 2018;15(4):40-43.
2. Poole C. Beyond the confidence interval. Am J Public Health. 1987;77(2):195-199.
3. Sullivan KM, Foster DA. Use of the confidence interval function. Epidemiology. 1990;1(1):39-42.
4. Rothman KJ, Greenland S, Lash TL. Modern epidemiology. 2012.
5. Singh K, Xie M, Strawderman WE. Confidence distribution (CD) – distribution estimator of a parameter. arXiv [mathST]. 2007.
6. Schweder T, Hjort NL. Confidence and Likelihood*. Scand J Stat. 2002;29(2):309-332.
7. Amrhein V, Trafimow D, Greenland S. Inferential Statistics as Descriptive Statistics: There is No Replication Crisis if We Don’t Expect Replication. Am Stat. 2018
8. Greenland S. Valid P-values Behave Exactly As They Should. Some misleading criticisms of P-values and their resolution with S-values. Am Stat. 2018;18(136).