A new method to correct for host star variability in multi-epoch observations of exoplanet transmission spectra