Decomposing feature-level variation with covariate Gaussian process latent variable models