“…In addition, numerical methods are needed for computing the function in (1), on which the probabilities and the likelihood function are based. For the Thomas process, the inner integral in , that is, , may be computed using an efficient numerical method described in DiDonato and Jarnagin (), which is implemented in, for example, the pmvnEll function in the package shotGroups (Wollschlaeger, ) written for use in r (R Core Team, ). If the point process is a Matérn cluster process, may be computed analytically (Appendix S2).…”