aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorKen Kellner <ken@kenkellner.com>2023-12-05 15:13:52 -0500
committerKen Kellner <ken@kenkellner.com>2023-12-05 15:13:52 -0500
commit9edbc8ce1fa292aa51f2418f1c1c0d9f298b5e82 (patch)
tree6dce0ec45685381cd8952fe0b3f4bc02204bc90c
parent0ea3e15296e68a51531e9dca142dbe9a1a376b3e (diff)
Use new input checking with jagsbasic
-rw-r--r--R/jagsbasic.R24
1 files changed, 15 insertions, 9 deletions
diff --git a/R/jagsbasic.R b/R/jagsbasic.R
index 17d465a..bd7e3f0 100644
--- a/R/jagsbasic.R
+++ b/R/jagsbasic.R
@@ -2,13 +2,20 @@
jags.basic <- function(data,inits=NULL,parameters.to.save,model.file,n.chains,n.adapt=NULL,n.iter,n.burnin=0,n.thin=1,
modules=c('glm'),factories=NULL,parallel=FALSE,n.cores=NULL,DIC=TRUE,seed=NULL,save.model=FALSE,verbose=TRUE){
- #Pass input data and parameter list through error check / processing
- data.check <- process.input(data,parameters.to.save,inits,n.chains,n.iter,n.burnin,n.thin,n.cores,DIC=DIC,
- verbose=verbose,parallel=parallel,seed=seed)
- data <- data.check$data
- parameters.to.save <- data.check$params
- inits <- data.check$inits
- if(parallel){n.cores <- data.check$n.cores}
+ if(!is.null(seed)){
+ stop("The seed argument is no longer supported, use set.seed() instead", call.=FALSE)
+ }
+
+ # Check input data
+ inps_check <- process_input(data=data, params=parameters.to.save, inits=inits,
+ n_chains=n.chains, n_adapt=n.adapt, n_iter=n.iter,
+ n_burnin=n.burnin, n_thin=n.thin, n_cores=n.cores,
+ DIC=DIC, quiet=!verbose, parallel=parallel)
+ data <- inps_check$data
+ parameters.to.save <- inps_check$params
+ inits <- inps_check$inits
+ mcmc.info <- inps_check$mcmc.info
+ if(parallel) n.cores <- inps_check$mcmc.info$n.cores
#Save start time
start.time <- Sys.time()
@@ -48,8 +55,7 @@ jags.basic <- function(data,inits=NULL,parameters.to.save,model.file,n.chains,n.
}
#Get more info about MCMC run
- end.time <- Sys.time()
- time <- round(as.numeric(end.time-start.time,units="mins"),digits=3)
+ time <- round(as.numeric(Sys.time()-start.time,units="mins"),digits=3)
if(verbose){cat('MCMC took',time,'minutes.\n')}
if(save.model){