[R] How to create a 'fit' plot

(Ted Harding) Ted.Harding at nessie.mcc.ac.uk
Fri Mar 18 17:17:12 CET 2005

On 18-Mar-05 R_xprt_wannabe wrote:
> Dear List,
> As someone who is in the process of trying to migrate
> from Excel, I'd appreciate any help on this question:
> I have a data set and want to fit, say, three
> distributions to it.  I would like to create a plot
> that shows my data points against all three fitted
> curves (estimated d.f.).  Basically, I lookint to
> creat a plot that looks like the one presented in the
> attached paper (Figure 5, page 12):
> http://www.math.ethz.ch/~mcneil/ftp/astin.pdf
> Could you please show me, or point me to example code
> showing, how that can be done?  

You can do something on the lines of:

  plot(x,empDF,pch=".") ## for the empirical plot of x

  lines(x0,GPD,lty="solid") ## for the GPD curve

  lines(x0,Pareto,lty="dashed") ## for the Pareto curve

  lines(x0,LogNormal,lty="dotdash") ## for the LogNormal curve

where x0 is a vector of (fairly finely spaced) x-values,
and LogNormal, Pareto and GPD are the corresponding y-values
(at each x0) for the corresponding curves.

Or the latter could be functions which computed the y-values,
given the x-values, in which case you might use

  lines(x0,GPD(x0),lty="solid") ## for the GPD curve


Since you're apparently beginning with R, you've encountered
one of the more deeply buried (yet commonly required) aspects
of R, namely how the details of a plot are set up.

You can consult the help on 'plot' and 'lines' with


but it takes a bit of poking around to find that you need to
look up


to find the details of things like 'lty' as used above.

You may also need to make sure that you get the axes spanning
over good ranges by making your first command


where a and b are the lower and upper limits for x, and c and d
are those for y, since otherwise 'plot' will choose these limits
in its own way depending on the ranges of values in x and in empDF,
which may not be suitable for proper display of the other graphs.
You can omit either or both of 'xlim' and 'ylim'.

Whatever you do about the limits, however, they will be fixed
once and for all once the first plot has been drawn.

Hoping this helps,

E-Mail: (Ted Harding) <Ted.Harding at nessie.mcc.ac.uk>
Fax-to-email: +44 (0)870 094 0861
Date: 18-Mar-05                                       Time: 16:17:12
------------------------------ XFMail ------------------------------

More information about the R-help mailing list