Fast Bayesian Estimation of Point Process Intensity as Function of Covariates