Relation between f_esc, Lyα and Lyα EW₀ for z ∼ 2.2 (stacks; see Sobral et al. 2017), z ∼ 2.6 (binning; Trainor et al. 2015) and comparison with z ∼ 0−0.3 samples (e.g. Cardamone 2009; Hayes et al. 2013; Henry et al. 2015; Yang et al. 2016, 2017; Verhamme et al. 2017), estimated from dust-corrected Hα luminosities (Eq. (1)). We show the 1σ and 2σ range for the fits at z ∼ 2.2−2.6 and z ∼ 0−0.3 separately, and find them to be consistent within those uncertainties, albeit with a potential steeper relation at higher redshift. We find a combined best fitting relation given by f_esc, Lyα = 0.0048 EW₀ ± 0.05. The observed relation is significantly away from what would be predicted based on observed UV slopes between β ≈ −2 and β ≈ −1 for LAEs (see DW10) and would require β ≈ +5 for a good fit using Eq. (5). Such red β slopes are not observed for LAEs. Modifying Eq. (5) to include the effect of ξ_ion and dust reveals that those physical parameters likely play an important role; see Sects. 3.2 and 3.3.