I've been using the mppm function of the spatstat R package in order to estimate the intensity of a non-homogeneous Poisson process with multiple realizations.
However, similarly to providing weights in glm, I need to estimate the intensity of the NHPP using different weights for the different point patterns.
Is it possible to do so using the mppm function of the spatstat package? If so, how? If not, are you aware of any other R packages that have implemented this?