Plot facet-equivalence results

Usage

plot_facet_equivalence(
  x,
  diagnostics = NULL,
  facet = NULL,
  type = c("forest", "rope"),
  draw = TRUE,
  ...
)

Arguments

x: Output from analyze_facet_equivalence() or an eligible MML fit_mfrm() object. Legacy equivalence bundles must be recomputed.
diagnostics: Optional matching output from diagnose_mfrm() when x is an mfrm_fit object.
facet: Facet to analyze when x is an mfrm_fit object.
type: Plot type: "forest" (default) or "rope".
draw: If TRUE (default), draw the plot. If FALSE, return the prepared plotting data.
...: Additional graphical arguments passed to base plotting functions.

Value

Invisibly returns the plotting data and inference/covariance basis. With draw = FALSE, returns the data without drawing.

Details

Fit inputs use the same eligibility checks as analyze_facet_equivalence(). Bundle inputs display the already calculated results. Both routes require the current inference and covariance basis, including when draw = FALSE.

Plot types

"forest" shows each level's deviation from the equally weighted facet mean, with covariance-aware deviation intervals and the practical region around zero. The raw marginal measure intervals remain in the data table.
"rope" shows the normal confidence-distribution mass within that region.

Interpreting output

Both plots describe grand-mean proximity. Colors in the forest plot indicate whether the deviation interval is inside, outside, or overlaps the practical region. Neither plot establishes pairwise equivalence or a Bayesian probability. Read the pairwise TOST results for pair-specific conclusions.

Typical workflow

Run analyze_facet_equivalence() with a prespecified practical bound.
Use type = "forest" to inspect deviations and their uncertainty.
Use type = "rope" for a descriptive proximity view.

Examples

# \donttest{
toy <- load_mfrmr_data("example_core")
fit <- fit_mfrm(toy, "Person", c("Rater", "Criterion"), "Score",
                method = "MML", quad_points = 31, maxit = 150)
eq <- analyze_facet_equivalence(fit, facet = "Rater")
pdat <- plot_facet_equivalence(eq, type = "forest", draw = FALSE)
c(pdat$facet, pdat$type)
#> [1] "Rater"  "forest"
# }