We numerically estimate the divergence of several two-vertex diagrams that contribute to the radiative corrections for the Lorentzian EPRL spin foam propagator. We compute the amplitudes as functions of a homogeneous cutoff over the bulk quantum numbers, fixed boundary data, and different Immirzi parameters, and find that for a class of two-vertex diagrams, those with fewer than six internal faces are convergent. The calculations are done with the numerical framework sl2cfoam-next.