# Load and aggregate the 'seizure' database
data(seizure)
aggCounts <- aggregate(x = cbind(seizure$count, seizure$trx), by =
list(seizure$id), FUN = sum)
# Estimate the parameters
mleSeizure <- mlePoissonTweedie(x = aggCounts[,2], a.ini = 0, D.ini =
10)
# Print
mleSeizure
# Extract loglikelihood
logLik(mleSeizure)
# Compute confidence inerval
confint(mleSeizure)
Run the code above in your browser using DataLab