r - 对广义 beta 素数分布的两个 pdf 进行卷积

Question

我有两个来自广义 beta 素数分布的 pdf，并希望通过卷积将它们组合成一个 pdf。在R中这样的事情可能吗？例如，我如何结合下面给定的后验分布 Y1 和 Y2？

require(extraDistr)
X <- seq(2,10,length.out = 100)
Y1 <- dnsbeta(x=X, shape1=10, shape2=5, min = min(X), max = max(X))
Y2 <- dnsbeta(x=X, shape1=7, shape2=4, min = min(X), max = max(X))

plot(X,Y1,type="l",ylim=c(range(Y1,Y2)))
lines(X,Y2,col="red")

如果不是卷积，蒙特卡洛方法会更好吗？