On 24/05/2009, at 3:29 PM, echo_july wrote:
> i have trouble in using spatstat package.
> i want to simulate a community under the Strauss process,which
> has a parameter gamma that controls interaction strength between
> points,and Strauss process is defined only for 0 ??gamma ?? 1 and
> is a model for inhibition between points.
> my problem is that in my data, many species's estimated gamma
> is larger than one.so if i still wanna simulate the point process
> of that species with interaction between points,what can i do then?
> is there other method i can use?
You might consider trying the Geyer process. Use cif="geyer" in
forming
the model to be used by rmh(). In fitting Geyer models to data, one way
to estimate the saturation parameter is via profile pseudolikelihood.
See ?profilepl.
cheers,
Rolf Turner
######################################################################
Attention:
This e-mail message is privileged and confidential. If you are not the
intended recipient please delete the message and notify the sender.
Any views or opinions presented are solely those of the author.
This e-mail has been scanned and cleared by MailMarshal
www.marshalsoftware.com
######################################################################