R/model-code.R

In poissonconsulting/jaggernaut: Bayesian Analyses using 'JAGS'

#' @title Get model code
#'
#' @description
#' Get the model_code component of a JAGS object.  
#' 
#' @param object a JAGS object.
#' @param ... further arguments passed to or from other methods.
#' @return The model_code component of a JAGS object.
#' @seealso \code{\link{jaggernaut}}  
#' @export
model_code <- function (object, ...) {
  UseMethod("model_code", object)
}

#' @title Set model code
#'
#' @description
#' Set the model_code component of a JAGS object.  
#' 
#' @usage
#' model_code(object) <- value
#' @param object a JAGS object.
#' @param value a character element defining the model code in the JAGS dialect
#' of the BUGS language.
#' @seealso \code{\link{jaggernaut}}  
#' @export
"model_code<-" <- function (object, value) {
  UseMethod("model_code<-", object)
}

model_code.jagr_model <- function (object, ...) {
  object$model_code
}

model_code_jagr_model <- function (object, ...) {
  stopifnot(is.jagr_model(object))
  model_code(object, ...)
}

#' @method model_code jags_model
#' @export
model_code.jags_model <- function (object, ...) {
  
  if(is_one_model(object))
    return (model_code(model(object), ...))
  
  lapply(models(object), model_code_jagr_model, ...) 
}

#' @method model_code jags_analysis
#' @export
model_code.jags_analysis <- function (object, ...) {
  model_code(as.jags_model(object), ...)
}

"model_code<-.jagr_model" <- function (object, value) {
  
  assert_that(jg_check(value))
  
  object$model_code <- value
  
  return (object)
}

#' @method model_code<- jags_model
#' @export
"model_code<-.jags_model" <- function (object, value) {
  
  if(is.list(value) && length(value) != nmodels(object))
    stop("if value is a list it must be the same length as the number of models in object")
  
  if(is.list(value))
    names(value) <- NULL
  
  models <- models(object)
  
  for (i in 1:length(models)) {
    if(!is.list(value)) {
      model_code(models[[i]]) <- value
    } else
      model_code(models[[i]]) <- value[[i]]
  }
  
  models(object) <- models
  object
}