quantile.BLC {BayesMortalityPlus}

BLC: Sample quantiles

Description

Compute the quantiles based on the resulting chains from a fitted BLC model.

Usage

## S3 method for class 'BLC'
quantile(x, q, name, ...)

Arguments

x

A BLC object, result of a call to blc() function.

q

A real number that represents the probability of the quantiles.

name

A character with a parameter name of the blc model that should be returned. It can be one of these: "alpha", "beta", "kappa", "phiv", "theta", "phiw".

...

Further arguments passed to or from other methods.

Value

A data.frame with the quantiles of the selected parameter.

See Also

quantile.PredBLC() for PredBLC method.

Examples

## Importing log-mortality data from Portugal:
data(PT)
Y <- PT

## Fitting the model
fit = blc(Y = Y, M = 100, bn = 20)

## Parameters' median and quantiles 0.05, 0.95
quantile(fit, c(0.05, 0.5, 0.95), "alpha")
quantile(fit, c(0.05, 0.5, 0.95), "beta")
quantile(fit, c(0.05, 0.5, 0.95), "kappa")
quantile(fit, c(0.05, 0.5, 0.95), "phiv") ## random error precision
quantile(fit, c(0.05, 0.5, 0.95), "theta") ## drift parameter
quantile(fit, c(0.05, 0.5, 0.95), "phiw")