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

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

Take a look at: http://sepwww.stanford.edu/software/ratfor.html 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:    www.econ.uiuc.edu/~roger            Roger Koenker
email    rkoenker_at_uiuc.edu            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
>
> ______________________________________________
> R-help_at_stat.math.ethz.ch 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.



R-help_at_stat.math.ethz.ch 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. 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 https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.