This function samples from a specific model inside a PMP object.
sample_multinomial(PMP_object, shrink = TRUE, data, which = "MPM",
discreteSurv = TRUE)formula of the model including all potential variables
should the coefficients be shrunken towards their prior mean?
the (training) data frame with all the information
which model should be sampled from? either an integer, 'MPM' or 'MAP'
Boolean variable telling us whether a 'simple' multinomial regression is looked for or if the goal is a discrete survival-time model for multiple modes of failure is needed.
returns an object with the model coefficients and supplementary information