Re: [R] Metropolis-Hastings Markov Chain Monte Carlo in Spatstat

From: roger koenker <>
Date: Wed, 06 Jun 2007 17:09:46 -0500

Take a look at: and in particular the link there to the original paper by Brian Kernighan describing ratfor; it is only 14 pages, but it is a model of clarity of exposition and design.

I wouldn't worry too much about the makefile -- it probably knows exactly what to do with ratfor provided you have the ratfor preprocessor available from the above link, and the rest of the tools to build from source.

url:            Roger Koenker
email            Department of Economics
vox:     217-333-4558                University of Illinois
fax:       217-244-6678                Champaign, IL 61820

On Jun 6, 2007, at 4:42 PM, Kevin C Packard wrote:

> I'm testing some different formulations of pairwise interaction
> point processes
> in Spatstat (version 1.11-6) using R 2.5.0 on a Windows platform
> and I wish to
> simulate them using the Metropolis-Hastings algorithm implemented
> with Spatstat.
> Spatstat utilizes Fortran77 code with the preprocessor RatFor to do
> the
> Metropolis-Hastings MCMC, but the Makefile is more complicated than
> any I have
> worked with.
> Any suggestions on how I could get started working with the Fortran
> code in
> conjunction with RatFor is appreciated.
> Sincerely,
> Kevin
> Kevin Packard
> Department of Forestry, PhD student
> Department of Statistics, MS student
> Virginia Polytechnic Institute and State University
> Blacksburg, Virginia, USA
> ______________________________________________
> mailing list
> PLEASE do read the posting guide
> guide.html
> and provide commented, minimal, self-contained, reproducible code. mailing list PLEASE do read the posting guide and provide commented, minimal, self-contained, reproducible code. Received on Wed 06 Jun 2007 - 22:22:57 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 06 Jun 2007 - 22:31:31 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.