We discover a family $A$ of sixteen statistics on fillings of any given Young diagram and prove new combinatorial formulas for modified Macdonald polynomials, that is, $$\tilde{H}_λ(X;q,t)=\sum_{σ\in T(λ)}x^σq^{maj(σ)}t^{η(σ)}$$ for each statistic $η\in A$. Building upon this new formula, we establish four compact formulas for the modified Macdonald polynomials, namely, $$\tilde{H}_λ(X;q,t)=\sum_σd_{\varepsilon}(σ)x^σq^{maj(σ)}t^{η(σ)}$$ which is summed over all canonical or dual canonical fillings of a Young diagram and $d_{\varepsilon}(σ)$ is a product of $t$-multinomials. Finally, the compact formulas enable us to derive four explicit expressions for the monomial expansion of modified Macdonald polynomials, one of which coincides with the formula given by Garbali and Wheeler (2020).