Quantcast
Channel: Active questions tagged r - Stack Overflow
Viewing all articles
Browse latest Browse all 205491

what is the Beta valuel when simulate a hybrid Geyer model using spatstat, R package

$
0
0

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.


Viewing all articles
Browse latest Browse all 205491

Trending Articles



<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>