nma.compare: Consistency vs Inconsistency plot

In audrey-b/BUGSnet: Bayesian network meta-analyses in compliance with best practice and reporting guidelines

Consistency vs Inconsistency plot

Description

Plots the posterior mean deviance of a consistency model vs an inconsistency model. This plot can help identify loops where inconsistency is present. Ideally, both models will contribute approximately 1 to the posterior mean deviance.

Usage

nma.compare(consistency.model.fit, inconsistency.model.fit, ...)

Arguments

consistency.model.fit	Results of nma.fit() of an consistency model.
inconsistency.model.fit	Results of nma.fit() of an inconsistency model.
...	Graphical arguments such as main=, ylab=, and xlab= may be passed in plot().

See Also

nma.run

Examples

data(diabetes.sim)
diabetes.slr <- data.prep(
  arm.data = diabetes.sim, 
  varname.t = "Treatment", 
  varname.s = "Study"
)

#Random effects, consistency model.
#Binomial family, cloglog link. This implies that the scale will be the Hazard Ratio.
diabetes.re.c <- nma.model(
  data = diabetes.slr,
  outcome = "diabetes",
  N = "n",
  reference = "Placebo",
  family = "binomial",
  link = "cloglog",
  effects = "random",
  type = "consistency",
  time = "followup"
)
 
diabetes.re.c.res <- nma.run(
  model = diabetes.re.c,
  n.adapt = 100,
  n.burnin = 0,
  n.iter = 100
)

#Random effects, inconsistency model.
#Binomial family, cloglog link. This implies that the scale will be the Hazard Ratio.
diabetes.re.i <- nma.model(
  data = diabetes.slr,
  outcome = "diabetes", 
  N = "n",
  reference = "Placebo",
  family = "binomial",
  link = "cloglog",
  effects = "random",
  type = "inconsistency",
  time = "followup"
)
 
diabetes.re.i.res <- nma.run(
  model = diabetes.re.i,
  n.adapt = 100,
  n.burnin = 0,
  n.iter = 100
)  

# Assess model fit for a both an inconsistency model and consistency model using nma.fit()
assess.consistency <- nma.fit(diabetes.re.c.res)
assess.inconsistency <- nma.fit(diabetes.re.i.res)

#Plot the results against each other to assess inconsistency
nma.compare(assess.consistency, assess.inconsistency)

audrey-b/BUGSnet documentation built on Feb. 2, 2025, 5:10 p.m.