R/statpsych1.R
ci.bayes.normal.Rd
Computes an approximate Bayesian credible interval for a normal prior distribution. This function can be used with any parameter estimator (e.g., mean, mean difference, linear contrast of means, slope coefficient) that has a normal sampling distribution. The mean and standard deviation of the posterior normal distribution are also reported.
ci.bayes.normal(alpha, prior.mean, prior.sd, est, se)
alpha level for 1-alpha credibility interval
mean of prior Normal distribution
standard deviation of prior Normal distribution
sample estimate
standard error of sample estimate
Returns a 1-row matrix. The columns are:
Posterior mean - posterior mean of Normal distributoin
Posterior SD - posterior standard deviation of Normal distributoin
LL - lower limit of the credible interval
UL - upper limit of the credible interval
Gelman ACJBSHS, Rubin DB (2004). Bayesian Data Analysis, 2nd edition. Chapman & Hall.
ci.bayes.normal(.05, 30, 2, 24.5, 0.577)
#> Posterior mean Posterior SD LL UL
#> [1,] 24.9226 0.5543895 23.83602 26.00919
# Should return:
# Posterior mean Posterior SD LL UL
# [1,] 24.9226 0.5543895 23.83602 26.00919