[R] Is ther a possibility to display the log(frequency) instead of the frequency in the hist plot?

Witold Eryk Wolski wolski at molgen.mpg.de
Thu Jan 18 09:27:30 CET 2001




Uwe Ligges wrote:

> Witold Eryk Wolski wrote:
> >
> >  Is there a possibility to display the log(frequency) instead of the
> > frequency in the hist?
>
> I think you want something like the following example:
>
>  x <- runif(100) * 10
>  my.hist <- hist(x)
>  my.hist   # have a look at the values!
>  my.hist$counts <- log(my.hist$counts, 10)  # logarithm with base 10
>  plot(my.hist)  # now plot it again with log(frequency)
>
> Uwe Ligges
>

Hi Uwe!
First thanks for your answer.
Some questions still remain to me.
1. How take an log if one of the my.hist$counts equals 0?
2. You code expample produces an error message which i doesnt understand.
> plot(my.hist)  # now plot it again with log(frequency)
Error in plot.window(xlim, ylim, log, asp, ...) :
        invalid xlim
3. I tried as workaround something similar.
plot(my.hist$mids,log(my.hist$counts),type="s") And this would work quite good for
my needs.
But problem with log(0) = -Inf still remains.

I think I have to formulate my question more general. Why, and this is an question
which I am asking myself since
I am started using  R becouse im using the hist plot quite often,  the hist plot in
the base package has no similar
options like eg. the histogram plot in XMGR?


THX
Eryk



> -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
> 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
> _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

--
Witold Eryk Wolski
Max Plank Institut fuer Molekulare Genetik
Ihnestr73 14195 Berlin
Germany
tel.:
Work 0049-30-84131426
http://www.molgen.mpg.de/~wolski


-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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