Use common axial limits for scatter matrix matrix panels (lattice distribution function)

Question

If several variables are on the same scale, it is useful to use the same axes for each panel of the scatter chart matrix to facilitate comparison. By default, it splomuses a more spatially efficient method of using scales defined by the range of each variable separately, and I cannot figure out how to override this. (I can’t get the arguments xlimand ylimto have any effect whatsoever, so I clearly don’t understand how the various arguments interact splom.)

For example, all measurements irisare in centimeters, so it would be reasonable to build them against each other with the same limits x and y, if we wanted to. How to extend this example on the help page splom? Simplified version:splom(~iris[1:4], groups = Species, data = iris)