在R中拟合β分布与flexmix的混合物

Question

给定数据（0和1之间的一维，非重复观察），拟合一个均匀的模型（其仅是固定的β）和具有变化k的k beta分布。

尝试了很多方法（我自己的EM实现）后，我转向CRAN上的flexmix包。然而，似乎没有内置的Beta分布族（如Poisson，Binomial和Gamma）。但是，可以选择为任意混合物编写M步骤。

鉴于插图中的示例M-driver，mymclust我刚刚与对象和依赖项纠缠在一起，甚至不知道从哪里开始。

如何为我的任务编写此驱动程序？