SSS_postmean_Laplace

Compute posterior means of data points for the Laplace(lambda) slab

Contains fast functions to calculate the exact Bayes posterior
for the Sparse Normal Sequence Model, implementing the algorithms
described in Van Erven and Szabo (2021,
<doi:10.1214/20-BA1227>). For general hierarchical
priors, sample sizes up to 10,000 are feasible within half an hour
on a standard laptop. For beta-binomial spike-and-slab priors, a
faster algorithm is provided, which can handle sample sizes of
100,000 in half an hour. In the implementation, special care has
been taken to assure numerical stability of the methods even for
such large sample sizes.

Tim van Erven

SequenceSpikeSlab

Exact Bayesian Model Selection Methods for the Sparse Normal
Sequence Model

Steven de Rooij

Botond Szabo

SSS_postmean_Laplace function

<dl><dt>x</dt>
<dd>Data vector of length n</dd>
<dt>logpsi</dt>
<dd>Vector of length n that represents a preprocessed version of the data.
It should contain the logs of the psi densities of the data points, as produced 
by <code>SSS_log_phi_psi_Laplace</code>.</dd>
<dt>postprobs</dt>
<dd>Vector of marginal posterior slab probabilities that \(x[i]\) has
non-zero mean for \(i=1,...,n\).</dd>
<dt>sigma</dt>
<dd>standard deviation of observations</dd>
<dt>lambda</dt>
<dd>parameter of Laplace slab density</dd></dl>

Arguments

Compute posterior means of data points for the Laplace(lambda) slab — SSS_postmean_Laplace

<dl>

<dt>x</dt>
<dd>Data vector of length n</dd>


<dt>logpsi</dt>
<dd>Vector of length n that represents a preprocessed version of the data.
It should contain the logs of the psi densities of the data points, as produced 
by <code>SSS_log_phi_psi_Laplace</code>.</dd>


<dt>postprobs</dt>
<dd>Vector of marginal posterior slab probabilities that \(x[i]\) has
non-zero mean for \(i=1,...,n\).</dd>


<dt>sigma</dt>
<dd>standard deviation of observations</dd>


<dt>lambda</dt>
<dd>parameter of Laplace slab density</dd>

</dl>

SSS_postmean_Laplace: Compute posterior means of data points for the Laplace(lambda) slab

Description

Usage

Value

Arguments