summary.crossnma {crossnma}R Documentation

Summary function for crossnma object


This function creates posterior summary statistics for the fitted cross network meta-analysis / meta-regression model


## S3 method for class 'crossnma'
summary(object, quantiles = c(0.025, 0.5, 0.975), exp = TRUE, ...)



An object generated by the crossnma.


A numeric vector of probabilities to present posterior summaries. The default value is c(0.025, 0.5, 0.975) for the 95% credible interval and the median.


A logical value indicating whether to exponentiate the parameters of relative treatment effect and covariate effect. Default is TRUE.


Additional arguments to be passed to summary() function


crossnma.summary returns a matrix containing the following summary statistics (in columns) for each estimated parameter:

Mean the mean of the posterior distribution

SD the standard deviation of the posterior distribution

2.5% (default) the 2.5% quantile of the posterior distribution (the lower bound of the 95% credible interval)

50% (default) the median of the posterior distribution

97.5% (default) the 97.5% quantile of the posterior distribution (the upper bound of the 95% credible interval)

Rhat Gelman-Rubin statistic. The further the value of Rhat from 1, the worse the mixing of chains and so the convergence.

n.eff An estimate of the effective sample size. The smaller the value of n.eff the greater the uncertainty associated with the corresponding parameter.


Tasnim Hamza

See Also



# We conduct a network meta-analysis assuming a random-effects
# model.
# The data comes from randomized-controlled trials and
# non-randomized studies (combined naively)
head(ipddata) # participant-level data
head(stddata) # study-level data

# Create a JAGS model
mod <- crossnma.model(treat, id, relapse, n, design, = ipddata, = stddata,
  reference = "A", trt.effect = "random", method.bias = "naive")

# Fit JAGS model
# (suppress warning 'Adaptation incomplete' due to n.adapt = 20)
fit <-
  suppressWarnings(crossnma(mod, n.adapt = 20,
    n.iter = 50, thin = 1, n.chains = 3))

# Display the output

[Package crossnma version 1.0.1 Index]