post2 <- stan("ebola2.stan", refresh = 0, # suppresses intermediate output
data = list(exposed = 7, survived = 5, alpha = 0.25, low = 0.4,
median = 0.55, high = 0.7))
post2
## Inference for Stan model: ebola2.
## 4 chains, each with iter=2000; warmup=1000; thin=1;
## post-warmup draws per chain=1000, total post-warmup draws=4000.
##
## mean se_mean sd 2.5% 25% 50% 75% 97.5% n_eff Rhat
## p 0.63 0.01 0.21 0.19 0.49 0.67 0.80 0.94 1433 1
## pi 0.64 0.00 0.14 0.36 0.54 0.65 0.74 0.87 1431 1
## lp__ -3.27 0.02 0.80 -5.39 -3.46 -2.97 -2.78 -2.72 1445 1
##
## Samples were drawn using NUTS(diag_e) at Sat Sep 8 07:42:49 2018.
## For each parameter, n_eff is a crude measure of effective sample size,
## and Rhat is the potential scale reduction factor on split chains (at
## convergence, Rhat=1).