Ограничения на параметры подобранной модели в spatstat

В Rпакете spatstatя хотел бы наложить некоторые ограничения на параметры модели, подходящей для моих данных.

Мои данные - это многотипный pppобъект, типы которого называются BRPи GluR. Как видно на графике ниже, точки (примерно) организованы в концентрические круги. BRPобъекты локализуются на кольце с радиусом 128,4, а GluRорганизованы на внешнем кольце с радиусом 261,8.

Using the ppm function, I would thus like to fit a model where the intensity function for type BRP scales with ~abs(sqrt(x^2+y^2)-128.4), and where the intensity function for type GluR scales with ~abs(sqrt(x^2+y^2)-261.8). However, when using the following trend formula:

fit_NMJ <- ppm(NMJ,~marks*(abs(sqrt(x^2+y^2)-128.4)+abs(sqrt(x^2+y^2)-261.8)))

I obtain cross-interactions between GluR and ~abs(sqrt(x^2+y^2)-radius_BRP ), and between BRP and ~abs(sqrt(x^2+y^2)-radius_GluR) (which I want to avoid):

> radius_GluR <- 261.8
> radius_BRP <- 128.4
> fit_NMJ <- ppm(NMJ,~marks*(abs(sqrt(x^2+y^2)-radius_BRP )+abs(sqrt(x^2+y^2)-radius_GluR )))
> fit_NMJ
Nonstationary multitype Poisson process

Possible marks: ‘BRP’ and ‘GluR’

Log intensity:  ~marks * (abs(sqrt(x^2 + y^2) - radius_BRP) + abs(sqrt(x^2 + y^2) - radius_GluR))

Fitted trend coefficients:
                                 (Intercept)                                    marksGluR 
                                -9.850558359                                  0.231315028 
           abs(sqrt(x^2 + y^2) - radius_BRP)           abs(sqrt(x^2 + y^2) - radius_GluR) 
                                -0.044039428                                  0.009055454 
 marksGluR:abs(sqrt(x^2 + y^2) - radius_BRP) marksGluR:abs(sqrt(x^2 + y^2) - radius_GluR) 
                                 0.039652754                                 -0.018198305 

                                                 Estimate       S.E.       CI95.lo      CI95.hi Ztest       Zval
(Intercept)                                  -9.850558359 1.92806275 -13.629491901 -6.071624817   *** -5.1090445
marksGluR                                     0.231315028 2.00684749  -3.702033781  4.164663836        0.1152629
abs(sqrt(x^2 + y^2) - radius_BRP)            -0.044039428 0.01749010  -0.078319402 -0.009759454     * -2.5179626
abs(sqrt(x^2 + y^2) - radius_GluR)            0.009055454 0.01392686  -0.018240694  0.036351602        0.6502150
marksGluR:abs(sqrt(x^2 + y^2) - radius_BRP)   0.039652754 0.01784817   0.004670991  0.074634517     *  2.2216710
marksGluR:abs(sqrt(x^2 + y^2) - radius_GluR) -0.018198305 0.01478904  -0.047184295  0.010787685       -1.2305263

i.e. I would like to constraint the parameters marksGluR:abs(sqrt(x^2 + y^2) - radius_BRP) and marksBRP:abs(sqrt(x^2 + y^2) - radius_GluR) to be zero. Is it possible to specify this in the trend formula ? Does it also work for multiple independent datasets (when using the mppm function) ?

У меня также есть возможность подобрать эти два типа по отдельности, но при подгонке разных непуассоновских моделей с точечными взаимодействиями я могу пропустить некоторые межтиповые взаимодействия (как, например, при подборе MultiStraussмодели).