lpost.gamma {ordgam}R Documentation

Posterior density function for the non-penalized parameters in an ordgam model

Description

Posterior density function for the non-penalized parameters in an ordgam model

Usage

lpost.gamma(model)

Arguments

model

An ordgam.object

Value

Log joint posterior density function for the non-penalized regression parameters.

Author(s)

Philippe Lambert p.lambert@uliege.be

References

Lambert, P. and Gressani, 0. (2023) Penalty parameter selection and asymmetry corrections to Laplace approximations in Bayesian P-splines models. Statistical Modelling. <doi:10.1177/1471082X231181173>. Preprint: <arXiv:2210.01668>.

See Also

ordgam.

Examples

library(ordgam)
data(freehmsData)
mod = ordgam(freehms ~ s(eduyrs) + s(age), data=freehmsData, descending=TRUE)
print(mod$theta) ## Model regression parameters
gam.hat = mod$theta[1:4] ## Non-penalized parameter estimates
ordgam::lpost.gamma(mod)(gam.hat)
[Package ordgam version 0.9.1 Index]