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 do i substitute the return local of program in the if expression of the egen command?
From
"Nick Cox" <[email protected]>
To
<[email protected]>
Subject
RE: st: How do i substitute the return local of program in the if expression of the egen command?
Date
Tue, 30 Mar 2010 13:36:35 +0100
Several questions are bundled up together here at different levels. Let me try to separate them. Let's put aside questions about ` ' (right) and ' ' (wrong) on one side as understood.
1. The underlying problems
==========================
Antonis wants the row total of variables if and only if all variables are zero or positive. Otherwise put: if the row minimum is zero or positive. Given the second formulation, no program is needed:
egen rowmin = rowmin(<varlist>)
egen rowtotal = rowtotal(<varlist>) if rowmin >= 0
Antonis' other code looks like more -egen- calls once he's done this. Correct me if I'm wrong.
2. Why a call to r() didn't work initially
==========================================
Maarten explained this. Antonis wants the incarnation of his r() result as a local macro and must use ` ' to get its contents.
3. Why the second call to r() still does not work
=================================================
This is more subtle. The first call to -egen- internally calls -count-. A side-effect of that is that existing r() results are zapped. Thus by the time of his second -egen- call the r() result Antonis wants no longer exists. He is then calling -egen- with an empty -if- condition, which is illegal.
4. What could be done instead
=============================
If you wanted to program this, and didn't like the solution in #1 for some reason, it is better to join the two programs into one based on something like
gen rowtotal = 0
foreach v of local varlist {
replace rowtotal = rowtotal + cond(`v' >= 0, `v', .)
}
Adding one missing whenever a variable is negative should be sufficient to get the result desired.
Antonis had a helper program that given (say)
local varlist var1 var2 var3 var4 var5
produces
var1 >= 0 & var2 >= 0 & var3 >= 0 & var4 >= 0 & var5 >= 0
There are tricks to do that more directly.
local cond : subinstr local varlist " " ">= 0 & ", all
local cond `cond' >= 0
That is, replace each internal space by ">= 0 &" and then append a final ">= 0".
Nick
[email protected]
Martin Weiss
set trace on
then, and see where the thing goes awry, i.e. post the lines surrounding the
error...
A Loumiotis
Sorry.... I mistyped the wrong left single quotation mark in my email
reply, but in my code I used the correct left quotation mark "`" and I
still get for the !(`r(andpve)') if expression
invalid syntax
r(198);
Martin Weiss
> Seems the left single quote is not correct. Replace "'" with "`", w/o the
> quotation marks, of course...
A Loumiotis
> Thanks for your help!
>
> using `r(andpve)' resolves the first problem...
>
> but for !('r(andpve)') i get the following
>
> invalid syntax
> r(198);
Maarten buis
A Loumiotis wrote:
>>> I have written a program (with the name andpve) that
>>> generates an expression of the form:
>>>
>>> var1>=0 & var2>=0 & ... & varN>=0
>>>
>>> and returns this value as a local with name r(andpve).
>>>
>>> I then want to use r(andpve) in another progam (rowctotal)
>>> where U substitute it in the if expression of the egen command.
>>>
>>> egen `left'_a = rowtotal(`right') if r(andpve), missing
>>> egen `left'_b = rowmin(`right') if !(r(andpve))
>>
>> egen `left'_a = rowtotal(`right') if `r(andpve)', missing
>> egen `left'_b = rowmin(`right') if !(`r(andpve)')
>>> But the egen command does not seem understand the r(andpve)
>>> as an if
>>> expression or the !(r(andpve)) as an if expression.
>>>
>>> What did I do wrong? I'm new to Stata. The code is
>>> below.
>
>>> program rowctotal, rclass
>>> version 11
>>> gettoken left right: 0
>>> capture drop `left'
>>> andpve `right'
>>> di r(andpve)
>>> egen `left'_a = rowtotal(`right') if r(andpve),
>>> missing
>>> egen `left'_b = rowmin(`right') if !(r(andpve))
>>> egen `left' = rowtotal(`left'_*), missing
>>> *drop `left'_*
>>> end
>>>
>>> program andpve, rclass
>>> local k = 1
>>> while "``k''" != "" {
>>> local ++k
>>> }
>>> local --k
>>> local i = 1
>>> local andpve ""
>>> while `i' <= `k'-1 {
>>> local andpve "`andpve'``i''>=0 & "
>>> local ++i
>>> }
>>> local andpve "`andpve'``i''>=0"
>>> return local andpve `andpve'
>>> end
>>>
*
* 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/