Package: mcmcsae 0.8.0

mcmcsae: Markov Chain Monte Carlo Small Area Estimation

Fit multi-level models with possibly correlated random effects using Markov Chain Monte Carlo simulation. Such models allow smoothing over space and time and are useful in, for example, small area estimation.

Authors:Harm Jan Boonstra [aut, cre], Grzegorz Baltissen [ctb]

mcmcsae_0.8.0.tar.gz
mcmcsae_0.8.0.zip(r-4.7)mcmcsae_0.8.0.zip(r-4.6)mcmcsae_0.8.0.zip(r-4.5)
mcmcsae_0.8.0.tgz(r-4.6-x86_64)mcmcsae_0.8.0.tgz(r-4.6-arm64)mcmcsae_0.8.0.tgz(r-4.5-x86_64)mcmcsae_0.8.0.tgz(r-4.5-arm64)
mcmcsae_0.8.0.tar.gz(r-4.7-arm64)mcmcsae_0.8.0.tar.gz(r-4.7-x86_64)mcmcsae_0.8.0.tar.gz(r-4.6-arm64)mcmcsae_0.8.0.tar.gz(r-4.6-x86_64)
mcmcsae_0.8.0.tgz(r-4.6-emscripten)
manual.pdf |manual.html✨
card.svg |card.png
mcmcsae/json (API)
NEWS

# Install 'mcmcsae' in R:
install.packages('mcmcsae', repos = c('https://hjboonstra.r-universe.dev', 'https://cloud.r-project.org'))
Uses libs:
  • c++– GNU Standard C++ Library v3

On CRAN:

Conda:

This package does not link to any Github/Gitlab/R-forge repository. No issue tracker or development information is available.

cpp

2.78 score 8 scripts 561 downloads 86 exports 26 dependencies

Last updated from:8aea6b0e74. Checks:13 OK. Indexed: yes.

TargetResultTimeFilesSyslog
linux-devel-arm64OK297
linux-devel-x86_64OK317
source / vignettesOK300
linux-release-arm64OK283
linux-release-x86_64OK349
macos-release-arm64OK251
macos-release-x86_64OK485
macos-oldrel-arm64OK228
macos-oldrel-x86_64OK580
windows-develOK325
windows-releaseOK452
windows-oldrelOK319
wasm-releaseOK157

Exports:%m*v%acceptance_ratesaggrMatrixAR1brtCGCG_controlchol_controlcombine_chainscombine_iterscompute_DICcompute_WAICcomputeDesignMatrixcreate_cMVN_samplercreate_samplercreate_TMVN_samplercrossprod_mvcustomf_binomialf_gammaf_gaussianf_gaussian_gammaf_multinomialf_negbinomialf_poissongengen_controlgenerate_dataget_drawget_meansget_sdsglregGMRF_structureiidlabels<-m_directm_Gibbsm_HMCm_HMCZigZagm_softTMVNmaximize_log_lh_pmc_offsetmcmcsae_exampleMCMCsimmecmodel_matrixn_chainsn_drawsn_effn_varsnegbin_controlpar_namesplot_coefpoisson_controlpr_betapr_exppr_fixedpr_gammapr_gigpr_invchisqpr_invwishartpr_MLiGpr_normalpr_truncnormalpr_unifR_hatread_drawsregRW1RW2ssampler_controlSBC_testseasonset_constraintsset_MHsetup_clustersim_marg_varspatialsplinesstop_clusterto_draws_arrayto_mcmctransform_dcvfacvreg

Dependencies:abindbackportscheckmateclicollapsedistributionalgenericsGIGrvggluelatticelifecycleloomagrittrMatrixmatrixStatsnumDerivpillarpkgconfigposteriorRcppRcppEigenrlangtensorAtibbleutf8vctrs

Basic area-level model

Rendered fromarea_level.Rmdusingknitr::rmarkdownon May 10 2026.

Last update: 2025-06-04
Started: 2020-09-01

Basic unit-level models

Rendered fromunit_level.Rmdusingknitr::rmarkdownon May 10 2026.

Last update: 2025-06-04
Started: 2020-09-01

Linear regression, prediction, and survey weighting

Rendered fromlinear_weighting.Rmdusingknitr::rmarkdownon May 10 2026.

Last update: 2025-06-04
Started: 2020-09-01

Readme and manuals

Help Manual

Help pageTopics
Markov Chain Monte Carlo Small Area Estimationmcmcsae-package mcmcsae
Return Metropolis-Hastings acceptance ratesacceptance_rates
Utility function to construct a sparse aggregation matrix from a factoraggrMatrix
Create a model component object for a BART (Bayesian Additive Regression Trees) component in the linear predictorbrt
Set options for the conjugate gradient (CG) samplerCG_control
Set options for Cholesky decompositionchol_control
Combine multiple mcdraws objects into a single one by combining their chainscombine_chains
Combine multiple mcdraws objects into a single one by combining their drawscombine_iters
Compute a list of design matrices for all terms in a model formula, or based on a sampler environmentcomputeDesignMatrix
Correlation factor structures in generic model componentsAR1 correlation custom iid RW1 RW2 season spatial splines
Set up a function for direct sampling from a constrained multivariate normal distributioncreate_cMVN_sampler
Create a sampler objectcreate_sampler
Set up a sampler object for sampling from a possibly truncated and degenerate multivariate normal distributioncreate_TMVN_sampler
Specify a binomial sampling distributionf_binomial
Specify a Gamma sampling distributionf_gamma
Specify a Gaussian sampling distributionf_gaussian
Specify a Gaussian-Gamma sampling distributionf_gaussian_gamma
Specify a multinomial sampling distributionf_multinomial
Specify a negative binomial sampling distributionf_negbinomial
Specify a Poisson sampling distributionf_poisson
Create a model component object for a generic random effects component in the linear predictorgen
Set computational options for the sampling algorithms used for a 'gen' model componentgen_control
Generate a data vector according to a modelgenerate_data
Extract a list of parameter values for a single drawget_draw
Create a model object for group-level regression effects within a generic random effects component.glreg
Set up a GMRF structure for a generic model componentGMRF_structure
Get and set the variable labels of a draws component object for a vector-valued parameterlabels labels.dc labels<-
Fast matrix-vector multiplications%m*v% crossprod_mv matrix-vector
Maximise the log-likelihood or log-posterior as defined by a sampler closuremaximize_log_lh_p
Create a model component object for an offset, i.e. fixed, non-parametrised term in the linear predictormc_offset
Compute MCMC diagnostic measuresMCMC-diagnostics n_eff R_hat
Convert a draws component object to another formatas.array.dc as.matrix.dc MCMC-object-conversion to_draws_array to_mcmc
Generate artificial data according to an additive spatio-temporal modelmcmcsae_example
Run a Markov Chain Monte Carlo simulationMCMCsim
Create a model component object for a regression (fixed effects) component in the linear predictor with measurement errors in quantitative covariatesmec
Compute possibly sparse model matrixmodel_matrix
Compute DIC, WAIC and leave-one-out cross-validation model measurescompute_DIC compute_WAIC loo.mcdraws model-information-criteria waic.mcdraws
Get the number of chains, samples per chain or the number of variables in a simulation objectn_chains n_chains-n_draws-n_vars n_draws n_vars
Set computational options for the sampling algorithmsnegbin_control
Get the parameter names from an mcdraws objectpar_names
Plot a set of model coefficients or predictions with uncertainty intervals based on summaries of simulation results or other objects.plot_coef
Trace, density and autocorrelation plots for (parameters of a) draws component (dc) objectplot.dc
Trace, density and autocorrelation plotsplot.mcdraws
Set computational options for the sampling algorithmspoisson_control
Get means or standard deviations of parameters from the MCMC output in an mcdraws objectget_means get_sds posterior-moments
Create an object representing beta prior distributionspr_beta
Create an object representing exponential prior distributionspr_exp
Create an object representing a degenerate prior fixing a parameter (vector) to a fixed valuepr_fixed
Create an object representing gamma prior distributionspr_gamma
Create an object representing Generalised Inverse Gaussian (GIG) prior distributionspr_gig
Create an object representing inverse chi-squared priors with possibly modelled degrees of freedom and scale parameterspr_invchisq
Create an object representing an inverse Wishart prior, possibly with modelled scale matrixpr_invwishart
Create an object representing a Multivariate Log inverse Gamma (MLiG) prior distributionpr_MLiG
Create an object representing a possibly multivariate normal prior distributionpr_normal
Create an object representing truncated normal prior distributionspr_truncnormal
Create an object representing uniform prior distributionspr_unif
Generate draws from the predictive distributionpredict.mcdraws
Display a summary of a 'dc' objectprint.dc_summary
Print a summary of MCMC simulation resultsprint.mcdraws_summary
Read MCMC draws from a fileread_draws
Specify a regression (fixed effects) component in the linear predictorreg
Extract draws of fitted values or residuals from an mcdraws objectfitted.mcdraws residuals-fitted-values residuals.mcdraws
Specify a smooth term component of the linear predictors
Set computational options for the sampling algorithmssampler_control
Simulation based calibrationSBC_test
Set up a system of linear equality and/or inequality constraintsset_constraints
Set options for Metropolis-Hastings samplingset_MH
Set up a cluster for parallel computingsetup_cluster
Compute a Monte Carlo estimate of the marginal variances of a (I)GMRFsim_marg_var
Stop a clusterstop_cluster
Select a subset of chains, samples and parameters from a draws component (dc) objectsubset.dc
Summarise a draws component (dc) objectsummary.dc
Summarise an mcdraws objectsummary.mcdraws
Functions for specifying the method and corresponding options for sampling from a possibly truncated and degenerate multivariate normal distributionm_direct m_Gibbs m_HMC m_HMCZigZag m_softTMVN TMVN-methods
Transform one or more draws component objects into a new one by applying a functiontransform_dc
Create a model component object for a variance factor component in the variance function of a gaussian sampling distributionvfac
Create a model component object for a regression component in the variance function of a gaussian sampling distributionvreg
Extract weights from an mcdraws objectweights.mcdraws