[R] How do I show real values on a log10 histogram

John Kane jrkrideau at inbox.com
Wed Mar 27 15:02:56 CET 2013


Well, I don't think that is what Shane wants either but my suggestion is clearly wrong.  I was reading the question as a dual axis not just a suplimentary , equivalent value axis.  DUH.

What I think he wants is what IIRC  Petr orginally suggested which is this..
 Create the data to be graphed
        x<-1:5
        y1<-x
        y2<-x^2

# Set the par values
       # op  <- par(las=1,xaxs="r",mai=c(1,1,1,1))

# Draw the  plot        
        plot(x,y1,xlim=c(0,5),ylim=c(0,5), ylab="y1", las = 1)
        title(main="Equvalent Axes in R")
        axis(3, labels = c("10", "100", "1,000", "10,000", "100,000"), at=1:5)
        text(12, 50, "y2", srt = 270, xpd = TRUE)
        par(op)  # reset par 



John Kane
Kingston ON Canada


> -----Original Message-----
> From: maechler at stat.math.ethz.ch
> Sent: Wed, 27 Mar 2013 14:37:58 +0100
> To: careyshan at gmail.com
> Subject: Re: [R] How do I show real values on a log10 histogram
> 
>>>>>> Shane Carey <careyshan at gmail.com>
>>>>>>     on Tue, 26 Mar 2013 11:03:20 +0000 writes:
> 
>     > Yup, Ive tried all these things and I think Johns might be
>     > the best approach,
> 
> well, :-)
> you have not yet seen the following one :
> 
> if(!require("sfsmisc")) install.packages("sfsmisc")
> require("sfsmisc")
> 
> ## the data:
>        set.seed(1); summary(x <- rlnorm(100, m = 2, sdl = 3))
> ## the plot (w/o  x-axis) :
>        r <- hist(log10(x), xaxt = "n", xlab = "x [log scale]")
> ## the nice axis:
>        axt <- axTicks(1)
>        eaxis(1, at = axt, labels = pretty10exp(10^axt, drop.1=TRUE))
> 
> 
> Martin Maechler,
> ETH Zurich
> 
> 
>     > thanks
> 
>     > On Tue, Mar 26, 2013 at 11:01 AM, PIKAL Petr
>     > <petr.pikal at precheza.cz> wrote:
> 
>     >> Hi
>     >>
>     >> maybe axis(3, (x)^10)
>     >>
>     >> Regards Petr
>     >>
>     >> > -----Original Message----- > From:
>     >> r-help-bounces at r-project.org [mailto:r-help-bounces at r- >
>     >> project.org] On Behalf Of Shane Carey > Sent: Tuesday,
>     >> March 26, 2013 11:19 AM > To: r-help at r-project.org >
>     >> Subject: [R] How do I show real values on a log10
>     >> histogram
>     >> >
>     >> > Hi,
>     >> >
>     >> > I have a histogram with values logged to the base 10
>     >> >
>     >> > hist(log10(x),breaks=60)
>     >> >
>     >> > How do I show the log values on the x-axis and a second
>     >> x-axis showing > the real values?
>     >> >
>     >> > Thanks
>     >> >
>     >> > --
>     >> > Shane
>     >> >
>     >> > [[alternative HTML version deleted]]
>     >> >
>     >> > ______________________________________________ >
>     >> R-help at r-project.org mailing list >
>     >> https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do
>     >> read the posting guide http://www.R-project.org/posting-
>     >> > guide.html > and provide commented, minimal,
>     >> self-contained, reproducible code.
>     >>
> 
> 
> 
>     > --
>     > Shane
> 
>     > 	[[alternative HTML version deleted]]
> 
>     > ______________________________________________
>     > R-help at r-project.org mailing list
>     > https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do
>     > read the posting guide
>     > http://www.R-project.org/posting-guide.html and provide
>     > commented, minimal, self-contained, reproducible code.
> 
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

____________________________________________________________
Receive Notifications of Incoming Messages
Easily monitor multiple email accounts & access them with a click.
Visit http://www.inbox.com/notifier and check it out!



More information about the R-help mailing list