spatstat (version 1.64-1)

plot.msr: Plot a Signed or Vector-Valued Measure


Plot a signed measure or vector-valued measure.


# S3 method for msr
plot(x, …,
                     add = FALSE,
                     how = c("image", "contour", "imagecontour"),
                     main = NULL,
                     do.plot = TRUE,
                     multiplot = TRUE,
                     massthresh = 0,
                     equal.markscale = FALSE,
                     equal.ribbon = FALSE)



The signed or vector measure to be plotted. An object of class "msr" (see msr).

Extra arguments passed to Smooth.ppp to control the interpolation of the continuous density component of x, or passed to or plot.ppp to control the appearance of the plot.


Logical flag; if TRUE, the graphics are added to the existing plot. If FALSE (the default) a new plot is initialised.


String indicating how to display the continuous density component.


String. Main title for the plot.


Logical value determining whether to actually perform the plotting.


Logical value indicating whether it is permissible to display a plot with multiple panels (representing different components of a vector-valued measure, or different types of points in a multitype measure.)


Threshold for plotting atoms. A single numeric value or NULL. If massthresh=0 (the default) then only atoms with nonzero mass will be plotted. If massthresh > 0 then only atoms whose absolute mass exceeds massthresh will be plotted. If massthresh=NULL, then all atoms of the measure will be plotted.


Logical value indicating whether different panels should use the same symbol map (to represent the masses of atoms of the measure).


Logical value indicating whether different panels should use the same colour map (to represent the density values in the diffuse component of the measure).


(Invisible) colour map (object of class "colourmap") for the colour image.


This is the plot method for the class "msr".

The continuous density component of x is interpolated from the existing data by Smooth.ppp, and then displayed as a colour image by

The discrete atomic component of x is then superimposed on this image by plotting the atoms as circles (for positive mass) or squares (for negative mass) by plot.ppp. By default, atoms with zero mass are not plotted at all.

To smooth both the discrete and continuous components, use Smooth.msr.

Use the argument clipwin to restrict the plot to a subset of the full data.

To remove atoms with tiny masses, use the argument massthresh.

See Also

msr, Smooth.ppp, Smooth.msr,, plot.ppp


   X <- rpoispp(function(x,y) { exp(3+3*x) })
   fit <- ppm(X, ~x+y)
   rp <- residuals(fit, type="pearson")
   rs <- residuals(fit, type="score")

   plot(rs, how="contour")
# }