pn.ICred {PowerNormal}

Credibility interval for alpha (PN distribution)

Description

Credibility interval for the power parameter of PN distribution

Usage

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

Arguments

x

the response vector

p

credibility level

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.ICred(x, 0.95)

pn.ICred(x, 0.95, prior = 'Uniform')

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