  The anomalies in the measurements of observables involving $b \to s \mu\mu$
decays, namely $R_K$, $R_{K^*}$, $P_5^{\prime}$, and $B_s^\phi$, may be
addressed by adding lepton-universality-violating new physics contributions to
the effective operators ${\cal O}_9, {\cal O}_{10}, {\cal O}^\prime_9, {\cal
O}^\prime_{10}$. We analyze all the scenarios where the new physics contributes
to a pair of these operators at a time. We perform a global fit to all relevant
data in the $b \to s$ sector to estimate the corresponding new Wilson
coefficients, $C_9^{\rm NP}, C_{10}^{\rm NP}, C_9^\prime, C_{10}^\prime$. In
the light of the new data on $R_K$ and $R_{K^*}$ presented in Moriond 2019, we
find that the scenarios with new physics contributions to the ($C_9^{\rm NP}$,
$C_9^\prime$) or ($C_9^{\rm NP}$, $C_{10}^\prime$) pair remain the most favored
ones. On the other hand, though the competing scenario ($C_9^{\rm NP}$,
$C_{10}^{\rm NP}$) remains attractive, its advantage above the SM reduces
significantly due to the tension that emerges between the $R_K$ and $R_{K^*}$
measurements with the new data. The movement of the $R_K$ measurement towards
unity would also result in the re-emergence of the one-parameter scenario
$C_9^{\rm NP} = -C_9^\prime$.
