Composite probability distribution of M₁₅₀₀ and SFR for the REF (top), DEC+NEB (centre) and RIS+NEB model (bottom) for the sample at z ~ 4 as determined for each galaxy from our 1000 Monte Carlo simulations. The overlaid points show the median value properties for each object in the sample with black dots for “weak” nebular emitters and white dots for “strong” nebular emitters. The overlaid contours indicate the 68% integrated probabilities on the ensemble properties measured from the centroid of the distribution. The dashed line represents the Kennicutt relation (Kennicutt 1998).