Совокупная функция распределения смеси двух бета-распределений в R

Question

У меня есть это распределение: X ~0,75* Бета (1,3)+0,25* Бета (5,1,75)

где функция плотности задается из этой функции:

dmist2b=function(y,alpha1=1,beta1=3,p=0.25,alpha2=5,beta2=1.75) {
   (1 - p) * dbeta(y, alpha1, beta1) + p * dbeta(y, alpha2, beta2)}

Я хочу вычислить CDF этого распределения (смесь двух бета-распределений) в точке R. Вот код функции, которую я написал:

pmist2b=function(y,alpha1=1,beta1=3,p=0.25,alpha2=5,beta2=1.75){
  (1 - p) * pbeta(y, alpha1, beta1) + p * pbeta(y, alpha2, beta2)
}

Это правильно? Любая помощь будет принята с благодарностью.