From: Gabor Grothendieck <ggrothendieck_at_gmail.com>

Date: Fri 17 Jun 2005 - 12:55:45 EST

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 Jun 17 12:59:59 2005

Date: Fri 17 Jun 2005 - 12:55:45 EST

On 6/16/05, Adaikalavan Ramasamy <ramasamy@cancer.org.uk> wrote:

> I am assuming that you want to do this empirically :

*>
**> x0 <- 0.899
**> x <- c( rnorm(6000), rnorm(4000, mean=3) )
**> plot( d <- density(x) )
**>
**> y0 <- approx( d$x, d$y, xout=x0 )$y # height at x0
**>
**> segments( x0, 0, x0, y0, col=2 )
**>
*

The following displays 4 arrows. Choose the one you like best:

# next 5 lines from previous post:

x0 <- 0.899

x <- c( rnorm(6000), rnorm(4000, mean=3) )
plot( d <- density(x) )

y0 <- approx( d$x, d$y, xout=x0 )$y # height at x0
# segments( x0, 0, x0, y0, col=2 )

points(x0,y0, pch=20, col="red") # red dot text(x0,y0, "\\da",vfont=vf, pos=3, col="green") text(x0,y0, "\\ua",vfont=vf, pos=1, col="red") yax <- par("usr")[3] # xaxis was drawn at this y level text(x0,yax, "\\ua", vfont=vf, pos=1, col="blue", xpd=TRUE, offset=0) text(x0,yax, "\\da", vfont=vf, pos=3, col="red", xpd=TRUE, offset=0)

The above uses xpd= which allows drawing outside of the main plot area and usr= which we need to get the position of the y coordinate of the boundary. It also uses the Hershey font for the arrow. Other arguments are described on the appropriate page. See:

?points ?text ?par # info on xpd= and usr= ?Hershey demo(Hershey) # examples of arrows ______________________________________________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 Jun 17 12:59:59 2005

*
This archive was generated by hypermail 2.1.8
: Fri 03 Mar 2006 - 03:32:44 EST
*