R/pi_MCSim_dirichlet.R

In MSRDT: Multi-State Reliability Demonstration Tests (MSRDT)

#' @title Dirichlet Prior Simulation for Multi-state RDT
#'
#' @description Define the simulation function to generate failure probability with Dirichlet prior distributions as conjugate prior to multinomial distributions (for multi-state RDT).
#'
#' @param M Simulation sample size
#' @param seed Random seed for random sample
#' @param par Parameters for dirichlet distribution
#' @return Vector of failure probability sample
#' @examples
#' pi <- pi_MCSim_dirichlet(M = 1000, seed = 10, par = c(1, 1, 1))
#' @export
#' @importFrom gtools rdirichlet
#' @family Prior distribution generation functions
#' @seealso \code{\link{pi_MCSim_beta}}

pi_MCSim_dirichlet <- function(M, seed, par){
  #requireNamespace("gtools")
  set.seed(seed)
  return(rdirichlet(M, par))
}