mbnma.nodesplit: Node-splitting model for testing consistency at the treatment...

Description Usage Arguments Details Value Methods (by generic) Examples

View source: R/inconsistency.functions.R

Description

Splits contributions for a given set of treatment comparisons into direct and indirect evidence. A discrepancy between the two suggests that the consistency assumption required for NMA and MBNMA may violated.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
mbnma.nodesplit(
  network,
  fun = "linear",
  user.fun = NULL,
  beta.1 = "rel",
  beta.2 = "rel",
  beta.3 = "rel",
  beta.4 = "rel",
  method = "common",
  knots = 3,
  comparisons = NULL,
  incldr = TRUE,
  ...
)

## S3 method for class 'nodesplit'
plot(x, plot.type = NULL, ...)

Arguments

network

An object of class mbnma.network.

fun

A character vector specifying a functional form to be assigned to the dose-response. Options are given in details.

user.fun

A formula specifying any relationship including dose and one/several of: beta.1, beta.2, beta.3, beta.4.

beta.1

Refers to dose-parameter(s) specified within the dose-response function(s). Can take either "rel", "common", "random", or be assigned a numeric value (see details).

beta.2

Refers to dose-parameter(s) specified within the dose-response function(s). Can take either "rel", "common", "random", or be assigned a numeric value (see details).

beta.3

Refers to dose-parameter(s) specified within the dose-response function(s). Can take either "rel", "common", "random", or be assigned a numeric value (see details).

beta.4

Refers to dose-parameter(s) specified within the dose-response function(s). Can take either "rel", "common", "random", or be assigned a numeric value (see details).

method

Can take either "common" or "random" to indicate whether relative effects should be modelled with between-study heterogeneity or not (see details).

knots

The number/location of knots if a restricted cubic spline dose-response function is fitted (fun="rcs"). If a single number is given it indicates the number of knots (they will be equally spaced across the range of doses). If a numeric vector is given it indicates the location of the knots. Minimum number of knots is 3.

comparisons

A matrix specifying the comparisons to be split (one row per comparison). The matrix must have two columns indicating each treatment for each comparison. Values can either be character (corresponding to the treatment names given in network) or numeric (corresponding to treatment codes within the network - note that these may change if drop.discon = TRUE).

incldr

A boolean object indicating whether or not to allow for indirect evidence contributions via the dose-response relationship. This can be used when node-splitting in dose-response MBNMA to allow for a greater number of potential loops in which to check for consistency.

...

Arguments to be sent to ggplot2::ggplot()

x

An object of class("nodesplit")

plot.type

A character string that can take the value of "forest" to plot only forest plots, "density" to plot only density plots, or left as NULL (the default) to plot both types of plot.

Details

The S3 method plot() on an nodesplit object generates either forest plots of posterior medians and 95\% credible intervals, or density plots of posterior densities for direct and indirect evidence.

Value

Plots the desired graph(s) and returns an object (or list of object if plot.type=NULL) of class(c("gg", "ggplot"))

Methods (by generic)

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
# Using the triptans data
network <- mbnma.network(HF2PPITT)

split <- mbnma.nodesplit(network, fun="emax", likelihood = "binomial", link="logit",
  method="common")


#### To perform nodesplit on selected comparisons ####

# Check for closed loops of treatments with independent evidence sources
# Including indirect evidence via the dose-response relationship
loops <- inconsistency.loops(network$data.ab, incldr=TRUE)

# This...
single.split <- mbnma.nodesplit(network, fun="exponential", likelihood = "binomial", link="logit",
             method="random", comparisons=rbind(c("sumatriptan_1", "almotriptan_1")))

#...is the same as...
single.split <- mbnma.nodesplit(network, fun="exponential", likelihood = "binomial", link="logit",
             method="random", comparisons=rbind(c(6, 12)))


# Plot results
plot(split, plot.type="density") # Plot density plots of posterior densities
plot(split, plot.type="forest") # Plot forest plots of direct and indirect evidence

# Print and summarise results
print(split)
summary(split) # Generate a data frame of summary results

MBNMAdose documentation built on Sept. 13, 2020, 5:08 p.m.