Galaxies are often modelled as disjoint composites of distinct spectral components, implying that different wavelength ranges are only weakly correlated. They are not. We construct a data-driven model to predict infrared emission from optical spectra, achieving almost lossless predictions. Traditional fitting methods are incapable of making predictions, being biased by model misspecifications.