Achieving device-to-device reproducibility is a critical bottleneck for scalable photonic integrated circuits, as subtle variations in bottom-up epitaxial growth and fabrication severely limit yield. We present a machine learning workflow for III-V multi-quantum well microring lasers that first optimizes growth and geometry parameters via multi-objective Bayesian optimization, then leverages variational autoencoders (VAEs) to attribute residual device-to-device variability to its underlying sources. By explicitly targeting threshold variance alongside absolute performance, we demonstrate 100% lasing yield across all designs. The optimized multi-quantum well microring laser fields achieved a median lasing threshold of $16~μ\mathrm{J}\,\mathrm{cm}^{-2}\,\mathrm{pulse}^{-1}$, a $73\%$ reduction in threshold variance relative to the previously reported best values, and a median emission wavelength of $1333~\mathrm{nm}$, in the telecommunications O-band. Furthermore, to diagnose residual performance dispersion under nominally identical conditions, VAEs were used to isolate the key components of device morphology that impact performance. This analysis successfully decoupled geometric from material disorder, quantitatively linking previously unmeasured morphological variations to population-level threshold fluctuations. This data-driven workflow bridges the gap between fundamental epitaxy and reliable manufacturing, establishing a generalizable blueprint for designing and yield-optimizing complex, non-linear optoelectronic devices.

24 pages, 4 figures