The dilepton invariant mass spectra and integrated branching ratios of the single and double Dalitz decays P → + − γ and P → + − + − (P = π 0 , η, η ; = e or µ) are predicted by means of a data-driven approach based on the use of rational approximants applied to π 0 , η and η transition form factor experimental data in the space-like region.