# This file takes too long testthat::skip_on_cran() data(tips) facethistBindwidth1 <- list(combo = wrap("facethist", binwidth = 1)) facethistBindwidth1Duo <- list( comboHorizontal = wrap("facethist", binwidth = 1), comboVertical = wrap("facethist", binwidth = 1) ) test_that("structure", { expect_null <- function(x) { expect_true(is.null(x)) } expect_obj <- function(x) { expect_s3_class(x$data, "data.frame") expect_type(x$plots, "list") expect_equal(length(x$plots), ncol(tips)^2) expect_null(x$title) expect_null(x$xlab) expect_null(x$ylab) expect_type(x$xAxisLabels, "character") expect_type(x$yAxisLabels, "character") expect_type(x$showXAxisPlotLabels, "logical") expect_type(x$showYAxisPlotLabels, "logical") expect_null(x$legend) expect_type(x$byrow, "logical") expect_null(x$gg) expect_true("gg" %in% names(x)) } expect_obj(ggduo(tips)) expect_obj(ggpairs(tips)) }) test_that("columns", { expect_obj <- function(pm, columnsX, columnsY) { expect_equal(length(pm$plots), length(columnsX) * length(columnsY)) expect_equal(pm$xAxisLabels, columnsX) expect_equal(pm$yAxisLabels, columnsY) expect_equal(pm$ncol, length(columnsX)) expect_equal(pm$nrow, length(columnsY)) } columnsUsed <- c("total_bill", "tip", "sex") pm <- ggpairs(tips, columns = columnsUsed) expect_obj(pm, columnsUsed, columnsUsed) columnsX <- c("total_bill", "tip", "sex") columnsY <- c("smoker", "day", "time", "size") pm <- ggduo(tips, columnsX, columnsY) expect_obj(pm, columnsX, columnsY) }) test_that("column labels", { expect_obj <- function(pm, columnLabelsX, columnLabelsY) { expect_equal(pm$xAxisLabels, columnLabelsX) expect_equal(pm$yAxisLabels, columnLabelsY) } columnTitles <- c("A", "B", "C") pm <- ggpairs(tips, 1:3, columnLabels = columnTitles) expect_obj(pm, columnTitles, columnTitles) columnTitles <- c("Total Bill %", "Tip 123456", "Sex ( /a asdf)") pm <- ggpairs(tips, 1:3, columnLabels = columnTitles) expect_obj(pm, columnTitles, columnTitles) columnLabelsX <- c("Total Bill %", "Tip 123456", "Sex ( /a asdf)") columnLabelsY <- c("Smoker !#@", "Day 678", "1", "NULL") pm <- ggduo(tips, 1:3, 4:7, columnLabelsX = columnLabelsX, columnLabelsY = columnLabelsY) expect_obj(pm, columnLabelsX, columnLabelsY) }) test_that("character", { expect_obj <- function(pm) { expect_true(is.factor(pm$data$sex)) expect_true(is.factor(pm$data$smoker)) } tips2 <- tips tips2$sex <- as.character(tips2$sex) tips2$smoker <- as.character(tips2$smoker) expect_obj(ggpairs(tips2)) expect_obj(ggduo(tips2)) }) test_that("upper/lower/diag = blank", { columnsUsed <- 1:3 au <- ggpairs(tips, columnsUsed, upper = "blank") ad <- ggpairs(tips, columnsUsed, diag = "blank") al <- ggpairs(tips, columnsUsed, lower = "blank") for (i in 1:3) { for (j in 1:3) { if (i < j) { expect_true(is_blank_plot(au[i, j])) expect_false(is_blank_plot(ad[i, j])) expect_false(is_blank_plot(al[i, j])) } if (i > j) { expect_false(is_blank_plot(au[i, j])) expect_false(is_blank_plot(ad[i, j])) expect_true(is_blank_plot(al[i, j])) } if (i == j) { expect_false(is_blank_plot(au[i, j])) expect_true(is_blank_plot(ad[i, j])) expect_false(is_blank_plot(al[i, j])) } } } a <- ggpairs(tips, columnsUsed) a[1, 1] <- ggplot(tips, aes(total_bill)) + geom_histogram() expect_false(is_blank_plot(a[1, 1])) }) test_that("stops", { expect_warning( { pm <- ggpairs(tips, axisLabels = "not_a_chosen", lower = facethistBindwidth1) }, "'axisLabels' not in " ) # nolint expect_warning( { pm <- ggduo(tips, axisLabels = "not_a_chosen", types = facethistBindwidth1Duo) }, "'axisLabels' not in " ) # nolint expect_warning( { pm <- ggpairs(tips, color = "sex") }, "Extra arguments: " ) # nolint expect_warning( { pm <- ggduo(tips, 2:3, 2:3, types = list(combo = "facetdensity")) }, "Setting:\n\ttypes" ) # nolint expect_error( { ggpairs(tips, columns = c("tip", "day", "not in tips")) }, "Columns in 'columns' not found in data" ) # nolint expect_error( { ggduo(tips, columnsX = c("tip", "day", "not in tips"), columnsY = "smoker") }, "Columns in 'columnsX' not found in data" ) # nolint expect_error( { ggduo(tips, columnsX = c("tip", "day", "smoker"), columnsY = "not in tips") }, "Columns in 'columnsY' not found in data" ) # nolint expect_warning( { pm <- ggpairs(tips, legends = TRUE) }, "'legends' will be deprecated" ) # nolint expect_error( { ggpairs(tips, params = c(size = 2)) }, "'params' is a deprecated" ) # nolint expect_error( { ggpairs(tips, columns = 1:10) }, "Make sure your numeric 'columns' values are less than or equal to" ) # nolint expect_error( { ggduo(tips, columnsX = 1:10) }, "Make sure your numeric 'columnsX' values are less than or equal to" ) # nolint expect_error( { ggduo(tips, columnsY = 1:10) }, "Make sure your numeric 'columnsY' values are less than or equal to" ) # nolint expect_error( { ggpairs(tips, columns = -5:5) }, "Make sure your numeric 'columns' values are positive" ) # nolint expect_error( { ggduo(tips, columnsX = -5:5) }, "Make sure your numeric 'columnsX' values are positive" ) # nolint expect_error( { ggduo(tips, columnsY = -5:5) }, "Make sure your numeric 'columnsY' values are positive" ) # nolint expect_error( { ggpairs(tips, columns = (2:10) / 2) }, "Make sure your numeric 'columns' values are integers" ) # nolint expect_error( { ggduo(tips, columnsX = (2:10) / 2) }, "Make sure your numeric 'columnsX' values are integers" ) # nolint expect_error( { ggduo(tips, columnsY = (2:10) / 2) }, "Make sure your numeric 'columnsY' values are integers" ) # nolint expect_error( { ggpairs(tips, columns = 1:3, columnLabels = c("A", "B", "C", "Extra")) }, "The length of the 'columnLabels' does not match the length of the 'columns'" ) # nolint expect_error( { ggduo(tips, columnsX = 1:3, columnLabelsX = c("A", "B", "C", "Extra")) }, "The length of the 'columnLabelsX' does not match the length of the 'columnsX'" ) # nolint expect_error( { ggduo(tips, columnsY = 1:3, columnLabelsY = c("A", "B", "C", "Extra")) }, "The length of the 'columnLabelsY' does not match the length of the 'columnsY'" ) # nolint expect_error( { ggpairs(tips, upper = c("not_a_list")) }, "'upper' is not a list" ) # nolint expect_error( { ggpairs(tips, diag = c("not_a_list")) }, "'diag' is not a list" ) # nolint expect_error( { ggpairs(tips, lower = c("not_a_list")) }, "'lower' is not a list" ) # nolint expect_error( { ggduo(tips, types = c("not_a_list")) }, "'types' is not a list" ) # nolint # # couldn't get correct error message # # variables: 'colour' have non standard format: 'total_bill + tip'. # expect_error({ # ggpairs(tips, mapping = ggplot2::aes(color = total_bill + tip)) # }, "variables\\: \"colour\" have non standard format") # nolint # expect_error({ # ggduo(tips, mapping = ggplot2::aes(color = total_bill + tip)) # }, "variables\\: \"colour\" have non standard format") # nolint errorString <- "'aes_string' is a deprecated element" expect_error( { ggpairs(tips, upper = list(aes_string = ggplot2::aes(color = day))) }, errorString ) # nolint expect_error( { ggpairs(tips, lower = list(aes_string = ggplot2::aes(color = day))) }, errorString ) # nolint expect_error( { ggpairs(tips, diag = list(aes_string = ggplot2::aes(color = day))) }, errorString ) # nolint expect_error( { ggduo(tips, types = list(aes_string = ggplot2::aes(color = day))) }, errorString ) # nolint expect_diag_warn <- function(key, value) { warnString <- str_c("Changing diag\\$", key, " from '", value, "' to '", value, "Diag'") diagObj <- list() diagObj[[key]] <- value expect_warning( { pm <- ggpairs(tips, diag = diagObj) }, warnString ) } # diag # continuous # densityDiag # barDiag # blankDiag # discrete # barDiag # blankDiag expect_diag_warn("continuous", "density") expect_diag_warn("continuous", "bar") expect_diag_warn("continuous", "blank") expect_diag_warn("discrete", "bar") expect_diag_warn("discrete", "blank") }) test_that("cardinality", { expect_silent(stop_if_high_cardinality(tips, 1:ncol(tips), NULL)) expect_silent(stop_if_high_cardinality(tips, 1:ncol(tips), FALSE)) expect_error( stop_if_high_cardinality(tips, 1:ncol(tips), "not numeric"), "'cardinality_threshold' should" ) expect_error( stop_if_high_cardinality(tips, 1:ncol(tips), 2), "Column 'day' has more levels" ) }) test_that("blank types", { columnsUsed <- 1:3 pmUpper <- ggpairs(tips, columnsUsed, upper = "blank", lower = facethistBindwidth1) pmDiag <- ggpairs(tips, columnsUsed, diag = "blank", lower = facethistBindwidth1) pmLower <- ggpairs(tips, columnsUsed, lower = "blank") for (i in columnsUsed) { for (j in columnsUsed) { if (i < j) { # upper expect_true(is_blank_plot(pmUpper[i, j])) expect_false(is_blank_plot(pmDiag[i, j])) expect_false(is_blank_plot(pmLower[i, j])) } else if (i > j) { # lower expect_false(is_blank_plot(pmUpper[i, j])) expect_false(is_blank_plot(pmDiag[i, j])) expect_true(is_blank_plot(pmLower[i, j])) } else { # diag expect_false(is_blank_plot(pmUpper[i, j])) expect_true(is_blank_plot(pmDiag[i, j])) expect_false(is_blank_plot(pmLower[i, j])) } } } columnsUsedX <- 1:3 columnsUsedY <- 4:5 pmDuo <- ggduo(tips, columnsUsedX, columnsUsedY, types = "blank") for (i in seq_along(columnsUsedX)) { for (j in seq_along(columnsUsedY)) { expect_true(is_blank_plot(pmDuo[j, i])) } } }) test_that("axisLabels", { expect_axis_labels <- function(pm, prefix, axisLabel) { expect_true(is.null(pm$showStrips)) if (axisLabel == "show") { expect_true(pm$showXAxisPlotLabels) expect_true(pm$showYAxisPlotLabels) expect_false(is.null(pm$xAxisLabels)) expect_false(is.null(pm$yAxisLabels)) } else if (axisLabel == "internal") { for (i in 1:(pm$ncol)) { p <- pm[i, i] expect_true(inherits(p$layers[[1]]$geom, "GeomText")) expect_true(inherits(p$layers[[2]]$geom, "GeomText")) expect_equal(length(p$layers), 2) } expect_false(pm$showXAxisPlotLabels) expect_false(pm$showYAxisPlotLabels) expect_true(is.null(pm$xAxisLabels)) expect_true(is.null(pm$yAxisLabels)) } else if (axisLabel == "none") { expect_false(pm$showXAxisPlotLabels) expect_false(pm$showYAxisPlotLabels) expect_false(is.null(pm$xAxisLabels)) expect_false(is.null(pm$yAxisLabels)) } vdiffr::expect_doppelganger(paste0("axisLabels-", prefix, "-", axisLabel), pm) } fn <- function(axisLabels) { pm <- ggpairs( iris, c(3, 4, 5, 1), upper = "blank", lower = facethistBindwidth1, axisLabels = axisLabels, title = str_c("axisLabels = ", axisLabels), progress = FALSE ) pm } for (axisLabels in c("show", "internal", "none")) { expect_axis_labels(fn(axisLabels), "ggpairs", axisLabels) } plots <- ggpairs(iris, 1:3)$plots for (val in c(TRUE, FALSE)) { pm <- ggmatrix( plots, 3, 3, showAxisPlotLabels = val ) expect_equal(pm$showXAxisPlotLabels, val) expect_equal(pm$showYAxisPlotLabels, val) } fn <- function(axisLabels) { a <- ggduo( iris, c(4, 5), c(5, 1), types = facethistBindwidth1Duo, axisLabels = axisLabels, title = str_c("axisLabels = ", axisLabels) ) a } for (axisLabels in c("show", "none")) { expect_axis_labels(fn(axisLabels), "ggduo", axisLabels) } }) test_that("strips and axis", { # axis should line up with left side strips pm <- ggpairs( tips, c(3, 1, 4), showStrips = TRUE, title = "Axis should line up even if strips are present", lower = list(combo = wrap("facethist", binwidth = 1)) ) vdiffr::expect_doppelganger("show-strips", pm) # default behavior. tested in other places # expect_silent({ # pm <- ggpairs(tips, c(3, 1, 4), showStrips = FALSE) # print(pm) # }) }) test_that("dates", { startDt <- as.POSIXct("2000-01-01", tz = "UTC") endDt <- as.POSIXct("2000-04-01", tz = "UTC") dts <- seq(startDt, endDt, 86400) # 86400 = as.numeric(ddays(1)) x <- data.frame( date = dts, x1 = rnorm(length(dts)), x2 = rnorm(length(dts)), cat = sample(c("a", "b", "c"), length(dts), replace = TRUE) ) class(x) <- c("NOT_data.frame", "data.frame") a <- ggpairs( x, c(2, 1, 4, 3), mapping = ggplot2::aes(color = cat), lower = "blank", diag = list(continuous = "densityDiag"), upper = list(continuous = "cor") ) p <- a[1, 2] expect_true(inherits(p$layers[[1]]$geom, "GeomText")) expect_true(inherits(p$layers[[2]]$geom, "GeomText")) expect_equal(length(p$layers), 2) a <- ggpairs( x, c(2, 1, 4, 3), mapping = ggplot2::aes(color = cat), lower = "blank", diag = list(continuous = "barDiag"), upper = list(continuous = "cor") ) p <- a[1, 1] expect_true(inherits(p$layers[[1]]$geom, "GeomBar")) expect_equal(length(p$layers), 1) }) test_that("mapping", { pm <- ggpairs(tips, mapping = 1:3) expect_equal(pm$xAxisLabels, names(tips)[1:3]) pm <- ggpairs(tips, columns = 1:3) expect_equal(pm$xAxisLabels, names(tips)[1:3]) expect_error( { ggpairs(tips, columns = 1:3, mapping = 1:3) }, "'mapping' should not be numeric" ) # nolint }) test_that("user functions", { p0 <- ggally_points(tips, ggplot2::aes(x = total_bill, y = tip)) pm1 <- ggpairs(tips, 1:2, lower = list(continuous = "points")) p1 <- pm1[2, 1] pm2 <- ggpairs(tips, 1:2, lower = list(continuous = ggally_points)) p2 <- pm2[2, 1] expect_equal_plots <- function(x, y) { expect_equal(length(x$layers), 1) expect_equal(length(y$layers), 1) expect_true( "GeomPoint" %in% class(x$layers[[1]]$geom) ) expect_true( "GeomPoint" %in% class(y$layers[[1]]$geom) ) expect_equal(x$labels, list(x = "total_bill", y = "tip")) expect_equal(x$labels, y$labels) } expect_equal_plots(p0, p1) expect_equal_plots(p0, p2) }) test_that("NA data", { expect_is_na_plot <- function(p) { expect_true(identical(as.character(p$data$label), "NA")) expect_true(inherits(p$layers[[1]]$geom, "GeomText")) expect_equal(length(p$layers), 1) } expect_not_na_plot <- function(p) { expect_false(identical(as.character(p$data$label), "NA")) } expect_is_blank <- function(p) { expect_true(is_blank_plot(p)) } dd <- data.frame(x = c(1:5, rep(NA, 5)), y = c(rep(NA, 5), 2:6), z = 1:10, w = NA) pm <- ggpairs(dd) test_pm <- function(pm, na_mat) { for (i in 1:4) { for (j in 1:4) { if (na_mat[i, j]) { expect_is_na_plot(pm[i, j]) } else { if (j == 3 && i < 3) { expect_warning( { p <- pm[i, j] }, "Removed 5 rows" ) } else { p <- pm[i, j] } expect_not_na_plot(p) } } } } na_mat <- matrix(FALSE, ncol = 4, nrow = 4) na_mat[1, 2] <- TRUE na_mat[2, 1] <- TRUE na_mat[1:4, 4] <- TRUE na_mat[4, 1:4] <- TRUE test_pm(pm, na_mat) }) test_that("strip-top and strip-right", { data(tips) double_strips <- function(data, mapping, ...) { dt <- plyr::count(data, c(mapping_string(mapping$x), mapping_string(mapping$y))) ggplot(dt, aes(xmin = 0.25, xmax = 0.75, ymin = 1, ymax = freq)) + geom_rect() + ggplot2::facet_grid(paste0(mapping_string(mapping$y), " ~ ", mapping_string(mapping$x))) + ggplot2::scale_x_continuous(breaks = 0.5, labels = NULL) } pm <- ggpairs( tips, 3:6, lower = "blank", diag = "blank", upper = list(discrete = double_strips), progress = FALSE ) vdiffr::expect_doppelganger("nested-strips-default", pm) pm <- ggpairs( tips, 3:6, lower = "blank", diag = "blank", upper = list(discrete = double_strips), showStrips = TRUE, progress = FALSE ) vdiffr::expect_doppelganger("nested-strips-true", pm) }) test_that("subtypes", { testthat::skip_on_cran() testthat::skip_if_not_installed("Hmisc") # list of the different plot types to check # continuous # points # smooth # smooth_loess # density # cor # blank # combo # box # dot plot # facethist # facetdensity # denstrip # blank # discrete # ratio # facetbar # blank gn <- function(x) { fnName <- attr(x, "name") ifnull(fnName, x) } ggpairs_fn1 <- function(title, types, diag, ...) { ggpairs( tips, 1:4, axisLabels = "show", title = paste( "upper = c(cont = ", gn(types$continuous), ", combo = ", gn(types$combo), ", discrete = ", gn(types$discrete), "); diag = c(cont = ", gn(diag$continuous), ", discrete = ", gn(diag$discrete), ")", sep = "" ), upper = types, lower = types, diag = diag, progress = FALSE, ... ) + ggplot2::theme(plot.title = ggplot2::element_text(size = 9)) } ggpairs_fn2 <- function(...) { ggpairs_fn1(..., mapping = ggplot2::aes(color = day), legend = c(1, 3)) } ggduo_fn1 <- function(title, types, diag, ...) { types$comboHorizontal <- types$combo types$comboVertical <- types$combo types$combo <- NULL ggduo( tips, 1:3, 1:4, axisLabels = "show", title = paste( "types = c(cont = ", gn(types$continuous), ", combo = ", gn(types$comboHorizontal), ", discrete = ", gn(types$discrete), ")", sep = "" ), types = types, progress = FALSE, ... ) + ggplot2::theme(plot.title = ggplot2::element_text(size = 9)) } ggduo_fn2 <- function(...) { ggduo_fn1(..., mapping = ggplot2::aes(color = day), legend = 3) + theme(legend.position = "bottom") } # re ordered the subs so that density can have no binwidth param conSubs <- list( "autopoint", "density", "points", "smooth", "smooth_lm", "smooth_loess", "cor", "blank" ) comSubs <- list( "autopoint", "box", "dot", "box_no_facet", "dot_no_facet", wrap("facethist", binwidth = 1), "facetdensity", "facetdensitystrip", # "summarise_by", # Issues with grid printing wrap("denstrip", binwidth = 1), "blank" ) disSubs <- list( "autopoint", "colbar", "count", "cross", "crosstable", "facetbar", "ratio", "rowbar", "table", # "trends", # Issues with grid printing "blank" ) conDiagSubs <- c("autopointDiag", "densityDiag", wrap("barDiag", binwidth = 1), "blankDiag") disDiagSubs <- c("autopointDiag", "barDiag", "countDiag", "tableDiag", "blankDiag") # for (fn in list(ggpairs_fn1, ggpairs_fn2, ggduo_fn1, ggduo_fn2)) { for (fn_info in list( list(fn = ggpairs_fn1, title = "ggpairs"), list(fn = ggpairs_fn2, title = "ggpairs_color"), list(fn = ggduo_fn1, title = "ggduo"), list(fn = ggduo_fn2, title = "ggduo_color") )) { fn <- fn_info$fn fn_name <- fn_info$title for (i in 1:max(c( length(conSubs), length(comSubs), length(disSubs), length(conDiagSubs), length(disDiagSubs) ))) { conSub <- if (i <= length(conSubs)) conSubs[[i]] else "blank" comSub <- if (i <= length(comSubs)) comSubs[[i]] else "blank" disSub <- if (i <= length(disSubs)) disSubs[[i]] else "blank" diagConSub <- if (i <= length(conDiagSubs)) conDiagSubs[[i]] else "blankDiag" diagDisSub <- if (i <= length(disDiagSubs)) disDiagSubs[[i]] else "blankDiag" # print(list( # fn_num = fn_num, # types = list( # continuous = conSub, # combo = comSub, # discrete = disSub # ), # diag = list( # continuous = diagConSub, # discrete = diagDisSub # ) # )) # expect_silent({ pm <- fn( types = list( continuous = conSub, combo = comSub, discrete = disSub ), diag = list( continuous = diagConSub, discrete = diagDisSub ) ) }) type_name <- function(x) { if (is.function(x)) { sub("ggally_", "", attr(x, "name")) } else { x } } type_names <- vapply(c(conSub, comSub, disSub, diagConSub, diagDisSub), type_name, character(1)) if (all(grepl("blank", type_names))) { # vdiffr can't handle blank plots next } pm_name <- paste0(type_names, collapse = "-") pm_name <- paste0(fn_name, "-", pm_name) tryCatch( { set.seed(123456) # keep jitter consistent suppressWarnings({ built_pm <- ggmatrix_gtable(pm) }) vdiffr::expect_doppelganger(pm_name, built_pm) }, error = function(e) { message("failed to create doppelganger: ", pm_name) print(e) barret <<- pm expect_silent(print(c("failed to create doppelganger", pm_name))) } ) } } expect_error({ ggpairs(tips, 1:2, lower = "blank", diag = "blank", upper = list(continuous = "BAD_TYPE")) }) }) # pm <- ggpairs(tips, upper = "blank") # # pm # # Custom Example # pm <- ggpairs( # tips[, c(1, 3, 4, 2)], # upper = list(continuous = "density", combo = "box"), # lower = list(continuous = "points", combo = "dot") # ) # # pm # # Use sample of the diamonds data # data(diamonds, package = "ggplot2") # diamonds.samp <- diamonds[sample(1:dim(diamonds)[1], 200), ] # # Custom Example # pm <- ggpairs( # diamonds.samp[, 1:5], # upper = list(continuous = "density", combo = "box"), # lower = list(continuous = "points", combo = "dot"), # color = "cut", # alpha = 0.4, # title = "Diamonds" # ) # # pm # # Will plot four "Incorrect Plots" # bad_plots <- ggpairs( # tips[, 1:3], # upper = list(continuous = "wrongType1", combo = "wrongType2"), # lower = list(continuous = "IDK1", combo = "IDK2", discrete = "mosaic"), # ) # # bad_plots # # Only Variable Labels on the diagonal (no axis labels) # pm <- ggpairs(tips[, 1:3], axisLabels = "internal") # # pm # # Only Variable Labels on the outside (no axis labels) # pm <- ggpairs(tips[, 1:3], axisLabels = "none") # # pm # # Custom Examples # custom_car <- ggpairs(mtcars[, c("mpg", "wt", "cyl")], upper = "blank", title = "Custom Example") # #' # ggplot example taken from example(geom_text) # #' plot <- ggplot2::ggplot(mtcars, ggplot2::aes(x = wt, y = mpg, label = rownames(mtcars))) # #' plot <- plot + # #' ggplot2::geom_text(ggplot2::aes(colour = factor(cyl)), size = 3) + # #' ggplot2::scale_colour_discrete(l = 40) # #' custom_car <- putPlot(custom_car, plot, 1, 2) # #' personal_plot <- ggally_text( # #' "ggpairs allows you\nto put in your\nown plot.\nLike that one.\n <---" # #' ) # #' custom_car <- putPlot(custom_car, personal_plot, 1, 3) # #' # custom_car