Plot robust median ATS line for censored data

Function used by other functions to plot the Akritas-Theil-Sen (ATS) line for censored data. Only one x variable allowed. Both Y and X variables may be censored.

kenplot(
  y1,
  ycen,
  x1,
  xcen,
  atsline = FALSE,
  xnam = NULL,
  ynam = NULL,
  Title = "Akritas - Theil - Sen line",
  ylim = NULL,
  xlim = NULL,
  pch = NULL,
  cex = NULL,
  xaxs = "r",
  yaxs = "r",
  ...
)

Arguments

y1

The column of y (response variable) values plus detection limits

ycen

The y-variable indicators, where 1 (or TRUE) indicates a detection limit in the y.var column, and 0 (or FALSE) indicates a detected value in y.var.

x1

The column of x (explanatory variable) values plus detection limits

xcen

The x-variable indicators, where 1 (or TRUE) indicates a detection limit in the x.var column, and 0 (or FALSE) indicates a detected value in x.var.

atsline

Indicator of whether to draw the ATS line or not. Default is FALSE.

xnam

Custom label for the x axis of plots. Default is x variable column name.

ynam

Custom label for the y axis of plots. Default is y variable column name.

Title

Custom title for plots. Default is "Akritas - Theil - Sen line".

ylim

argument consistent with base plotting to adjust y-axis limits; see par for more details

xlim

argument consistent with base plotting to adjust x-axis limits; see par for more details

pch

argument consistent with base plotting to adjust point type, see par for more details

cex

argument consistent with base plotting to adjust point size; see par for more details

xaxs

argument consistent with base plotting to adjust x-axis type; see par for more details

yaxs

argument consistent with base plotting to adjust y-axis type; see par for more details

...

argument to adjust other base plotting functions; see par for more details

Value

Scatterplot of data plus ATS line. Censored values are drawn for both X and Y variables as dashed lines up to the detection limits.

References

Helsel, D.R., 2011. Statistics for Censored Environmental Data using Minitab and R, 2nd ed. John Wiley & Sons, USA, N.J.

See also

NADA::cenken

Examples

# \donttest{
# Both y and x are censored
data(PbHeron)
with(PbHeron, kenplot(Blood, BloodCen, Kidney, KidneyCen))


# x is not censored
data(Brumbaugh)
with(Brumbaugh, kenplot(Hg, HgCen, PctWetland,rep(0, times=length(PctWetland))))

# }