[R] finding the plot limits generated by default
marc_schwartz at me.com
Fri May 14 17:45:59 CEST 2010
On May 14, 2010, at 9:59 AM, Justin Fincher wrote:
> I have two datasets that I would like to plot in a single figure. The first
> plot is generated by a function that then takes a subset of the data. (It
> is biological data so it is usually by chromosome e.g.
> function(data1,subset="chr8") ) Since not only are the chromosomes different
> sizes, but across different datasets there may be different numbers of
> points for a single chromosome, I do not know the coordinates on the x-axis
> prior to running the function.
> I want the scale of the second plot to align to that of the first, but I
> don't know how to find the default xlim for the graph generated by the
> function. Is there any way to access this information so I can then pass it
> to the plot() of the second plot? Thanks!
> - Fincher
The result of par("usr") will give you the actual limits of the x and y axes (plot region) once the first graphic is drawn to the device. See ?par for more information. Note that by default, R will extend the data ranges by 4% (see 'xaxs' in ?par).
However, you are better off getting the common ?range of the relevant values in both subsets of data before plotting and then explicitly set the same x and y axis ranges in each plot by using the 'xlim' and 'ylim' arguments to plot().
More information about the R-help