We investigate the effect of random columnar disorder on the superconducting phase transition of a type-II superconductor in zero applied magnetic field using numerical simulations of three dimensional XY and vortex loop models. We consider both an unscreened model, in which the bare magnetic penetration length is approximated as infinite, and a strongly screened model, in which the magnetic penetration length is of order the vortex core radius. We consider both equilibrium and dynamic critical exponents. We show that, as in the disorder free case, the equilibrium transitions of the unscreened and strongly screened models lie in the same universality class, however scaling is now anisotropic. We find for the correlation length exponent ν = 1.2 ± 0.1, and for the anisotropy exponent ζ = 1.3 ± 0.1. We find different dynamic critical exponents for the unscreened and strongly screened models.