From: Haynes, Maurice (NIH/NICHD) [E] <haynesm_at_cfr.nichd.nih.gov>

Date: Fri 23 Dec 2005 - 23:32:13 EST

R-help@stat.math.ethz.ch mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

http://www.R-project.org/posting-guide.html

R-help@stat.math.ethz.ch mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Fri Dec 23 23:51:41 2005

Date: Fri 23 Dec 2005 - 23:32:13 EST

Hi,

I believe Jim Lemon's solution requires his plotrix package.

Maurice Haynes

National Institutes of Health

Child and Family Research Section

-----Original Message-----

From: Jim Lemon [mailto:bitwrit@ozemail.com.au]
Sent: Friday, December 23, 2005 10:12 PM
To: Ronnie Babigumira

Cc: R Help

Subject: Re: [R] Plot problems: xlim

Ronnie Babigumira wrote:

*> Hi,
*

> Still fresh in R, tried to figure this out, now on my second day

running with no luck (and a pile of hair on my desk) so

> I have thrown in the towel and would like to ask for some help.

*>
**> Here is what I am trying to do. I am trying to plot a distribution, I
*

have 99 points, bound in the range

*>
*

> xlim.min: -0.0173

*> xlim.max: 0.02103
**>
**> However, I have a value outside this range (0.2454959) which I would
*

like to add to the plot as a line and to do this I

> use abline(v=0.2454959)

*>
**> This is what I write
**>
**> >xlim = c(-0.02, 0.3)
**> >denz <- density(morp)
**> >plot.density(denz, xlim = xlim, ylim = c(0,70))
**> >hist(morp, freq=F, add= T)
**> >abline(v=0.2454959)
**>
**> Without any options, plot.density spreads out nicely, however,
*

naturally, the line I want to add is not plotted since it

> is well outside the range automatically determined by plot.density

hence the need to add xlim however this produces

> something I dont find aesthetically appealing. The plot is squeezed

out into a very lean "bell" shape.

*>
*

> So (finally) my question, how can i widen the spread of my plot and

yet also be able to add my xline.

*>
*

Hi Ronnie,

For only one line, it is probably easiest to stick in an axis break and label the line on the x axis. Notice that the position of the line is arbitrarily set to be far enough beyond the end of the density curve to allow room for the axis break.

testdata<-rnorm(50,sd=0.01)

denz<-density(testdata)

plot(denz,xlim=c(-0.02,0.04),axes=FALSE)
box()

axis(1,at=c(-0.02,0,0.02,0.039),labels=c(-0.02,0,0.02,0.2454959))
abline(v=0.039)

axis.break(1,breakpos=0.037)

This is probably a common enough problem for inclusion in the axis.break

example.

Jim

R-help@stat.math.ethz.ch mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

http://www.R-project.org/posting-guide.html

R-help@stat.math.ethz.ch mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Fri Dec 23 23:51:41 2005

*
This archive was generated by hypermail 2.1.8
: Fri 03 Mar 2006 - 03:41:40 EST
*