Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
RE: st: How to set a year index?
From
Nick Cox <[email protected]>
To
Nick Cox <[email protected]>, "'[email protected]'" <[email protected]>
Subject
RE: st: How to set a year index?
Date
Fri, 20 May 2011 16:37:07 +0100
Sorry; Oliver's code was not
local mean =
but
local basevalue =
so my point included a distortion of what he wrote. I think there remains scope for discussion about style, but I should have not misquoted Oliver's code.
Nick
[email protected]
-----Original Message-----
From: Nick Cox
Sent: 20 May 2011 16:25
To: '[email protected]'
Subject: RE: st: How to set a year index?
We're talking style, and long ago wise people warned de gustibus non est disputandum.
But, sincerely...
How does the extra code make it easier to understand? The presumption is that a reader needs the clue
local mean = r(mean)
to realise that r(mean) contains the mean. Really??? Reminds me of those places that give you warnings about hot drinks being hot. Also, the price is confusion for anyone who has no idea what a local is. So, the target readership is: people who know about locals, but not about means.
The -scalar- is no more necessary than the -local-.
I agree: -meanonly- is likely to have a trivial efficiency effect here. But it's worth learning about, and telling other people about, for problems where it matters.
Nick
[email protected]
Oliver Jones
On 20.05.2011 16:58, Nick Cox wrote:
> The logic here is fine, but why do it this way?
>
> summarize price if year == 2005, meanonly
> gen price_index = price/r(mean)
>
> Note how
>
> 1. The middle macro is completely unnecessary.
I used the macro to make the code more readable/understandable.
>
> 2. Although it will not usually bite, you are probably losing asmidgen of precision by using a local.
Good point! I'm not sure, but I think the use of a scalar can solve the precision
problem of the local. So my code would be
scalar base_value = r(mean)
gen price_index = price/base_value
>
> 3. Specifying -meanonly- is efficient if all you need is the mean.
>
I'm not an expert on computing efficiency, but I wouldn't expect to
see a difference in execution time when -summarize-ing one single
observation.
But I admit that it is good style to use the meanonly option.
> Small points indeed, but in computing the devil is often in the details.
>
> Nick
> [email protected]
*
* 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/