fix package

caravagnalab · Nov 14, 2024 · 77c92a8 · 77c92a8
1 parent 1487914
commit 77c92a8
Showing 1 changed file with 2 additions and 5 deletions.
diff --git a/R/fit_h.R b/R/fit_h.R
@@ -90,13 +90,10 @@ fit_h = function(x, max_attempts=2, INIT=TRUE, tolerance = 0.01, possible_k = c(
     tau_and_w_draws <- res$draws(variables = vars)
     tau_and_w_summary <- res$summary(variables = vars)
 
-    # implement it differenctly without stanfit?
-    # stanfit <- rstan::read_stan_csv(res$output_files())
     # Check log likelihood values  POSSO ESTRARRE LA LOG LIK DIRETTAMENTE DAL MODELLO E NON DALLE GENERATED QUANTITIES?
 
-    # log_lik_matrix <- rstan::extract_log_lik(stanfit, parameter_name = "log_lik", merge_chains = TRUE) # merge chains potrebbe non servire  #getter?
-    log_lik <- res$draws(variables = "log_lik")
-    log_lik_matrix <- res$draws(variables = "log_lik_matrix")
+    # log_lik_contributions <- res$draws(variables = "log_lik_matrix")
+    log_lik_matrix <- res$draws(variables = "log_lik")
     log_lik_matrix_list[[K]] <- log_lik_matrix
     result_single <- list(draws = tau_and_w_draws, summary = tau_and_w_summary)
     draws_and_summary[[K]] = result_single