Running examples in ‘COMIX-Ex.R’ failed
The error most likely occurred in:

> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: effectiveSampleSize
> ### Title: This function creates an object that summarizes the effective
> ###   sample size for the parameters of the model.
> ### Aliases: effectiveSampleSize
> 
> ### ** Examples
> 
> library(COMIX)
> # Number of observations for each sample (row) and cluster (column):
> njk <- 
+   matrix(
+     c(
+       150, 300,
+       250, 200
+     ),
+     nrow = 2,
+     byrow = TRUE
+   )
> 
> # Dimension of data:
> p <- 3
> 
> # Scale and skew parameters for first cluster:
> Sigma1 <- matrix(0.5, nrow = p, ncol = p) + diag(0.5, nrow = p)
> alpha1 <- rep(0, p)
> alpha1[1] <- -5
> # location parameter for first cluster in first sample:
> xi11 <- rep(0, p)
> # location parameter for first cluster in second sample (aligned with first):
> xi21 <- rep(0, p)
> 
> # Scale and skew parameters for second cluster:
> Sigma2 <- matrix(-1/3, nrow = p, ncol = p) + diag(1 + 1/3, nrow = p)
> alpha2 <- rep(0, p)
> alpha2[2] <- 5
> # location parameter for second cluster in first sample:
> xi12 <- rep(3, p)
> # location parameter for second cluster in second sample (misaligned with first):
> xi22 <- rep(4, p)
> 
> # Sample data:
> set.seed(1)
> Y <- 
+   rbind(
+     sn::rmsn(njk[1, 1], xi = xi11, Omega = Sigma1, alpha = alpha1),
+     sn::rmsn(njk[1, 2], xi = xi12, Omega = Sigma2, alpha = alpha2),
+     sn::rmsn(njk[2, 1], xi = xi21, Omega = Sigma1, alpha = alpha1),
+     sn::rmsn(njk[2, 2], xi = xi22, Omega = Sigma2, alpha = alpha2)
+   )
> 
> C <- c(rep(1, rowSums(njk)[1]), rep(2, rowSums(njk)[2]))
> 
> prior <- list(zeta = 1, K = 10)
> pmc <- list(naprt = 5, nburn = 200, nsave = 200) # Reasonable usage
> pmc <- list(naprt = 5, nburn = 2, nsave = 5) # Minimal usage for documentation
> # Fit the model:
> res <- comix(Y, C, pmc = pmc, prior = prior)
initializing all particles...
Done
> 
> # Relabel to resolve potential label switching issues:
> res_relab <- relabelChain(res)
K=10
T=5
N=900
J=2
p=3
> effssz <- effectiveSampleSize(res_relab, "w")
Error in `summarize()`:
ℹ In argument: `Freq = table(.data$k)/sum(table(.data$k))`.
ℹ In group 1: `j = 1`.
Caused by error:
! `Freq` must be size 1, not 9.
ℹ To return more or less than 1 row per group, use `reframe()`.
Backtrace:
     ▆
  1. ├─COMIX::effectiveSampleSize(res_relab, "w")
  2. │ └─COMIX::tidyChain(res, params)
  3. │   └─... %>% select(-.data$Freq)
  4. ├─dplyr::select(., -.data$Freq)
  5. ├─dplyr::mutate(., x = (P + 1)/2, k = names(.data$Freq), frq_t = as.numeric(.data$Freq))
  6. ├─dplyr::summarize(., Freq = table(.data$k)/sum(table(.data$k)))
  7. ├─dplyr:::summarise.grouped_df(., Freq = table(.data$k)/sum(table(.data$k)))
  8. │ └─dplyr:::summarise_cols(.data, dplyr_quosures(...), by, "summarise")
  9. │   └─base::withCallingHandlers(...)
 10. ├─dplyr:::dplyr_internal_error(...)
 11. │ └─rlang::abort(class = c(class, "dplyr:::internal_error"), dplyr_error_data = data)
 12. │   └─rlang:::signal_abort(cnd, .file)
 13. │     └─base::signalCondition(cnd)
 14. └─dplyr (local) `<fn>`(`<dpl:::__>`)
 15.   └─dplyr (local) handler(cnd)
 16.     └─rlang::abort(message, class = error_class, parent = parent, call = error_call)
Execution halted