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
######################################################################
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
