Current version of this R package conducts mediation path analysis for multiple mediators in two stages.
