barplotMD: Barplot of frequencies, densities or both of the two...

In UncertainInterval: Uncertain Interval Methods for Three-Way Cut-Point Determination in Test Results

Description

Barplot of frequencies, densities or both of the two distributions of patients with and without the targeted condition.

Usage

barplotMD(
  ref,
  test,
  name.test = "",
  fixed.range = c(NULL, NULL),
  plot = c("frequencies", "densities", "both", "none"),
  target.condition = "Target Condition",
  position.legend = "top",
  cex.legend = 1
)

Arguments

ref	Vector of patient status with two ordered values. The first indicates the patients without the targeted condition (for instance 0), the second indicates the patients with the targeted condition (for instance 1). This order is relevant.
test	Vector of ordinal test scores. The values range from min(test) to max(test) and need to be ordered. Missing values in between are shown in the plots as gaps.
name.test	Name used for the title of the x axis.
fixed.range	Default = NULL. If test has numeric values you can set the values to cover a fixed range. This may enable the comparison of different samples with truncated test values. Default: use min(test):max(test) as the range of values.
plot	Default: 'frequencies'. Which plots to create: 'frequencies', 'densities' or 'both'.
target.condition	Default: 'Target Condition'. Name of target condition.
position.legend	Default: 'top'. Position of the legend. Most used values: "topleft", "top", "topright".
cex.legend	Default: 1. Relative size of the legend.

Value

named matrix with 2 rows and max(test)-min(test)+1 columns that provide for the position on the x-axis of each of the test values. These values can be used to draw vertical lines to indicate cutoff scores (see example).

See Also

plotMD

Examples

data(tostbegg2)
barplotMD(ref=tostbegg2$d, test=tostbegg2$y, name='Metastatic Rating', cex=1)
x.axis = barplotMD(ref=tostbegg2$d, test=tostbegg2$y, plot='densities',
name='Metastatic Rating', cex=1)
# Use x.axis to plot vertical line between test score 3 and 4
segments(x0=(x.axis[2,4]+x.axis[1,3])/2, y0=0, y1=.4, col='red')
# include zero score (in this sample empty)
barplotMD(ref=tostbegg2$d, test=tostbegg2$y, fixed.range = c(0, 5),
          plot='densities',name='Metastatic Rating', cex=1)
op = par(mfrow=c(2,1))
barplotMD(ref=tostbegg2$d, test=tostbegg2$y, plot='both',
          name='Metastatic Rating', cex.legend=.6, pos='top')
par(mfrow=op)

UncertainInterval documentation built on March 3, 2021, 1:10 a.m.