Kevin C Packard
2007-Jun-06 21:42 UTC
[R] Metropolis-Hastings Markov Chain Monte Carlo in Spatstat
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
roger koenker
2007-Jun-06 22:09 UTC
[R] Metropolis-Hastings Markov Chain Monte Carlo in Spatstat
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.