data(pgm)
interv <- intervals(composite(pgm))
print(interv)
mean(interv) # Mean interval
# Now fit log-normal distribution instead of Weibull.
intervals(composite(pgm), densfun = "lognormal")
## Not run: ------------------------------------
# # Boxplot of fire interval distribution.
# boxplot(intervals(composite(pgm))$intervals)
## ---------------------------------------------
Run the code above in your browser using DataLab