Joint models for longitudinal and time-to-event data are useful in situations where an association exists between a longitudinal marker and an event time. These models are typically complicated due to the presence of shared random effects and multiple submodels. As a consequence, software implementation is warranted that is not prohibitively time consuming. While methodological research in this area continues, several statistical software procedures exist to assist in the fitting of some joint models. We review the available implementation for frequentist and Bayesian models in the statistical programming languages R, SAS, and Stata. A description of each procedure is given including estimation techniques, input and data requirements, available options for customization, and some available extensions, such as competing risks models. The software implementations are compared and contrasted through extensive simulation, highlighting their strengths and weaknesses. Data from an ongoing trial on adrenal cancer patients is used to study different nuances of software fitting on a practical example.