# 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