I am having some trouble figuring out how to use profilepl() in the R spatstat package for more complicated hybrid Gibbs models. For example, I am trying to build a hybrid model with a hardcore and multiple Geyer interactions.