Package: susieR 0.16.5

Peter Carbonetto

susieR: Sum of Single Effects Linear Regression

Implements methods for variable selection in linear regression based on the "Sum of Single Effects" (SuSiE) model, as described in Wang et al (2020) <doi:10.1101/501114> and Zou et al (2021) <doi:10.1101/2021.11.03.467167>. These methods provide simple summaries, called "Credible Sets", for accurately quantifying uncertainty in which variables should be selected. The methods are motivated by genetic fine-mapping applications, and are particularly well-suited to settings where variables are highly correlated and detectable effects are sparse. The fitting algorithm, a Bayesian analogue of stepwise selection methods called "Iterative Bayesian Stepwise Selection" (IBSS), is simple and fast, allowing the SuSiE model be fit to large data sets (thousands of samples and hundreds of thousands of variables).

Authors:Gao Wang [aut], Yuxin Zou [aut], Alexander McCreight [aut], Kaiqian Zhang [aut], William R.P. Denault [aut], Peter Carbonetto [aut, cre], Matthew Stephens [aut]

susieR_0.16.5.tar.gz
susieR_0.16.5.zip(r-4.7)susieR_0.16.5.zip(r-4.6)susieR_0.16.5.zip(r-4.5)
susieR_0.16.5.tgz(r-4.6-x86_64)susieR_0.16.5.tgz(r-4.6-arm64)susieR_0.16.5.tgz(r-4.5-x86_64)susieR_0.16.5.tgz(r-4.5-arm64)
susieR_0.16.5.tar.gz(r-4.7-arm64)susieR_0.16.5.tar.gz(r-4.7-x86_64)susieR_0.16.5.tar.gz(r-4.6-arm64)susieR_0.16.5.tar.gz(r-4.6-x86_64)
susieR_0.16.5.tgz(r-4.6-emscripten)
manual.pdf |manual.html
DESCRIPTION
card.svg |card.png
susieR/json (API)

# Install 'susieR' in R:
install.packages('susieR', repos = c('https://stephenslab.r-universe.dev', 'https://cloud.r-project.org'))

Bug tracker:https://github.com/stephenslab/susier/issues

Uses libs:
  • openblas– Optimized BLAS
  • c++– GNU Standard C++ Library v3
  • openmp– GCC OpenMP (GOMP) support library
Datasets:

On CRAN:

Conda:

openblascppopenmp

14.23 score 260 stars 5 packages 1.7k scripts 3.7k downloads 2 mentions 54 exports 31 dependencies

Last updated from:b10d47decb. Checks:8 WARNING, 2 OK, 3 ERROR. Indexed: yes.

TargetResultTimeFilesSyslog
linux-devel-arm64WARNING308
linux-devel-x86_64WARNING316
source / vignettesOK630
linux-release-arm64WARNING292
linux-release-x86_64WARNING308
macos-release-arm64WARNING243
macos-release-x86_64WARNING374
macos-oldrel-arm64WARNING294
macos-oldrel-x86_64WARNING403
windows-develERROR316
windows-releaseERROR330
windows-oldrelERROR307
wasm-releaseOK132

Exports:absolute.orderblock_coordinate_ascentcalc_zcoef.mr.ashcoef.susiecompute_marginal_bhat_shatcompute_suff_statestimate_s_rssget_cs_correlationget_objectiveget.full.posterioribss_finalizeibss_initializeis_symmetric_matrixkriging_rssmr.ashmr.ash.rsspath.orderpost_loglik_prior_hookpre_loglik_prior_hookpredict.mr.ashpredict.susieprint.summary.susieprint.summary.susie_post_outcome_configurationslot_prior_betabinomslot_prior_poissonsummary.susiesummary.susie_post_outcome_configurationsusiesusie_autosusie_get_cssusie_get_cs_attainablesusie_get_lfsrsusie_get_nitersusie_get_objectivesusie_get_pipsusie_get_posterior_meansusie_get_posterior_samplessusie_get_posterior_sdsusie_get_prior_variancesusie_get_residual_variancesusie_init_coefsusie_plotsusie_plot_changepointsusie_plot_iterationsusie_post_outcome_configurationsusie_rsssusie_rss_lambdasusie_sersusie_sssusie_trendfiltersusie_workhorseunivar.orderunivariate_regression

Dependencies:clicpp11cpp11armadillocrayonfarverggplot2gluegtableirlbaisobandlabelinglatticelifecycleMatrixmatrixStatsmixsqpplyrR6RColorBrewerRcppRcppArmadilloRcppParallelreshapeRfastrlangS7scalesvctrsviridisLitewithrzigg

News and Updates
Release notes | Version 2.0.0 - Major Release | Release Overview | New Features | Enhancements

Last update: 2026-06-12
Started: 2025-11-20

Fine-mapping with summary statistics
The data-set | Summary statistics from simple regression | Fine-mapping with susieR using summary statistics | Single-effect summary-statistics analysis without LD | Fine-mapping with susieR using LD matrix from reference panel | Session information

Last update: 2026-05-24
Started: 2018-10-16

Modeling and Accounting for LD Reference Mismatch in Summary Statistics Fine-mapping
Recap: a "preventable" allele coding artifact | Example using real GWAS data and a 1000 Genomes European LD reference panel | Ordinary SuSiE-RSS | Kriging diagnostic | Finite LD reference panel correction | Empirical Bayes mismatch correction | Diagnostic metrics | Diagnostic with one effect | Summary | Caution: allele coding artifacts are outside this model | Session information | Reference

Last update: 2026-05-13
Started: 2026-05-02

Post-hoc credible-set filtering without an LD reference
Example | References | Session information

Last update: 2026-05-13
Started: 2026-05-13

Compare susie_rss variants

Last update: 2026-05-03
Started: 2022-04-05

Diagnostic for fine-mapping with summary statistics
LD information from the original genotype data | LD information from the reference panel | Session information

Last update: 2026-05-03
Started: 2021-05-31

Fine-mapping example
The data-set | Simple regression summary statistics | Fine-mapping with susieR | Credible sets | Posterior inclusion probabilities | Choice of priors | A note on covariate adjustment | Sufficient statistics: compute_suff_stat and susie_ss | Session information

Last update: 2026-04-25
Started: 2018-06-27

Accounting for uncertainty in residual variances for small sample studies
Data | Baseline SuSiE fit | SuSiE with Servin-Stephens SER | References

Last update: 2026-04-21
Started: 2025-09-24

Fine-mapping with SuSiE-ash and SuSiE-inf
Data | Summary Statistics and Z-Scores | Step 1: Standard SuSiE and False Positives | Step 2: Increasing Purity to Reduce False Positives | Step 3: Fitting SuSiE-inf | Step 4: SuSiE-ash Achieves the Middle Ground | Summary | What if we increase L for standard SuSiE? | Session Information

Last update: 2026-03-01
Started: 2025-12-04

Refine SuSiE model
Session information

Last update: 2025-11-20
Started: 2021-04-11

SuSiE with L0Learn initialization example
Simulate data | Fit L0Learn | Build an initialization object | Run susieR with initialization | References

Last update: 2025-11-20
Started: 2018-07-19

Evaluation of sparse version of SuSiE
Set up environment | Overview | Simulate data | X in a dense form | X in a sparse form | Further step

Last update: 2025-11-20
Started: 2018-09-18

Trend filtering
Introduction | Examples

Last update: 2025-11-20
Started: 2018-06-27

A minimal example
Session information

Last update: 2021-03-23
Started: 2018-06-27

Readme and manuals

Help Manual

Help pageTopics
Ordering of Predictors from Coefficient Estimatesabsolute.order
Block coordinate ascent for iterative model refinement.block_coordinate_ascent
Extract Regression Coefficients from Mr.ASH Fitcoef.mr.ash
Extract regression coefficients from susie fitcoef.susie
Per-Position Marginal OLS Regression of `Y` on Each Column of `X`compute_marginal_bhat_shat
Compute sufficient statistics for input to 'susie_ss'compute_suff_stat
Simulated Small-sample eQTL Data.data_small
Estimate s in 'susie_rss' Model Using Regularized LDestimate_s_rss
Simulated Fine-mapping Data with Convergence Problem.FinemappingConvergence
Get Correlations Between CSs, using Variable with Maximum PIP From Each CSget_cs_correlation
Approximation Posterior Expectations from Mr.ASH Fitget.full.posterior
Compute Distribution of z-scores of Variant j Given Other z-scores, and Detect Possible Allele Switch Issuekriging_rss
Multiple Regression with Adaptive Shrinkagemr.ash
Bayesian Multiple Regression with Mixture-of-Normals Prior (RSS)mr.ash.rss
Simulated Fine-mapping Data with Two Effect VariablesN2finemapping
Simulated Fine-mapping Data with Three Effect Variables.N3finemapping
Ordering of Predictors by Regularization Pathpath.order
Predict Outcomes or Extract Coefficients from Mr.ASH Fitpredict.mr.ash
Predict outcomes or extract coefficients from susie fit.predict.susie
Print a summary.susie_post_outcome_configuration objectprint.summary.susie_post_outcome_configuration
Real-data SuSiE-RSS example with R-reference mismatch.rss_mismatch_example
Slot Activity Prior for SuSiEslot_prior_betabinom slot_prior_poisson
Summarize Susie Fit.print.summary.susie summary.susie
Summarise a susie_post_outcome_configuration resultsummary.susie_post_outcome_configuration
Simulated Fine-mapping Data with LD matrix From Reference Panel.SummaryConsistency
Sum of Single Effects (SuSiE) Regressionsusie
Attempt at Automating SuSiE for Hard Problemssusie_auto
Inferences From Fitted SuSiE Modelsusie_get_cs susie_get_cs_attainable susie_get_lfsr susie_get_niter susie_get_objective susie_get_pip susie_get_posterior_mean susie_get_posterior_samples susie_get_posterior_sd susie_get_prior_variance susie_get_residual_variance
Initialize a susie object using regression coefficientssusie_init_coef
SuSiE Plots.susie_plot susie_plot_iteration
Plot changepoint data and susie fit using ggplot2susie_plot_changepoint
Post-hoc causal-configuration probabilities for one or more SuSiE-class fitssusie_post_outcome_configuration
SuSiE with Regression Summary Statistics (RSS)susie_rss
Sum of Single Effects Regression using the RSS-lambda likelihoodsusie_rss_lambda
Single-effect regression from summary statisticssusie_ser
SuSiE using Sufficient Statisticssusie_ss
Apply susie to trend filtering (especially changepoint problems), a type of non-parametric regression.susie_trendfilter
Ordering of Predictors from Univariate Regressionunivar.order
Perform Univariate Linear Regression Separately for Columns of Xcalc_z univariate_regression
Simulated Fine-mapping Data with Sparse, Oligogenic and Polygenic Effects.unmappable_data