2021-12-08から1日間の記事一覧
結論 この関数を使おう ※stanファイルの対数尤度の変数名がlog_likとして推定している事を仮定 waic2 <- function(fit){ # log_lik <- rstan::extract(fit)$log_lik # rtan用 log_lik <- fit$draws("log_lik") %>% as_draws_df() %>% select(contains("log_l…
結論 この関数を使おう ※stanファイルの対数尤度の変数名がlog_likとして推定している事を仮定 waic2 <- function(fit){ # log_lik <- rstan::extract(fit)$log_lik # rtan用 log_lik <- fit$draws("log_lik") %>% as_draws_df() %>% select(contains("log_l…