Multiplicity graphs

Introduction

With gsDesign version 3.1, we have added a function to support graphical multiplicity methods using ggplot2, this function is migrated to the gMCPLite package now.

The graphical method is introduced nicely in Bretz et al. (2011) and originally supported by the gMCP packages (Rohmeyer and Klinglmueller (2020)), and now gMCPLite package. It was extended to group sequential design by Maurer and Bretz (2013). While the gMCPLite package supports graphics, here we add the hGraph() function to create multiplicity graphs using the ggplot2 package as a convenience for those desiring this format. We demonstrate basic formatting in this article and demonstrate use with gMCPLite in the other article.

Use of hGraph() is organized as follows:

  • The basic graph layout
  • Specifying hypothesis names and α-allocation
  • Text formatting and location, ellipse size
  • Specifying the transition matrix between hypotheses
  • Using colors and legends

The basic graph layout

We begin with the default plot to demonstrate the basic formatting below. Hypotheses and initial α-allocation or weighting are presented in shaded ellipses with the first hypothesis specified in the upper-left part of the plot. While the user has full control of ellipse placement and color, the default is to present clockwise on a larger ellipse with gray shading. One advantage of default placement is that generally transition lines between hypotheses will not cross hypothesis ellipses, including when the graph is updated when some hypotheses are rejected. Transition weights between hypotheses are specified on directional lines between the hypothesis ellipses. Try this in a plot window. Note the effect if you change the size of the plot window. In a similar fashion, if you are using R Markdown, parameters like fig.width, fig.height, and fig.asp will affect formatting.

hGraph()

Specifying hypothesis names and α-allocation

Next we specify the number of hypotheses, hypothesis names and α-allocation or weighting. Whereas the default plot above used α-allocation adding to 0.025, here we use weights adding to 1. Note the \n character used to insert a carriage return in hypothesis name text. The character used for weights (specified in wchar) is by default α under Windows and w otherwise. Note the clockwise placement of hypotheses.

hGraph(
  nHypotheses = 3,
  nameHypotheses = c("HA\n First", "HC\n Second", "HB\n Third"),
  alphaHypotheses = c(.2, .3, .5),
  wchar = "w" # Character before weights
)

Text formatting and location, ellipse size

You can specify location of hypothesis ellipses in two ways.

  • For the first, specify in radians where on the large ellipse where the first hypothesis is placed. Setting π/2 for radianStart places the first hypothesis at the center and top of the plot (left plot).
  • Specifying x and y coordinates, right graph, allows custom placement; note that the differences in minimum and maximum x- and y-values you use as well as window size will impact formatting.

Ellipse and the text size as well as maximum significant digits for hypotheses are controlled as follows:

  • Size of text in hypothesis ellipses is controlled by the size parameter (left graph). This text is always placed at the center of the ellipse with no user control.
  • Size of ellipses is controlled by halfWid (left graph) and halfHgt (not shown; default 0.5) parameters.
  • digits specifies the maximum digits for the α-allocation (not shown here; default 5)

Box and text size as well as maximum significant digit display for transition weights are controlled as follows:

  • boxtextsize (right graph) controls text size for transition weights
  • trhw (left graph) controls half-width for transition box size
  • trhh (not shown; default 0.075) controls half-height for transition box size
  • trprop (right graph) specifies the proportional placement of transition weight boxes along the transition lines At present, this is the same for all transitions
  • arrowsize (left plot) controls the size of the arrows between ellipses

The other parameter shown here is offset (left graph). This is in radians; it increases and decreases offset of transition lines between hypothesis ellipses. If there are not any transition arrows in both directions between any pair of hypotheses, offset = 0 is a reasonable option.

grid.arrange(
  # Left graph in figure
  hGraph(
    nHypotheses = 3,
    size = 5, # Decrease hypothesis text size from default 6
    halfWid = 1.25, # Increase ellipse width from default 0.5
    trhw = 0.25, # Increase transition box sizes from default 0.075
    radianStart = pi / 2, # First hypothesis top middle
    offset = pi / 20, # Decrease offset between transition lines
    arrowsize = .03 # Increase from default 0.02
  ),
  # Right graph in figure
  hGraph(
    nHypotheses = 3,
    x = c(-1, 1, -1), # Custom placement using x and y
    y = c(1, 1, -1),
    halfWid = 0.7, # Increase ellipse width from default 0.5
    boxtextsize = 3, # Decrease box text size from default 4
    trprop = .15 # Slide transition boxes closer to initiating hypothesis
  ),
  nrow = 1
)

Using colors and legends

The default color palette for ellipse shading is monochrome (left plot below). To limit the parameters used, we have not adjusted text, ellipse and box size as above. Colors are specified as follows:

  • A fill color category; here we group the first 2 hypotheses and last 2 hypotheses (both graphs)
  • The palette specifies a palette for the colors (right graph)
grid.arrange(
  # Left graph in figure
  hGraph(
    fill = c(1, 1, 2, 2),
    alphaHypotheses = c(.2, .2, .2, .4) * .025
  ),
  # Right graph in figure
  hGraph(
    fill = c(1, 1, 2, 2),
    palette = c("pink", "lightblue"),
    alphaHypotheses = c(.2, .2, .2, .4) * .025
  ),
  nrow = 1
)

Next, we add a legend.

  • The legend is added when the default of legend.positon = "none" is replaced. Here we use a custom position specifying the relative positioning on the (0,1) scale in the x- and y-range of the graph. Commonly used options are "left", "right", "top", or ’“bottom”`
  • legend.name (default is none) specifies a name for the legend
  • labels specifies labels corresponding to the fill parameter
hGraph(
  nHypotheses = 3,
  fill = c(1, 1, 2),
  palette = c("yellow", "lightblue"),
  legend.name = "Color scheme",
  labels = c("Primary", "Secondary"),
  legend.position = c(.75, .25)
)

Specifying the transition matrix between hypotheses

Transition weights are specified in matrix format in the variable m. The matrix must be square of dimension nHypotheses and have 0 on the diagonal. Rows represent the hypotheses from which transitions start; row values should sum to 1. Columns represent the hypotheses to which transition arrows go. Any transition weight of 0 has no corresponding transition line in the figure.

hGraph(
  nHypotheses = 3,
  m = matrix(c(
    0, 1, 0,
    .2, 0, .8,
    .3, .7, 0
  ),
  nrow = 3, byrow = TRUE
  ),
)

References

Bretz, Frank, Martin Posch, Ekkehard Glimm, Florian Klinglmueller, Willi Maurer, and Kornelius Rohmeyer. 2011. “Graphical Approaches for Multiple Comparison Procedures Using Weighted Bonferroni, Simes or Parametric Tests.” Biometrical Journal 53 (6): 894–913.
Maurer, Willi, and Frank Bretz. 2013. “Multiple Testing in Group Sequential Trials Using Graphical Approaches.” Statistics in Biopharmaceutical Research 5: 311–20.
Rohmeyer, Kornelius, and Florian Klinglmueller. 2020. gMCP: Graph Based Multiple Test Procedures. https://cran.r-project.org/package=gMCP.