plotAreaInROPE {BEST}R Documentation

Area of the posterior density in the ROPE as a function of its width.


Calculates and (optionally) plots the posterior probability mass included in the Region of Practical Equivalence (ROPE: see plot.BEST) as a function of the width of the ROPE.


plotAreaInROPE(paramSampleVec, credMass = 0.95, compVal = 0, maxROPEradius,
  n = 201, plot = TRUE, ROPEColor = "darkred", ...)



A vector of samples drawn from the target distribution; see Examples.


The probability mass to include in credible intervals.


a value for comparison with those plotted.


The maximum value of the ROPE radius (ie. half-width) to include in the plot.


The number of equally spaced points at which the area in the ROPE is to be estimated.


If FALSE, the plot will be suppressed but the values will be returned.


an optional color name such as "darkred" or a RGB specification such as "#8B0000" that controls the color used to plot the ROPE.


Other graphical parameters.


Defining a Region of Practical Equivalence (ROPE) allows decisions on whether a parameter is, for practical purposes, equivalent to a hypothetical null value, given a posterior probability density for the parameter. The null value may be considered credible if (A) 95% (say) of the probability mass lies within the ROPE, or (B) the 95% highest density interval (95% HDI) lies entirely within the ROPE.

How wide should the ROPE be? Different people at different times will have different ideas on the range of values equivalent to the null. The function plotAreaInROPE plots the probability mass lying within the ROPE for a range of widths (or rather radii or half-widths). It also shows the radius at which the HDI falls entirely within the ROPE.


Returns invisibly a list with elements:


A vector of ROPE radii from 0 to maxROPEradius.


The corresponding proportion of the posterior density included in the ROPE.


John K. Kruschke, with minor modifications by Mike Meredith.



# Generate a fake MCMC posterior for effect size and plot it:
mcmcChain <- rnorm(50000,0.03,0.025)
plotPost(mcmcChain, compVal=0, ROPE=c(-0.1, 0.1))

# How does the mass within the ROPE vary with ROPE radius?
plotAreaInROPE(mcmcChain, credMass = 0.95, compVal = 0,
  maxROPEradius = 0.15)

# Generate real MCMC chains, takes up to 1 min:
y1 <- c(4.77, 4.33, 3.59, 3.33, 2.66, 3.48)
y2 <- c(3.88, 3.55, 3.29, 2.59, 2.33, 3.59)
BESTout <- BESTmcmc(y1, y2, parallel=FALSE)

meanDiff <- BESTout$mu1 - BESTout$mu2
plotAreaInROPE(meanDiff, credMass = 0.95, compVal = 0,
  maxROPEradius = 3)

[Package BEST version 0.5.4 Index]