#18 JAGS compatibility with BLAS / LAPACK 3.12.0

[Anonymous]

The current version of R 4.4.2 ships with LAPACK / BLAS v3.12.0. When testing some R package running JAGS under the, we've noticed some incompatibility with these newer versions of LAPACK / BLAS. Here's a specific example:

library(blrm)
library(testthat)

mean <- c(-3.068, 0.564)
se <- c(2.706, 0.728)
corr <- -0.917
prior <- list(mean=mean, se=se, corr=corr)

# parameters
seeds <- 1:2
nsamples <- 10000
burn_in <- 0.2
drug_name <- "DRUG-X"
dose_unit <- "mg"
prov_dose <- c(360, 480, 720, 1080, 1440)
ref_dose <- 720
category_bound <- c(0.16, 0.33)
category_name <- c("under-dosing", "targeted-toxicity", "over-dosing")
ewoc <- 0.25

# observed cohorts
dose <- 480
n_pat <- 3
dlt <- 0

# combine prior, parameters, and observed cohorts to a list
data <- list(seeds=seeds, nsamples=nsamples, burn_in=burn_in, drug_name=drug_name,
            dose_unit=dose_unit, prov_dose=prov_dose, ref_dose=ref_dose,
            dose=dose, n_pat=n_pat, dlt=dlt, category_bound=category_bound,
            category_name=category_name, ewoc=ewoc)

# run function
observed <- blrm_mono_ss(prior=prior, data=data, output_excel=FALSE, output_pdf=FALSE)

Error in `rjags::jags.model(textConnection(modelstring), data = mydata, 
    n.chains = 2, quiet = TRUE, inits = inits.list)`: BLAS/LAPACK routine 'DSYRK ' gave error code -7
Backtrace:
    ▆
1. └─blrm::blrm_mono_ss(...) at blrm.R:64:3
2.   └─blrm (local) single_cohort(prior_para)
3.     └─rjags::jag