On numerical methods for bifurcation analysis of delay differential equations