Two-stage multilevel latent class analysis with covariates in the presence of direct effects