netmeta (version 1.2-1)

netmatrix: Create a matrix with additional information for pairwise comparisons

Description

Auxiliary function to create a matrix with additional information for pairwise comparisons

Usage

netmatrix(
  x,
  var,
  levels,
  labels = levels,
  func = "mode",
  ties.method = "random"
)

Arguments

x

A netmeta object.

var

Variable with additional information.

levels

An optional vector of the values that var might have taken (see factor).

labels

An optional vector with labels for var (see factor).

func

A character string with the function name to summarize values within pairwise comparisons; see Details.

ties.method

A character string describing how ties are handled if func = "mode"; see Details.

Value

A matrix with the same row and column names as the adjacency matrix x$A.matrix.

Details

For each pairwise comparison, unique values will be calculated for the variable var based on the argument func: "mode" (most common value), "min" (minimum value), "max", "mean", "median", and "sum". In order to determine the most common value, the argument ties.method can be used in the case of ties with "first" meaning that the first / smallest value will be selected; similar for "last" (last / largest value) and "random" (random selection).

See Also

netmeta, netgraph.netmeta

Examples

# NOT RUN {
data(smokingcessation)
# Add variable with (fictious) risk of bias values
# with 1 = "low risk" and 2 = "high risk"
#
smokingcessation$rob <- rep(1:2, 12)

p1 <- pairwise(list(treat1, treat2, treat3),
               event = list(event1, event2, event3),
               n = list(n1, n2, n3),
               data = smokingcessation,
               sm = "OR")
net1 <- netmeta(p1, comb.fixed = FALSE, ref = "A")

# Generate network graph with information on risk of bias
#
col.rob <- netmatrix(net1, rob, ties.method = "last",
                     levels = 1:2,
                     labels = c("green", "yellow"))
#
netgraph(net1, col = col.rob,
         plastic = FALSE, thickness = "number.of.studies", multi = FALSE)

n.trts <- net1$n.trts
labs <- paste(net1$trts, " (n=", n.trts, ")", sep = "")
#
netgraph(net1, col = col.rob,
         plastic = FALSE, thickness = "number.of.studies", multi = FALSE,
         points = TRUE, col.points = "blue",
         cex.points = 6 * sqrt(n.trts / max(n.trts)),
         labels = labs)

# }