This study aims to estimate the latent interaction effect in the CLPM model through a two-step multiple imputation analysis. The estimation of within × within and between × within latent interaction under the CLPM model framework is compared between the one-step Bayesian LMS method and the two-step multiple imputation analysis through a simulation study. The analysis show that the two-step multiple imputation analysis can provide unbiased estimation parameter, similar to the one-step Bayesian LMS method. This study also uses self-esteem and depression data from NLSY79 to perform a two-step multiple imputation analysis of the CLPM model, as well as an empirical example of latent interaction analysis. Mplus syntax is provided for researchers' reference.