pn.bayes {PowerNormal}

Fit univariate PN distribution (Bayesian)

Description

Return the posterior mean, median and variance of power parameter for PN distribution

Usage

pn.bayes(x,prior= "Jeffreys", shape_0 = NULL, rate_0 = NULL)

Arguments

x

the response vector

prior

the prior distribution of power parameter: "Jeffreys" (default), "Uniform" and "Gamma"

shape_0, rate_0

shape and rate hyperparameters of the gamma distribution.

References

Agamez-Montavo, G. (2017). Modelos de mistura finita usando a classe de distribuicoes alpha potencia. Thesis (Doctoral), University of Sao Paulo.

Examples

x <- rpn(100, 25)

pn.bayes(x)

pn.bayes(x, prior = 'Uniform')

pn.bayes(x, prior = 'Gamma', 1/100, 1/100)