SubgrPlots (version 0.1.3)

plot_nightingale: Nightingale plot

Description

This function produces a nightingale plot with a binary response.

Usage

plot_nightingale(
  dat,
  covari.sel,
  trt.sel = NULL,
  resp.sel,
  outcome.type = "binary",
  seq_by = 100,
  title = NULL,
  lab.y = NULL,
  strip = "Response",
  palette_colors = c("#faa8d2", "#80b1d3")
)

Arguments

dat

a data set

covari.sel

a vector of indices of the two covariates

trt.sel

a covariate index specifying the treatment code

resp.sel

a covariate index specifying the response variable

outcome.type

a string specifying the type of the response variable, it can be "continuous", or "binary" or "survival".

seq_by

a number specifying the breaks in the y-axis

title

a string specifying the main title.

lab.y

a character value specifying the label of the y axis

strip

title for the legend

palette_colors

the colors for the graphics, a two character vector.

Examples

# NOT RUN {
library(dplyr)
dat <- prca
levels(dat$age_group) = c("Young","Middle-aged","Old")
levels(dat$weight_group)  = c("Low","Mid","High")
comb_levels = c("Young - Low", "Young - Mid", "Young - High",
                "Middle-aged - Low", "Middle-aged - Mid", "Middle-aged - High",
                "Old - Low", "Old - Mid", "Old - High")
dat %>%
  mutate(AgeWeight = factor(sprintf("%s - %s", age_group, weight_group),
                            levels = comb_levels))  %>%
  mutate(survival = factor(ifelse(survtime > 24 , "Yes", "No"),
                           levels = c("No", "Yes"))) %>%
  mutate(rx = factor(rx, labels = c("Control", "Treatment"))) -> dat

plot_nightingale(dat = dat, covari.sel = 16, resp.sel = 17,
                 strip = "2-year survival")
plot_nightingale(dat = dat, trt.sel = 3, covari.sel = 16,
                 resp.sel = 17,
                 seq_by = 50,
                 strip = "2-year survival")
# }