Summarize the reporting group assignment probabilities for each individual in the mixture. See vignette for details.
Examples
# set up input data and run multistage model
msgsi_dat <-
prep_msgsi_data(mixture_data = mix,
baseline1_data = base_templin, baseline2_data = base_yukon,
pop1_info = templin_pops211, pop2_info = yukon_pops50, sub_group = 3:5)
#> Compiling input data, may take a minute or two...
#> Time difference of 10.76239 secs
# run model
msgsi_out <- msgsi_mdl(msgsi_dat, nreps = 25, nburn = 15, thin = 1, nchains = 1)
#> Running model... and good things come to Femme Queen Vogue!
#> Time difference of 1.617084 secs
#> June-18-2026 20:31
# individual assignment summary
ind_iden <- indiv_assign(msgsi_out, msgsi_dat)
