Statistics for Joint modeling, variable selection and multiply robust estimation in mediation analysis with multiple mediators