@@ -593,22 +593,25 @@ mcmc_intervals_data <- function(x,
593593 prob = 0.5 ,
594594 prob_outer = 0.9 ,
595595 point_est = c(" median" , " mean" , " none" ),
596- rhat = numeric (),
597- .data_long = NULL ) {
596+ rhat = numeric ()) {
598597 check_ignored_arguments(... )
599598 probs <- check_interval_widths(prob , prob_outer )
600599 prob <- probs [1 ]
601600 prob_outer <- probs [2 ]
602601
603- if (is.null(.data_long )) {
604- x <- prepare_mcmc_array(x , pars , regex_pars , transformations )
605- x <- merge_chains(x )
606- data_long <- melt_mcmc(x ) %> %
607- dplyr :: as_tibble() %> %
608- rlang :: set_names(tolower )
609- } else {
610- data_long <- .data_long
611- }
602+ data_long <- melt_mcmc(
603+ merge_chains(prepare_mcmc_array(x , pars , regex_pars , transformations ))
604+ ) %> %
605+ dplyr :: as_tibble() %> %
606+ rlang :: set_names(tolower )
607+
608+ compute_intervals(data_long , prob , prob_outer , point_est , rhat )
609+ }
610+
611+ # Internal helper shared by mcmc_intervals_data() and mcmc_areas_data()
612+ compute_intervals <- function (data_long , prob , prob_outer ,
613+ point_est = c(" median" , " mean" , " none" ),
614+ rhat = numeric ()) {
612615
613616 probs <- c(0.5 - prob_outer / 2 ,
614617 0.5 - prob / 2 ,
@@ -702,9 +705,9 @@ mcmc_areas_data <- function(x,
702705 dplyr :: as_tibble() %> %
703706 rlang :: set_names(tolower )
704707
705- intervals <- mcmc_intervals_data( x , prob = probs [1 ], prob_outer = probs [ 2 ],
706- point_est = temp_point_est , rhat = rhat ,
707- .data_long = data_long )
708+ intervals <- compute_intervals( data_long , prob = probs [1 ],
709+ prob_outer = probs [ 2 ] ,
710+ point_est = temp_point_est , rhat = rhat )
708711
709712 # Compute the density intervals
710713 data_inner <- data_long %> %
0 commit comments