analyse_baSAR() crashes for various settings of variable.names

[mcol]

These crash:

analyse_baSAR(CWOSL.sub,
                source_doserate = c(0.04, 0.001),
                signal.integral = c(1, 2),
                background.integral = c(8, 10),
                method_control = list(n.chains = 1,
                                      variable.names = c('central_D', 'sigma_D', 'D', 'Q', 'a', 'b')),
                n.MCMC = 10, verbose = FALSE)

# Error in plot_matrix[4, i] : subscript out of bounds
# In addition: Warning messages:
# 1: [analyse_baSAR()] Only multiple grain data provided, automatic selection skipped
# 2: [analyse_baSAR()] [calc_OSLLxTxRatio()] Number of background channels for Tx < 25, error estimation might not be reliable # This warning occurred 3 times
# 3: In rjags::jags.model(file = textConnection(baSAR_models[[distribution]]),  :
#   Adaptation incomplete

analyse_baSAR(CWOSL.sub,
                source_doserate = c(0.04, 0.001),
                signal.integral = c(1, 2),
                background.integral = c(8, 10),
                method_control = list(n.chains = 1,
                                      variable.names = c('sigma_D')),
                n.MCMC = 10, verbose = FALSE)

# Error: [analyse_baSAR()] Plots for 'central_D' and 'sigma_D' could not be produced. You are probably monitoring the wrong variables

analyse_baSAR(CWOSL.sub,
                source_doserate = c(0.04, 0.001),
                signal.integral = c(1, 2),
                background.integral = c(8, 10),
                method_control = list(n.chains = 1,
                                      variable.names = c('central_D', 'sigma_D')),
                n.MCMC = 10, verbose = FALSE)

# Error in matrixStats::colQuantiles(x = plot_matrix, probs = c(0.25, 0.75)) :
#  Argument 'x' is of type ‘list’. Only 'integer', 'numeric', and 'logical' is supported.

This is not a crash, but the error message doesn't make sense to me, as the set of names is a subset of the default names:

analyse_baSAR(CWOSL.sub,
                source_doserate = c(0.04, 0.001),
                signal.integral = c(1, 2),
                background.integral = c(8, 10),
                method_control = list(n.chains = 1,
                                      variable.names = c('central_D', 'sigma_D', 'D')),
                n.MCMC = 10, verbose = FALSE)

# [analyse_baSAR()] Error: Wrong 'variable.names' monitored, dose responses curves could not be plotted

This is also not a crash, but we could avoid the error message with better validation:

analyse_baSAR(CWOSL.sub,
                source_doserate = c(0.04, 0.001),
                signal.integral = c(1, 2),
                background.integral = c(8, 10),
                method_control = list(n.chains = 1,
                                      variable.names = ""),
                n.MCMC = 10, verbose = FALSE)

# Error in parse.varnames(variable.names[type == t]) :
#   Invalid variable subset
¬¬¬