d_nbinom <- dist_negbinomial(size = 3.5, mu = 8.75)
x <- d_nbinom$sample(100)
d_emp <- dist_empirical(x)
plot_distributions(
empirical = d_emp,
theoretical = d_nbinom,
estimated = d_nbinom,
with_params = list(
estimated = inflate_params(
fitdistrplus::fitdist(x, distr = "nbinom")$estimate
)
),
.x = 0:max(x)
)
Run the code above in your browser using DataLab