# Fit a Tweedie density using tweedie family function from statmod
pretend <- data.frame( y = stats::rgamma(20, shape = 1, rate = 1) )
fit <- glm(y ~ 1, data = pretend,
family = statmod::tweedie(link.power = 0, var.power = 2.1))
# Compute the AIC
logLiktweedie(fit)
Run the code above in your browser using DataLab