We present a new Monte Carlo technique to investigate the signal delay of advanced BJTs featuring relevant non-local effects. The method is suited to analyze base and collector signal delays in presence of significant nonequilibrium transport effects and up to high currents, and to verify the physical meaning and applicability of delay expressions for compact models.