| # Bonus! Sample counts from predictionsP, take 95% HDPI | |
| hdpiPoisP <- apply(predictionsP, 2, HPDI, prob = .95) | |
| meanPoisP <- colMeans(predictionsP) | |
| plot(rangeA, meanPoisP, type = "l", ylim = c(0, 15), | |
| yaxp = c(0, 15, 5), xlab = "Min Age (standardized)", | |
| ylab = expression(paste(lambda, " / no. eggs laid"))) | |
| shade(hdpiPoisP, rangeA) | |
| poisample <- sapply(1:100, function(k){ | |
| rpois(nrow(predictionsP), predictionsP[,k]) | |
| }) | |
| hdpiSample <- apply(poisample, 2, HPDI, prob = .95) | |
| shade(hdpiSample, rangeA) |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment