Dirichlet Process Clustering with Dissimilarities

A Bayesian hierarchical model for clustering dissimilarity data using the Dirichlet process. The latent configuration of objects and the number of clusters are automatically inferred during the fitting process. The package supports multiple models which are available to detect clusters of various shapes and sizes using different covariance structures. Additional functions are included to ensure adequate model fits through prior and posterior predictive checks.

Kernel Density Estimation for Heaped and Rounded Data

In self-reported or anonymised data the user often encounters heaped data, i.e. data which are rounded (to a possibly different degree of coarseness). While this is mostly a minor problem in parametric density estimation the bias can be very large for non-parametric methods such as kernel density estimation. This package implements a partly Bayesian algorithm treating the true unknown values as additional parameters and estimates the rounding parameters to give a corrected kernel density estimate. It supports various standard bandwidth selection methods. Varying rounding probabilities (depending on the true value) and asymmetric rounding is estimable as well: Gross, M. and Rendtel, U. (2016) (). Additionally, bivariate non-parametric density estimation for rounded data, Gross, M. et al. (2016) (), as well as data aggregated on areas is supported.

Monte Carlo and Bayesian Stop-Signal Reaction Time Estimation

Estimates stop-signal reaction time (SSRT) in the stop-signal task using the integration and mean methods described by Verbruggen and colleagues (2019) . In addition to point estimates, the package provides Monte Carlo tools (nonparametric bootstrap confidence intervals, parametric ex-Gaussian simulation, minimum-trial-count and power analysis, and sensitivity analysis under violations of the horse-race assumptions) and Bayesian estimation via 'Stan', including single-subject and hierarchical ex-Gaussian horse-race models with an optional trigger-failure parameter following Matzke and colleagues (2013) , posterior inhibition functions, and posterior predictive checks. The Bayesian layer works with either the 'cmdstanr' or 'rstan' backend.

Bayesian Spectral Analysis Models using Gaussian Process Priors

Contains functions to perform Bayesian inference using a spectral analysis of Gaussian process priors. Gaussian processes are represented with a Fourier series based on cosine basis functions. Currently the package includes parametric linear models, partial linear additive models with/without shape restrictions, generalized linear additive models with/without shape restrictions, and density estimation model. To maximize computational efficiency, the actual Markov chain Monte Carlo sampling for each model is done using codes written in FORTRAN 90. This software has been developed using funding supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (no. NRF-2016R1D1A1B03932178 and no. NRF-2017R1D1A3B03035235).

Extracting and Visualizing Bayesian Graphical Models

Fit and visualize the results of a Bayesian analysis of networks commonly found in psychology. The package supports cross-sectional network models fitted using the packages 'BDgraph', 'bgms' and 'BGGM', as well as network comparison tests fitted using the packages 'bgms' and 'BBGM'. The package provides the parameter estimates, posterior inclusion probabilities, inclusion Bayes factor, and the posterior density of the parameters. In addition, for 'BDgraph' and 'bgms' it allows to assess the posterior structure space. Furthermore, the package comes with an extensive suite for visualizing results.

Bayesian Exponential Smoothing Models with Trend Modifications

An implementation of a number of Global Trend models for time series forecasting that are Bayesian generalizations and extensions of some Exponential Smoothing models. The main differences/additions include 1) nonlinear global trend, 2) Student-t error distribution, and 3) a function for the error size, so heteroscedasticity. The methods are particularly useful for short time series. When tested on the well-known M3 dataset, they are able to outperform all classical time series algorithms. The models are fitted with MCMC using the 'rstan' package.

Bayesian Augmented Control for Clinical Trials

Implements the Bayesian Augmented Control (BAC, a.k.a. Bayesian historical data borrowing) method under clinical trial setting by calling 'Just Another Gibbs Sampler' ('JAGS') software. In addition, the 'BACCT' package evaluates user-specified decision rules by computing the type-I error/power, or probability of correct go/no-go decision at interim look. The evaluation can be presented numerically or graphically. Users need to have 'JAGS' 4.0.0 or newer installed due to a compatibility issue with 'rjags' package. Currently, the package implements the BAC method for binary outcome only. Support for continuous and survival endpoints will be added in future releases. We would like to thank AbbVie's Statistical Innovation group and Clinical Statistics group for their support in developing the 'BACCT' package.

Generalized Additive Models for Location Scale and Shape

Functions for fitting the Generalized Additive Models for Location Scale and Shape introduced by Rigby and Stasinopoulos (2005), . The models use a distributional regression approach where all the parameters of the conditional distribution of the response variable are modelled using explanatory variables.

Spatial Implementation of Bayesian Networks and Mapping

Allows spatial implementation of Bayesian networks and mapping in geographical space. It makes maps of expected value (or most likely state) given known and unknown conditions, maps of uncertainty measured as coefficient of variation or Shannon index (entropy), maps of probability associated to any states of any node of the network. Some additional features are provided as well: parallel processing options, data discretization routines and function wrappers designed for users with minimal knowledge of the R language. Outputs can be exported to any common GIS format.

Bayesian Inference with Laplace Approximations and P-Splines

Laplace approximations and penalized B-splines are combined for fast Bayesian inference in latent Gaussian models. The routines can be used to fit survival models, especially proportional hazards and promotion time cure models (Gressani, O. and Lambert, P. (2018) ). The Laplace-P-spline methodology can also be implemented for inference in (generalized) additive models (Gressani, O. and Lambert, P. (2021) ). See the associated website for more information and examples.