Stata The Stata listserver
[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

st: RE: tempvar


From   Joly.Patrick@ic.gc.ca
To   statalist@hsphsun2.harvard.edu
Subject   st: RE: tempvar
Date   Fri, 7 Jun 2002 10:50:30 -0400

In order to help you, we need to know more about your problem:

- which line is causing you grief?  (We can't simply guess which one it is.)
- if you don't know which line Stata is choking on, you need to -set trace
on-, create a log of a sample run of your program, and then paste it in one
of your message. (I mentioned this yesterday.)
- don't send us just part of you program.  (Again, we can't guess what some
of your macros contain.  And, again, I mentioned this yesterday.)

Sometimes errors may be obvious to us, but you can't just assume that.  Help
us help you if you want any help at all.

Patrick Joly
joly.patrick@ic.gc.ca
pat.joly@utoronto.ca

> -----Original Message-----
> From: Anna-Lisa Pierre [mailto:annaf2001@hotmail.com]
> Sent: June 07, 2002 10:30 AM
> To: statalist@hsphsun2.harvard.edu
> Subject: st: tempvar
> 
> 
> Hi
> 
> I am trying to use the 'tempvar' command so that I do not 
> have to drop 
> variables at the end. However, I am having problems using 
> this within loops. 
> Are 'tempvar' and 'tempname' not supposed to be used within 
> loops? The part 
> of my programme that is giving trouble is as follows:
> Any insights will be gratefully appreciated.
> 
> tempvar wj_1 phiwj_1 PHIwj_1 diffqq2;
>   while `j' <= `depmax' { ;
> 
> 	local xless1 = `x'-1 ;
> 
> 	if `j' < `depmax' { ;
> 		ge `diffqq2' = ( normprob(_b[_cut`x']-`scor_uqq') - 
> normprob(_b[_cut`x']-`scor_nqq') )
> 			  - ( normprob(_b[_cut`xless1']-`scor_uqq') - 
> normprob(_b[_cut`xless1']-`scor_nqq') )
> 			  if `1'~=. ;
> 	} ;
> 	else { ;
> 		ge `diffqq2' = normprob(`scor_uqq'-_b[_cut`xless1']) - 
> normprob(`scor_nqq'-_b[_cut`xless1']) if `1'~=. ;
> 	} ;
> 
> 
> 	su `diffqq2' `if' ;
> 	local diff`x' = r(mean) ;
> 	scalar `check' = `check' + `diff`x'' ;
> 
> 
>         rename `wj' `wj_1' ;
> 	rename `phiwj' `phiwj_1' ;
> 	rename `PHIwj' `PHIwj_1' ;
> 
> 	if `j'~=`depmax' { ;
> 		ge `wj' = ( _b[_cut`x'] - `scoreqq' ) `if' ;
> 		ge `phiwj' = (1/sqrt(2*_pi)) * exp( (`wj'^2)/-2 ) `if' ;
> 		ge `PHIwj' = normprob(`wj') `if' ;
> 	} ;
> 	else { ;
> 		ge `wj' = 0 ;
> 		ge `phiwj' = 0 ;
> 		ge `PHIwj' = 1 ;
> 	} ;
> 
> 	replace `eqq'=( `phiwj_1' - `phiwj' ) / ( `PHIwj' - 
> `PHIwj_1' ) if `1'==`j' 
> & `ifexp' ;
> 	replace `e2qq'=( (`wj_1' * `phiwj_1') - (`wj' * `phiwj') )
> 			   / ( `PHIwj' - `PHIwj_1' )        if 
> `1'==`j' & `ifexp' ;
> 	replace `e3qq'=( ((`wj_1'^2) * `phiwj_1') - ((`wj'^2) * 
> `phiwj') )
> 			   / ( `PHIwj' - `PHIwj_1' )        if 
> `1'==`j' & `ifexp' ;
> 	replace `e4qq'=( ((`wj_1'^3) * `phiwj_1') - ((`wj'^3) * 
> `phiwj') )
> 			   / ( `PHIwj' - `PHIwj_1' )        if 
> `1'==`j' & `ifexp' ;
> 
> 
> 	tempvar y`x'iqq  y`xless1'iqq;
> 	if `j'~=`depmax' { ;
> 		capture drop `y`x'iqq' ;
> 		if _rc == 0 { local dropvar="`dropvar' `y`x'iqq`" } ;
> 		ge `y`x'iqq' = 0 `if' ;
> 		replace `y`x'iqq'=( `phiwj' ) / ( `PHIwj' - 
> `PHIwj_1' ) if `1'==`j' & 
> `ifexp' ;
> 	} ;
> 
> 	if `x'~=2 { ;
> 		replace `y`xless1'iqq'=( -1 * `phiwj_1' ) / ( 
> `PHIwj' - `PHIwj_1' )
> 						if `1'==`j' & `ifexp' ;
> 	} ;
> 	su `1' if `1'>`j' & `ifexp' ;
> 	local j = r(min) ;      /* replace j with next value of 
> dep var */
> 	local x=`x'+1 ;
> 	drop `wj_1' `phiwj_1' `PHIwj_1' ;
>   } ;
> 
> 
>   #delimit cr
> 
> _________________________________________________________________
> Join the world's largest e-mail service with MSN Hotmail. 
> http://www.hotmail.com
> 
> *
> *   For searches and help try:
> *   http://www.stata.com/support/faqs/res/findit.html
> *   http://www.stata.com/support/statalist/faq
> *   http://www.ats.ucla.edu/stat/stata/
> 
*
*   For searches and help try:
*   http://www.stata.com/support/faqs/res/findit.html
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/



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