flat_full_bayesian.Rmd: fix ddl parameter problem for Wishart distribution

dev/flat_full_bayesian.Rmd

@@ -210,7 +210,9 @@ EM_init_params <- function(
   params$mu <- list()
   params$sigma <- list()
   params$sigma_prior <- list()
-  params$alpha <- rep((sigma_weight + ncol(x)) / 2, length(params$gamma))
+  params$alpha <- (sigma_weight + ncol(x)) / 2
+  params$alpha <- max(c(params$alpha, ncol(x)))
+  params$alpha <- rep(params$alpha, length(params$gamma))
   for (k in 1:length(params$gamma)) {
     params$beta[[k]] <- sigma_weight / 2 * cov(x)
     params$sigma_prior[[k]] <- retry::retry(MCMCpack::riwish(
@@ -647,7 +649,6 @@ comparison <- data %>%
     dplyr::select(count_m, count_f) %>%
     as.matrix() %>% 
     compute_tpm() %>%
-    # log() %>%
     compare_models_constaint(nboot = 10, bootsize = 1000, core = 1)
 comparison %>% 
   ggplot(aes(x = name, y = loglik, fill = name)) +