  The Dirichlet Process (DP) mixture model has become a popular choice for
model-based clustering, largely because it allows the number of clusters to be
inferred. The sequential updating and greedy search (SUGS) algorithm (Wang and
Dunson, 2011) was proposed as a fast method for performing approximate Bayesian
inference in DP mixture models, by posing clustering as a Bayesian model
selection (BMS) problem and avoiding the use of computationally costly Markov
chain Monte Carlo methods. Here we consider how this approach may be extended
to permit variable selection for clustering, and also demonstrate the benefits
of Bayesian model averaging (BMA) in place of BMS. Through an array of
simulation examples and well-studied examples from cancer transcriptomics, we
show that our method performs competitively with the current state-of-the-art,
while also offering computational benefits. We apply our approach to
reverse-phase protein array (RPPA) data from The Cancer Genome Atlas (TCGA) in
order to perform a pan-cancer proteomic characterisation of 5,157 tumour
samples. We have implemented our approach, together with the original SUGS
algorithm, in an open-source R package named sugsvarsel, which accelerates
analysis by performing intensive computations in C++ and provides automated
parallel processing. The R package is freely available from:
https://github.com/ococrook/sugsvarsel
