From 10bbc366f10fe19409b625642cf67bb20f223ca1 Mon Sep 17 00:00:00 2001
From: Laurent Modolo <laurent.modolo@ens-lyon.fr>
Date: Wed, 30 Aug 2023 09:00:58 +0200
Subject: [PATCH] flat_full_bayesian.Rmd: fix ddl parameter problem for Wishart
distribution
---
dev/flat_full_bayesian.Rmd | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/dev/flat_full_bayesian.Rmd b/dev/flat_full_bayesian.Rmd
index 4837de5..562293d 100644
--- a/dev/flat_full_bayesian.Rmd
+++ b/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)) +
--
GitLab