Abstract
We generalize the log Gaussian Cox process (LGCP) framework to model multiplecorrelated point data jointly. The resulting log Gaussian Cox process network(LGCPN) considers the observations as realizations of multiple LGCPs, whose logintensities are given by linear combinations of latent functions drawn fromGaussian process priors. The coefficients of these linear combinations are alsodrawn from Gaussian processes and can incorporate additional dependencies apriori. We derive closed-form expressions for the moments of the intensityfunctions in our model and use them to develop an efficient variationalinference algorithm that is orders of magnitude faster than competingdeterministic and stochastic approximations of multivariate LGCP andcoregionalization models. Our approach outperforms the state of the art injointly estimating multiple bovine tuberculosis incidents in Cornwall, UK, andmultiple crime type intensities across New York city.