  Structural glasses feature quasilocalized excitations whose frequencies
$\omega$ follow a universal density of states ${\cal
D}(\omega)\!\sim\!\omega^4$. Yet, the underlying physics behind this
universality is not fully understood. Here we study a mean-field model of
quasilocalized excitations in glasses, viewed as groups of particles embedded
inside an elastic medium and described collectively as anharmonic oscillators.
The oscillators, whose harmonic stiffness is taken from a rather featureless
probability distribution (of upper cutoff $\kappa_0$) in the absence of
interactions, interact among themselves through random couplings (characterized
by strength $J$) and with the surrounding elastic medium (an interaction
characterized by a constant force $h$). We first show that the model gives rise
to a gapless density of states ${\cal D}(\omega)\!=\!A_{\rm g}\,\omega^4$ for a
broad range of model parameters, expressed in terms of the strength of
stabilizing anharmonicity, which plays a decisive role in the model. Then --
using scaling theory and numerical simulations -- we provide a complete
understanding of the non-universal prefactor $A_{\rm g}(h,J,\kappa_0)$, of the
oscillators' interaction-induced mean square displacement and of an emerging
characteristic frequency, all in terms of properly identified dimensionless
quantities. In particular, we show that $A_{\rm g}(h,J,\kappa_0)$ is a
nonmonotonic function of $J$ for a fixed $h$, varying predominantly
exponentially with $-(\kappa_0 h^{2/3}\!/J^2)$ in the weak interactions (small
$J$) regime -- reminiscent of recent observations in computer glasses -- and
predominantly decaying as a power-law for larger $J$, in a regime where $h$
plays no role. We discuss the physical interpretation of the model and its
possible relations to available observations in structural glasses, along with
delineating some future research directions.
