title: Fitting mixture of beta distributions with flexmix in R
body: <p>Given the data (one-dimensional, non-repeated observations between 0 and 1), fit the model of one uniform (which is just fixed beta) and k beta distributions with varying k. </p>
  
  <p>After trying many approaches (my own EM implementation as well), I turned to <code>flexmix</code> package on CRAN. However, it seems there is no inbuilt family for Beta distributions (like for Poisson, Binomial and Gamma). However there is an option to write an M-step for an arbitrary mixture.</p>
  
  <p>Given the example M-driver in the vignette, <code>mymclust</code> I just got horribly entangled with the objects and dependencies and do not even know where to start.</p>
  
  <p>How do I write this driver for my task?  </p>