Fits bootstrap with univariate spatial regression models using Bootstrap for Rapid Inference on Spatial Covariances (BRISC) for large datasets using nearest neighbor Gaussian processes detailed in Saha and Datta (2018) <doi:10.1002/sta4.184>.

Quantitative methods for benefit-risk analysis help to condense
    complex decisions into a univariate metric describing the overall benefit
    relative to risk.  One approach is to use the multi-criteria decision
    analysis framework (MCDA), as in Mussen, Salek, and Walker
    (2007) <doi:10.1002/pds.1435>.  Bayesian benefit-risk
    analysis incorporates uncertainty through posterior distributions which are
    inputs to the benefit-risk framework.  The brisk package provides functions
    to assist with Bayesian benefit-risk analyses, such as MCDA.
    Users input posterior samples, utility functions, weights, and the package
    outputs quantitative benefit-risk scores.  The posterior of the benefit-risk
    scores for each group can be compared.  Some plotting capabilities are also
    included.

Implementation of the record linkage methodology proposed by Sadinle (2017) <doi:10.1080/01621459.2016.1148612>.  It handles the bipartite record linkage problem, where two duplicate-free datafiles are to be merged.

Provides two main functions, il() and fil(). The il() function implements the EM algorithm developed by Ibrahim and Lipsitz (1996) <DOI:10.2307/2533068> to estimate the parameters of a logistic regression model with the missing response when the missing data mechanism is nonignorable. The fil() function implements the algorithm proposed by Maity et. al. (2017+) <https://github.com/arnabkrmaity/brlrmr> to reduce the bias produced by the method of Ibrahim and Lipsitz (1996) <DOI:10.2307/2533068>.

Fits novel models for the conditional relative risk, risk difference and odds ratio <doi:10.1080/01621459.2016.1192546>.

Fit Bayesian generalized (non-)linear multivariate multilevel models
    using 'Stan' for full Bayesian inference. A wide range of distributions
    and link functions are supported, allowing users to fit -- among others --
    linear, robust linear, count data, survival, response times, ordinal,
    zero-inflated, hurdle, and even self-defined mixture models all in a
    multilevel context. Further modeling options include both theory-driven and
    data-driven non-linear terms, auto-correlation structures, censoring and 
    truncation, meta-analytic standard errors, and quite a few more. 
    In addition, all parameters of the response distribution can be predicted 
    in order to perform distributional regression. Prior specifications are 
    flexible and explicitly encourage users to apply prior distributions that
    actually reflect their prior knowledge. Models can easily be evaluated and
    compared using several methods assessing posterior or prior predictions. 
    References: Bürkner (2017) <doi:10.18637/jss.v080.i01>; 
    Bürkner (2018) <doi:10.32614/RJ-2018-017>;
    Bürkner (2021) <doi:10.18637/jss.v100.i05>; Carpenter et al. (2017)
    <doi:10.18637/jss.v076.i01>.