Четырехпараметрическое бета-распределение в R

Question

Существует ли встроенная функция для расчета четырехпараметрического бета-распределения в R? Т.е. распределение с двумя параметрами формы и двумя граничными параметрами, чтобы оно не было ограниченным [0,1]?

Я сделал свой собственный, но любопытно, если эта функциональность уже существует. Не нужно изобретать велосипед.

shp1 <- 20
shp2 <- 5
X <- seq(0,1,length.out = 100)
Y <- dbeta(x = X,shape1 = shp1, shape2 = shp2)
plot(X,Y)
# scaled between two boundaries
dbeta4param <- function(x,shp1,shp2,bnd1,bnd2){
  mask = (x>=bnd1 & x<=bnd2)
  xScale = (x-bnd1)/(bnd2-bnd1)
  mask * dbeta(xScale,shape1=shp1,shape2=shp2)/(bnd2-bnd1)
}
bnd1 <- 2
bnd2 <- 4
X2 <- seq(bnd1,bnd2,length.out = 100)
Y2 <- dbeta4param(x=X2, shp1 = shp1, shp2 = shp2, bnd1 = bnd1, bnd2 = bnd2)
plot(X2,Y2)

Answer 1

CRAN Task View: вероятностные распределения - всегда хороший ресурс о том, где найти различные R распределения вероятностей. В соответствии с этим, вы должны использовать пакет extraDistr (extraDistr::rnsbeta например).