  Generalized approximate message passing (GAMP) is a computationally efficient
algorithm for estimating an unknown signal $w_0\in\mathbb{R}^N$ from a random
linear measurement $y= Xw_0 + \epsilon\in\mathbb{R}^M$, where
$X\in\mathbb{R}^{M\times N}$ is a known measurement matrix and $\epsilon$ is
the noise vector. The salient feature of GAMP is that it can provide an
unbiased estimator $\hat{r}^{\rm G}\sim\mathcal{N}(w_0, \hat{s}^2I_N)$, which
can be used for various hypothesis-testing methods. In this study, we consider
the bootstrap average of an unbiased estimator of GAMP for the elastic net. By
numerically analyzing the state evolution of \emph{approximate message passing
with resampling}, which has been proposed for computing bootstrap statistics of
the elastic net estimator, we investigate when the bootstrap averaging reduces
the variance of the unbiased estimator and the effect of optimizing the size of
each bootstrap sample and hyperparameter of the elastic net regularization in
the asymptotic setting $M, N\to\infty, M/N\to\alpha\in(0,\infty)$. The results
indicate that bootstrap averaging effectively reduces the variance of the
unbiased estimator when the actual data generation process is inconsistent with
the sparsity assumption of the regularization and the sample size is small.
Furthermore, we find that when $w_0$ is less sparse, and the data size is
small, the system undergoes a phase transition. The phase transition indicates
the existence of the region where the ensemble average of unbiased estimators
of GAMP for the elastic net norm minimization problem yields the unbiased
estimator with the minimum variance.
