From cc6fde7bec784ec5c44223749e2e6d8043b391d9 Mon Sep 17 00:00:00 2001 From: Nicholas Clark Date: Wed, 21 Jun 2023 15:20:05 +1000 Subject: [PATCH] fix the damn bins! --- R/ppc.mvgam.R | 10 ++++++---- 1 file changed, 6 insertions(+), 4 deletions(-) diff --git a/R/ppc.mvgam.R b/R/ppc.mvgam.R index 942e632a..d1cb954d 100644 --- a/R/ppc.mvgam.R +++ b/R/ppc.mvgam.R @@ -564,12 +564,14 @@ ppc.mvgam = function(object, newdata, data_test, series = 1, type = 'hist', } breaks <- seq(xlim[1], xlim[2], length.out = n_bins) - breaks <- sort(c(min(c(truths, as.vector(preds))), - breaks, - max(c(truths, as.vector(preds))))) + truths <- truths[truths<=xlim[2]] + truths <- truths[truths>=xlim[1]] + preds <- preds[preds<=xlim[2]] + preds <- preds[preds>=xlim[1]] ylim <- c(0, max(c(max(hist(truths, breaks = breaks, plot = F)$density, na.rm = TRUE), - max(hist(preds, breaks = breaks, plot = F)$density, na.rm = TRUE)))) + max(hist(preds, breaks = breaks, plot = F)$density, na.rm = TRUE)), + na.rm = TRUE)) if(missing(xlab)){ xlab <- paste0('Count')