# fitting con fitdistrplus
install.packages("fitdistrplus")
library(fitdistrplus)

set.seed(1)
n <- 100
x <- rgamma(n = n, shape = 2, rate = 0.5)


# max verosimiglianza

# metodo dei momenti

# metodo dei quantili

# confronto weibull, lognormale, gamma

# analisi e confronto