posteriorSamples | R Documentation |
Draw samples from the empirical Bayes posterior predictive distribution derived from unmarked models or ranef objects
## S4 method for signature 'unmarkedRanef'
posteriorSamples(object, nsims=100, ...)
## S4 method for signature 'unmarkedFit'
posteriorSamples(object, nsims=100, ...)
object |
An object inheriting class |
nsims |
Number of draws to make from the posterior predictive distribution |
... |
Other arguments |
unmarkedPostSamples
object containing the draws from the
posterior predictive distribution. The draws are in the @samples
slot.
Ken Kellner contact@kenkellner.com
ranef
,
predict
# Simulate data under N-mixture model
set.seed(4564)
R <- 20
J <- 5
N <- rpois(R, 10)
y <- matrix(NA, R, J)
y[] <- rbinom(R*J, N, 0.5)
# Fit model
umf <- unmarkedFramePCount(y=y)
fm <- pcount(~1 ~1, umf, K=50)
# Estimates of conditional abundance distribution at each site
(re <- ranef(fm))
#Draw from the posterior predictive distribution
(ppd <- posteriorSamples(re, nsims=100))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.