The paper considers the superposition of modified Omori functions as a conditional intensity function for a point process model used in the exploratory analysis of earthquake clusters. For the examples discussed, the maximum likelihood estimates converge well starting from appropriate initial values even though the number of parameters estimated can be large (though never larger than the number of observations). Three datasets are subjected to different analyses, showing the use of the model to discover and study individual clustering features.
