Function to perform an EM algorithm (LDA projection + 1D KDE (nonparametric, correlation-aware)) for a 2-component mixture where each observation iii has a p-value pip_ipi, and an auxiliary feature vector viv_ivi.