[R] histogram

Agustin Lobo alobo at paleo.ija.csic.es
Wed Jun 12 12:41:23 CEST 2002


If your initial values are 1:10,
a <- paste(paste("log(",1:10,sep=""),")",sep="")

would give you a vector string a:
> a
[1] "log(1)"  "log(2)"  "log(3)"  "log(4)"  "log(5)"  "log(6)"  "log(7)"
 [8] "log(8)"  "log(9)"  "log(10)"

You can use "a" for in the histogram in a number of ways.
Faster:

> b <- rnorm(100)
> b <- b[b>0]
> hist(b)
note tick labels in the x-axis. Here from 0 to 3 by 0.5:
> hist(b,labels=paste(paste("log(",seq(0,3,0.5),sep=""),")",sep=""))

Better:
> hist(b,axes=F)
> axis(1,labels=paste(paste("log(",seq(0,3,0.5),sep=""),")",sep=""))
> axis(2)

Agus

Dr. Agustin Lobo
Instituto de Ciencias de la Tierra (CSIC)
Lluis Sole Sabaris s/n
08028 Barcelona SPAIN
tel 34 93409 5410
fax 34 93411 0012
alobo at ija.csic.es


On Wed, 12 Jun 2002, Luca De Benedictis wrote:

> Dear friends,
> I have been following the help-list in the next two months and I would
> like to tell how much I admire the effort of such a great teem of R
> expert. I am learning a lot. Thank you.
> 
> My problem is the following. I work with a fairly large data set (700
> productive sectors, 70 countries, 25 years) and I am applying kernel
> density (using MASS or sm) to sectoral exports which have markedly
> skewed distributions, so I do the analysis on a log scale.
> What I would like when I plot histograms or densities is to have on the
> horizontal axis not the value y which is equal to log(x) but the
> expression log(x) itself (i.e. log(0.3), log(1), log(34),...)
> 
> I have played a little bit with names.arg :
> plot(density(log(x[,3]),width=0.4, names.arg="log(x)")
> 
> but I was not very successful. Any suggestion would be appreciated.
> Thanks.
> 
> Luca
> 

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list