Mediation analysis is a methodology used to understand how and why an independent variable (X) transmits its effect to an outcome (Y) through a mediator (M). New causal mediation methods based on the potential outcomes framework and counterfactual framework are a seminal advancement for mediation analysis, because they focus on the causal basis of mediation analysis. There are several programs available to estimate causal mediation effects, but these programs differ substantially in data set up, estimation, output, and software platform. To compare these programs, an empirical example is presented, and a single mediator model with treatment-mediator interaction was estimated with a continuous mediator and a continuous outcome in each program. Even though the software packages employ different estimation methods, they do provide similar causal effect estimates for mediation models with a continuous mediator and outcome. A detailed explanation of program similarities, unique features, and recommendations is discussed.