RELEASE NOTES spatstat 1.18-1 19 march 2010 We thank Andrew Bevan and Ege Rubak for contributions. OVERVIEW New 'hard core' interaction. Envelopes for 3D point patterns. Extended functionality for rmh. Bug fixes. NEW FUNCTIONS Hardcore Hard core interaction (for use in ppm). envelope.pp3 Simulation envelopes for 3D point patterns npoints Number of points in a point pattern of any kind. IMPROVEMENTS/ENHANCEMENTS rmh.ppm, rmhmodel.ppm rmh.default, rmhmodel.default Hard core models can be simulated rmh.ppm, rmhmodel.ppm It is now possible to simulate Gibbs point process models that are fitted to multitype point patterns using a non-multitype interaction, e.g. data(amacrine) fit <- ppm(amacrine, ~marks, Strauss(0.1)) rmh(fit, ...) rmh.default, rmhmodel.default The argument 'par' is now required to be a list, in all cases (previously it was sometimes a list and sometimes a vector). Fest Computation of 'envelope(X, Fest, ...)' has been accelerated. summary.pp3 now returns an object of class 'summary.pp3' containing useful summary information. It is plotted by 'plot.summary.pp3'. intersect.owin, union.owin These functions now handle any number of windows. F3est, G3est, K3est These functions now accept 'correction="best"' envelope.ppp, envelope.ppm, envelope.kppm Argument lists have changed slightly BUGS [1] Fest The result of Fest(X, correction="rs") had a slightly corrupted format, so that envelope(X, Fest, correction="rs") in fact computed the envelopes based on the "km" correction. (Spotted by Ege Rubak). Fixed. [2] rmh (rmh.ppm, rmhmodel.ppm) On some systems and some versions of R, rmh failed for non-stationary point process models, with a message about "missing value where TRUE/FALSE needed". (Spotted by Andrew Bevan). Fixed. [3] diagnose.ppm, lurking Calculations were not always correct if the model had conditional intensity equal to zero at some locations. Fixed. [4] ppm, profilepl If data points are illegal under the model (i.e. if any data points have conditional intensity equal to zero) the log pseudolikelihood should be -Inf but was sometimes returned as a finite value. Consequently, profilepl() did not always work correctly for models with a hard core. Fixed. [5] F3est, G3est Debug messages were printed unnecessarily. Fixed.