In the R-package spatstat there is a function called ppm. This package augments ppm with a Baysian fitting algorithm.