Function for plotting the cross-validated model selection/tuning profiles of a sbh object.
It uses the user's choice of cross-validation criterion statistics among the Log Hazard Ratio (LHR),
Log-Rank Test (LRT) or Concordance Error Rate (CER). The function plots (as it applies) both profiles
of cross-validation criterion as a function of variables screening size (cardinal subset of top-screened
variables in the PRSP variable screening procedure), and peeling length (number of peeling steps
of the peeling sequence in the inner loop of the PRSP or PRGSP algorithm).
plot_profile(object,
main = NULL,
xlim = NULL,
ylim = NULL,
add.sd = TRUE,
add.profiles = TRUE,
add.caption = TRUE,
text.caption = c("Mean","Std. Error"),
pch = 20,
col = 1,
lty = 1,
lwd = 0.5,
cex = 0.5,
device = NULL,
file = "Profile Plots",
path = getwd(),
horizontal = FALSE,
width = 8.5,
height = 5.0, ...)Object of class sbh as generated by the main function sbh.
Character vector. Main Title.
Defaults to NULL.
Numeric vector of length 2. The x limits [x1, x2] of the plot.
Defaults to NULL.
Numeric vector of length 2. The y limits [y1, y2] of the plot.
Defaults to NULL.
Logical scalar. Shall the standard error bars be plotted?
Defaults to TRUE.
Logical scalar. Shall the individual profiles (for all replicates) be plotted?
Defaults to TRUE.
Logical scalar. Should the caption be plotted?
Defaults to TRUE.
Character vector of caption content.
Defaults to {"Mean","Std. Error"}.
Integer scalar of symbol number for all the profiles.
Defaults to 20.
Integer scalar of line color of the mean profile.
Defaults to 1.
Integer scalar of line type of the mean profile.
Defaults to 1.
Numeric scalar of line width of the mean profile.
Defaults to 0.5.
Numeric scalar of symbol expansion for all the profiles.
Defaults to 0.5.
Graphic display device in {NULL, "PS", "PDF"}.
Defaults to NULL (standard output screen).
Currently implemented graphic display devices are "PS" (Postscript) or "PDF" (Portable Document Format).
File name for output graphic. Defaults to "Profile Plot".
Absolute path (without final (back)slash separator). Defaults to working directory path.
Logical scalar. Orientation of the printed image.
Defaults to FALSE, that is potrait orientation.
Numeric scalar. Width of the graphics region in inches.
Defaults to 8.5.
Numeric scalar. Height of the graphics region in inches.
Defaults to 5.0.
Generic arguments passed to other plotting functions.
Invisible. None. Displays the plot(s) on the specified device.
This work made use of the High Performance Computing Resource in the Core Facility for Advanced Research Computing at Case Western Reserve University. This project was partially funded by the National Institutes of Health NIH - National Cancer Institute (R01-CA160593) to J-E. Dazard and J.S. Rao.
Model tuning is done by applying the cross-validation criterion defined by the user's choice of specific statistic. The goal is to find the optimal value of model parameters by maximization of LHR or LRT, or minimization of CER. The parameters to optimize are (i) the cardinal of top-ranked variables subsets (if the "prsp" variable screening is chosen), and (ii) the number of peeling steps of the peeling sequence (inner loop of our PRSP algorithm) in any case of variable screening method.
Currently, this is done internally for visualization purposes, but it will ultimately offer the option to be done interactively with the end-user as well for parameter choosing/model selection.
Dazard J-E. and Rao J.S. (2018). "Variable Selection Strategies for High-Dimensional Survival Bump Hunting using Recursive Peeling Methods." (in prep).
Rao J.S., Huilin Y. and Dazard J-E. (2018). "Disparity Subtyping: Bringing Precision Medicine Closer to Disparity Science." (in prep).
Diaz-Pachon D.A., Saenz J.P., Dazard J-E. and Rao J.S. (2018). "Mode Hunting through Active Information." (in press).
Diaz-Pachon D.A., Dazard J-E. and Rao J.S. (2017). "Unsupervised Bump Hunting Using Principal Components." In: Ahmed SE, editor. Big and Complex Data Analysis: Methodologies and Applications. Contributions to Statistics, vol. Edited Refereed Volume. Springer International Publishing, Cham Switzerland, p. 325-345.
Yi C. and Huang J. (2017). "Semismooth Newton Coordinate Descent Algorithm for Elastic-Net Penalized Huber Loss Regression and Quantile Regression." J. Comp Graph. Statistics, 26(3):547-557.
Dazard J-E., Choe M., LeBlanc M. and Rao J.S. (2016). "Cross-validation and Peeling Strategies for Survival Bump Hunting using Recursive Peeling Methods." Statistical Analysis and Data Mining, 9(1):12-42.
Dazard J-E., Choe M., LeBlanc M. and Rao J.S. (2015). "R package PRIMsrc: Bump Hunting by Patient Rule Induction Method for Survival, Regression and Classification." In JSM Proceedings, Statistical Programmers and Analysts Section. Seattle, WA, USA. American Statistical Association IMS - JSM, p. 650-664.
Dazard J-E., Choe M., LeBlanc M. and Rao J.S. (2014). "Cross-Validation of Survival Bump Hunting by Recursive Peeling Methods." In JSM Proceedings, Survival Methods for Risk Estimation/Prediction Section. Boston, MA, USA. American Statistical Association IMS - JSM, p. 3366-3380.
Dazard J-E. and J.S. Rao (2010). "Local Sparse Bump Hunting." J. Comp Graph. Statistics, 19(4):900-92.