I fit a hybrid geyer model including a hardcore and a geyer saturation components, now I want to simulate the fitted model by rmh() in spatstat, the code is following:
mo.X<-rmhmodel(cif=c("hardcore","geyer"),
par=list(list(beta=0.0001746418,hc=30.31542),
list(beta=1, gamma=0.4195612,r=122,sat=2)), w=Window(A1a2.unmark))
there are two beta in simulating the hybrid model: beta for hardcore(), and beta for geyer (), however, when fitting real data to get a hybrid model, only one beta value provided in the result:
unmark.hybrid<-ppm(A1a2.unmark~1,Mo.hybrid, correction="bord")
unmark.hybrid
#beta = 0.0001052206
#Hard core distance: 33.66719
#gamma: 0.3578263
is it right to set the second beta=1? thanks for any suggestion.