Extends 'multiverse' package
    (Sarma A., Kale A., Moon M., Taback N., Chevalier F., Hullman J., Kay M., 2021)
    <doi:10.31219/osf.io/yfbwm>, which allows users perform to create explorable
    multiverse analysis in R. This extension provides an additional level of
    abstraction to the 'multiverse' package with the aim of creating
    user friendly syntax to researchers, educators, and students in statistics.
    The 'mverse' syntax is designed to allow piping and takes hints from
    the 'tidyverse' grammar. The package allows users to define and inspect
    multiverse analysis using familiar syntax in R.

Multivariate estimation and testing, currently a package for testing parametric data. To deal with parametric data, various multivariate normality tests and outlier detection are performed and visualized using the 'ggplot2' package. Homogeneity tests for covariance matrices are also possible, as well as the Hotelling's T-square test and the multivariate analysis of variance test. We are exploring additional tests and visualization techniques, such as profile analysis and randomized complete block design, to be made available in the future and making them easily accessible to users.

Fit Bayesian Dynamic Generalized Additive Models to sets of time series. Users can build dynamic nonlinear State-Space models that can incorporate semiparametric effects in observation and process components, using a wide range of observation families. Estimation is performed using Markov Chain Monte Carlo with Hamiltonian Monte Carlo in the software 'Stan'. References: Clark & Wells (2022) <doi:10.1111/2041-210X.13974>.

Generates multivariate subgaussian stable probabilities using the QRSVN
    algorithm as detailed in Genz and Bretz (2002) <DOI:10.1198/106186002394>
    but by sampling positive stable variates not chi/sqrt(nu).

Maximum likelihood estimates are obtained via an EM algorithm with either a first-order or a fully exponential Laplace approximation as documented by Broatch and Karl (2018) <doi:10.48550/arXiv.1710.05284>,
    Karl, Yang, and Lohr (2014) <doi:10.1016/j.csda.2013.11.019>, and by 
	Karl (2012) <doi:10.1515/1559-0410.1471>. Karl and Zimmerman <doi:10.1016/j.jspi.2020.06.004> use this package to illustrate how the home field effect estimator from a mixed model can be biased under nonrandom scheduling.