I am wondering if I want to model a space-time non-homogeneous Poisson process that has more probability of occurrence in some areas in space, how should I do that? For example we want to model the intensity function of earthquake occurrence that is more probable to happen in some areas (around earthquake line) and therefore the intensity function is a function of space and time. Thanks, [[alternative HTML version deleted]]