通过估计参数重新创建分布

Question

让我们说我有一个我想要近似的发行版：

dist1<-sample(c(rbeta((tasks-tasks/2),10,1,1), rbeta((tasks/2),1,10,1)),100)
dist2<-rbeta(n,10,10,1)
combined<-dist1+dist2[1]

我们假设我不知道'combined'是如何构建的。 如何估算'combined'的参数以获取β分布的形状值？

换句话说，我想使用模仿rbeta()的{​​{1}}生成测试版发布。