The $p_{\rm T}$ integrated yield (d$N/\mathrm{d}y$) (upper panels) and mean transverse momentum $\left(\langle\pt\rangle\right)$ (bottom panels) for ${\rm K^{*0}}$ (left) and $\phi$ (right) mesons as a function of $y$, divided by the d$N/\mathrm{d}y$ and $\langle p_{\mathrm{T}} \rangle$ at $y = 0$ for the multiplicity class 0$-$100$\%$ in p$-$Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV. The predictions from EPOS-LHC, EPOS3 with and without UrQMD, DPMJET, HIJING, and PYTHIA8/Angantyr are shown as different curves. The statistical uncertainties are represented as bars whereas the boxes indicate total systematic uncertainties.