I'm looking for an equivalent Python function to perform finite mixtures of regression models using the EM algorithm with an R function on PCs to analyse some latent clusters. Does anyone know if such a function exists? Inparticular an equivalent package of the following code #STEP FLEX MIX library(flexmix) scores_PC1_2or<- as.data.frame(score_PC1_2[,1:2]) SF_MIX <- stepFlexmix( PC2 ~ PC1, data = score_PC1_2[,1:2], k = 3, nrep = 50) Continue reading...