Estimate vaccine efficacy (VE) using immunogenicity data.
    The inclusion of immunogenicity data in regression models can increase precision in VE. 
    The methods are described in the publications "Elucidating vaccine efficacy using a correlate of protection, demographics, and logistic regression" and "Improving precision of vaccine efficacy evaluation using immune correlate data in time-to-event models" by Julie Dudasova, Zdenek Valenta, and Jeffrey R. Sachs (2024).

Position adjustments for 'ggplot2' to implement "visualize as
    you randomize" principles, which can be especially useful when
    plotting experimental data.

Computes the Gaussian variational approximation of the Bayesian 
    empirical likelihood posterior. This is an implementation of the function 
    found in Yu, W., & Bondell, H. D. (2023) <doi:10.1080/01621459.2023.2169701>.

The shared random effects joint model is one of the most widely used approaches to study the associations between longitudinal biomarkers and a survival outcome and make dynamic risk predictions using the longitudinally measured biomarkers. 
    One major limitation of joint models is that they could be computationally expensive for complex models where the number of the shared random effects is large.  
    This package can be used to fit complex multivariate joint models using our newly developed algorithm Jieqi Tu and Jiehuan Sun (2023) <doi:10.1002/sim.9619>, which is based on Gaussian variational approximate inference and is computationally efficient.

Fit and simulate latent position and cluster models for network data, using a fast Variational Bayes approximation developed in Salter-Townshend and Murphy (2013) <doi:10.1016/j.csda.2012.08.004>.

A Variational Bayesian algorithm for high-dimensional multi-source 
    heterogeneous linear models. More details have been written up in a paper 
    submitted to the journal Statistics in Medicine, and the details of variational
    Bayesian methods can be found in Ray and Szabo (2021) <doi:10.1080/01621459.2020.1847121>.
    It simultaneously performs parameter estimation and variable selection. The 
    algorithm supports two model settings: (1) local models, where variable selection
    is only applied to homogeneous coefficients, and (2) global models, where variable
    selection is also performed on heterogeneous coefficients. Two forms of 
    Spike-and-Slab priors are available: the Laplace distribution and the Gaussian 
    distribution as the Slab component.