BnNsir function

Bayesian analysis of n >= 2 Normal means with standard improper reference priors.

Several exploratory plots are produced, after which this function calculates and plots the posterior densities of the treatment means and their differences. Pooled or separate variances can be specified. Note that this function uses Standard Improper Reference (SIR) priors.

BnNsir(formula, data, var.equal = TRUE, alpha = 0.05, plotit = TRUE,
 r = 10000)

Arguments

• formula: the usual formula interface: response ~ factor
• data: a data.frame containing the response and the factor variables
• var.equal: should a pooled variance be used? Specify var.equal = FALSE if you want separate variances to be fitted
• alpha: 1 - level of credibility, so that for alpha = 0.05 (the default) credible intervals will have 95% credibility
• plotit: are plots desired?
• r: the number of samples of the posterior that should be taken

Returns

none returned: the function is used for the plots and the printed information it produces

References

van Hulst, R. 2018. Evaluating Scientific Evidence. ms.

Author(s)

Robert van Hulst

See Also

B2Nsir

Examples

data(PlantGrowth)
BnNsir(weight ~ group, PlantGrowth)