d_beta <- dist_beta(shape1 = 2, shape2 = 2, ncp = 0)
x <- d_beta$sample(100)
d_emp <- dist_empirical(x)
plot_distributions(
empirical = d_emp,
theoretical = d_beta,
estimated = d_beta,
with_params = list(
estimated = inflate_params(
fitdistrplus::fitdist(x, distr = "beta")$estimate
)
),
.x = seq(0, 2, length.out = 100)
)
Run the code above in your browser using DataLab