if (requireNamespace("fitdistrplus", quietly = TRUE)) {
cases <- rnbinom(n = 100, mu = 5, size = 0.7)
pois_fit <- fitdistrplus::fitdist(data = cases, distr = "pois")
geom_fit <- fitdistrplus::fitdist(data = cases, distr = "geom")
nbinom_fit <- fitdistrplus::fitdist(data = cases, distr = "nbinom")
ic_tbl(pois_fit, geom_fit, nbinom_fit)
}
Run the code above in your browser using DataLab