# Estimating two-component mixture of Weibull distributions?

Is there any existing package (preferably in Python or Matlab) to estimate the parameters of a two-component Weibull mixture model? And failing that, I am hoping to get some pointers towards rolling out a fairly simple implementation of my own (perhaps something EM-like, a la mixture of Gaussians?).

Note: I considered using mixtures of Gaussians as an easier alternative, but my data are definitely not Gaussian, and being extremal values are quite well-modeled by Weibull.

One link for helping to roll your own implementation that came up quickly on a web search is provided by Razali and Al-Wakeel who use maximum likelihood and provide a formula for the log-likelihoods of Weibull mixtures. Each choice of combinations of Weibull components and numbers of parameters is analyzed separately. At first glance, however, their approach for choosing among those competing models seems like it might be prone to over-fitting as it appears to be based on $$R^2$$ values in the data set at hand rather than a more generalizable approach like bootstrapping.