We present a systematic numerical study of the propagation of astrophysical magnetized relativistic jets, in the context of resistive relativistic magnetohydrodynamics (RRMHD) simulations. Simulations are obtained with the PLUTO code. The Taub equation of state is combined here for the first time with IMplicit-EXplicit Runge-Kutta routines for time-stepping, allowing a proper treatment of stiff terms in the evolutionary equation for the electric field. We investigated different values and models for the plasma resistivity coefficient, assessing their impact on the level of turbulence, the formation of current sheets and reconnection plasmoids and the electromagnetic energy content.