Statalist


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

RE: st: twoway scatter ..., yscale(log) with axis on log and arithmetic scales


From   Joana Cunha-Cruz <[email protected]>
To   [email protected]
Subject   RE: st: twoway scatter ..., yscale(log) with axis on log and arithmetic scales
Date   Mon, 20 Oct 2008 11:07:57 -0700 (PDT)

Dear Nick,
Thank you! 
It worked well. I'll use log10.
Thanks again.
Joana



--- On Sun, 10/19/08, Nick Cox <[email protected]> wrote:

> From: Nick Cox <[email protected]>
> Subject: RE: st: twoway scatter ..., yscale(log) with axis on log and arithmetic scales
> To: [email protected]
> Date: Sunday, October 19, 2008, 8:30 AM
> If the readership is unfamiliar with logarithms, they will
> probably feel
> more comfortable with base 10 logarithms. Clearly, using
> -log10()- not
> -log()- would answer that. 
> 
> Nick 
> [email protected] 
> 
> Maarten buis
> 
> --- Joana Cunha-Cruz <[email protected]> wrote:
> > What I want with the graph is to show the values of
> the data on both
> > scales (for the reader not familiar with log scale,
> they can see how
> > the numbers correspond from one scale to the other).
> 
> Notice that by default the numbers on the axes are still on
> the
> original (not log) scale, so I think adding extra axes with
> the log
> values just adds confusion to those unfamiliar with them,
> rather than
> solves it. If you still want to do that than this is how
> you do it:
> 
> *--------------------- begin example
> --------------------------
> sysuse auto, clear
> twoway scatter price mpg, yscale(log) xscale(log)
> forvalues n= 5000(5000)15000{
> 	local ylab `ylab' `n'
> "`=string(log(`n'),"%9.2f")'"
> }
> forvalues n = 10(10)40 {
> 	local xlab `xlab' `n'
> "`=string(log(`n'),"%9.2f")'"
> }
> 
> twoway scatter price mpg, yaxis(1) xaxis(1)        ///
>        yscale(log axis(1)) xscale(log axis(1)) ||  ///
>        scatter price mpg, yaxis(2) xaxis(2)        ///
>        yscale(log axis(2)) xscale(log axis(2))     ///
>        ylab(`ylab', axis(2)) xlab(`xlab', axis(2))
> ///
>        msymbol(i) legend(off)                      ///
>        ytitle("ln(price)", axis(2))              
>  ///
>        xtitle("ln(mpg)", axis(2))
> *--------------------- end example
> ----------------------------
> (For more on how to use examples I sent to the Statalist,
> see
> http://home.fsw.vu.nl/m.buis/stata/exampleFAQ.html )
> 
> *
> *   For searches and help try:
> *   http://www.stata.com/help.cgi?search
> *   http://www.stata.com/support/statalist/faq
> *   http://www.ats.ucla.edu/stat/stata/

__________________________________________________
Do You Yahoo!?
Tired of spam?  Yahoo! Mail has the best spam protection around 
http://mail.yahoo.com 
*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/



© Copyright 1996–2024 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index