R/nof1.prior.R

In MikeJSeo/nof1: Analysis of N of 1 Trials

nof1.prior.default <- function(prior.param){

  with(prior.param, {

  if(is.null(alpha.prior)){
    alpha.prior <- list("dnorm", 0, 0.001)
  }

  if(is.null(beta.prior)){
    if(response == "ordinal"){
      beta.prior <- list("dnorm", 0, 0.25)
    } else{
      beta.prior <- list("dnorm", 0, 0.001)    
    }
  }

  if(is.null(gamma.prior)){
    gamma.prior <- list("dnorm", 0, 0.001)
  }

  if(is.null(rho.prior)){
    rho.prior <- list("dunif", -1, 1)
  }

  if(is.null(hy.prior)){
    if(response == "binomial"){
      hy.prior <- list("dunif", 0, 5)
    } else if (response == "normal"){
      hy.prior <- list("dgamma", 0.001, 0.001)
    #  hy.prior <- list("dunif", 0, 20)
    } else{
      hy.prior <- list("dgamma", 1, 0.001)
    }

  }

  prior.param.default <- list(alpha.prior = alpha.prior, beta.prior = beta.prior, gamma.prior = gamma.prior, rho.prior = rho.prior, hy.prior = hy.prior)

  if(response == "ordinal"){
    if(is.null(dc.prior)){
      dc.prior <- list("dunif", 0, 20)
    }
    if(is.null(c1.prior)){
      c1.prior <- list("dunif", -20, 20)
    }
    prior.param.default$dc.prior <- dc.prior
    prior.param.default$c1.prior <- c1.prior
  }

  return(prior.param.default)
  })
}