#定义 simmr_in <- simmr_load( mixtures = mix, source_names = s_names, source_means = s_means, source_sds = s_sds, correction_means = c_means, correction_sds = c_sds, concentration_means = conc ) #运行 #1.马尔可夫链蒙特运行模型的代码 carlo（mcmc） simmr_out <- simmr_mcmc(simmr_in) #检查算法收敛（2不需要检查 prior_viz(simmr_out) #检查模型拟合 post_pred <- posterior_predictive(simmr_out) #比较来源 compare_sources(simmr_out, source_names = c( "gw", "s1", "s2", "s3", "s4", "s5", "s6" ) )