library(ordbetareg)
# generate 100 random variates with an average of 0.7
# all will be in the closed interval \[0,1\]
<- rordbeta(n=100, mu=0.7, phi=2)
ordbeta_var
# Will be approx mean = 0.7 with high positive skew
summary(ordbeta_var)
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.0000 0.6471 0.9866 0.7582 1.0000 1.0000