diff options
author | Ken Kellner <ken@kenkellner.com> | 2023-12-05 15:13:52 -0500 |
---|---|---|
committer | Ken Kellner <ken@kenkellner.com> | 2023-12-05 15:13:52 -0500 |
commit | 9edbc8ce1fa292aa51f2418f1c1c0d9f298b5e82 (patch) | |
tree | 6dce0ec45685381cd8952fe0b3f4bc02204bc90c | |
parent | 0ea3e15296e68a51531e9dca142dbe9a1a376b3e (diff) |
Use new input checking with jagsbasic
-rw-r--r-- | R/jagsbasic.R | 24 |
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){ |