Statalist


[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

RE: st: RE: Forvalues, syntax error


From   "Gauri Khanna" <gwkhanna@hotmail.com>
To   statalist@hsphsun2.harvard.edu
Subject   RE: st: RE: Forvalues, syntax error
Date   Thu, 05 Jul 2007 14:27:12 +0000

Dear Neil, Rajesh and Stata list,

I think why there is an error-my observations (see below) are not getting picked up when I -summarize i- and hence there can no -local max1-. I have the uneviable task of am using someone else's unanotated code!

Rajesh, the -finaltemp-had been generated but I did not copy it in the code I sent.

Thank you both and I plod away.


set more off
set matsize 5000
set trace on /*Suggestion from Stata list*/
set tracedepth 2 /*Suggestion from Stata list*/

use `allsvys1', clear
egen i=group(un_region2 sex_code) if excind!=1 & sex_code!=3 & total_age_group!=1 & agegrp==0
summarize i /*PROBLEM*/
local max1=r(max)
gen error1=.
gen error2=.
save `allsvys', replace

OUTPUT for summarize i :
summarize i

Variable | Obs Mean Std. Dev. Min Max
-------------+--------------------------------------------------------
i | 0


Gauri



From: "Rajesh Tharyan" <R.Tharyan@exeter.ac.uk>
Reply-To: statalist@hsphsun2.harvard.edu
To: <statalist@hsphsun2.harvard.edu>
Subject: st: RE: Forvalues, syntax error
Date: Thu, 5 Jul 2007 15:15:25 +0100

Hi,


One guess is with the lines for append and save

Try using the full path

Something like (replace C:\path with your drive and path)

append using "C:\path\\`file'"
save "C:\path\\`file'",replace

Also, I cannot see where the `finaltemp' is being created in the first
place..

Hope this helps
rajesh


-----Original Message-----
From: owner-statalist@hsphsun2.harvard.edu
[mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Gauri Khanna
Sent: 05 July 2007 14:25
To: statalist@hsphsun2.harvard.edu
Subject: st: Forvalues, syntax error

Dear Stata List,

I have run into a syntax error r(198) after running the command -
forvalues-.

The strange thing is that the code works when I use my entire dataset but I
get this syntax error when I use a subset of the data. Since this is a
syntax problem, my guess is that it is not a substantive issue but I am
quite puzzled as to the problem with syntax. I have copied the code and the
output for you to see. I am grateful for your help.

COMMANDS
***Regression***
set more off
set matsize 5000

use `allsvys1', clear
egen i=group(un_region2 sex_code) if excind!=1 & sex_code!=3 &
total_age_group!=1 & agegrp==0
summarize i
local max1=r(max)
gen error1=.
gen error2=.
save `allsvys', replace


forvalues svy = 1/`max1' {
use `allsvys', clear
keep if i==`svy'

summarize smk_d
local maxsd=r(max)
if `maxsd'>=95 local maxsd=95
local minsd=r(min)
if `minsd'<=5 local minsd=5

summarize cig_d
local maxcd=r(max)
if `maxcd'>=95 local maxcd=95
local mincd=r(min)
if `mincd'<=5 local mincd=5

replace lgsmkd=ln((smk_d-(`minsd'-5))/((`maxsd'+5)-smk_d))
replace lgcigd=ln((cig_d-(`mincd'-5))/((`maxcd'+5)-cig_d))

replace agesmkd=midage*lgsmkd
replace agecigd=midage*lgcigd

list un_region2 in 1
list sex_code in 1

sw regress lgcigd lgsmkd midage agesmkd [aweight=sample_size] if
excind!=1
& sex_code!=3 & total_age_group!=1 & totalind!=1 & agegrp==0 &
source_code!=101692 & source_code!=101731 & source_code!=101712 &
source_code!=101014 & source_code!=101711 & source_code!=101724, pr(0.05)
lockterm1 hier
mark type if excind!=1 & sex_code!=3 & total_age_group!=1 &
totalind!=1 &
agegrp==0 & lgcigd!=. & lgsmkd!=. & source_code!=101014 &
source_code!=101711 & source_code!=101724 & source_code!=101692 &
source_code!=101731 & source_code!=101712
predict yhat1
replace r2noc=e(r2_a)
replace cig_d_est2 =
((`mincd'-5)+(`maxcd'+5)*exp(yhat1))/(1+exp(yhat1))
replace cig_d=min(cig_d_est2, cig_c, smk_d) if cig_d == . &
cig_d_est2 !=
. & agegrp==0

****************************************************************************
******************************************************
****************************************************************************
******************************************************

sw regress lgsmkd lgcigd midage agecigd [aweight=sample_size] if
type==1,
pr(0.05) lockterm1 hier
predict yhat3
replace r2nos=e(r2_a)
replace smk_d_est2 =
((`minsd'-5)+(`maxsd'+5)*exp(yhat3))/(1+exp(yhat3))

replace trace22="typeconv cigd->smkd ok" if cig_d!=. & cig_d!=0 &
smk_d
== . & smk_d_est2 != . & smk_c==. & smk_d_est2>cig_d & agegrp==0
replace smk_d=smk_d_est2 if cig_d!=. & cig_d!=0 &
smk_d
== . & smk_d_est2 != . & smk_c==. & smk_d_est2>cig_d & agegrp==0
replace trace22="typeconv cigd->smkd cigd" if cig_d!=. & cig_d!=0 &
smk_d
== . & smk_d_est2 != . & smk_c==. & smk_d_est2<=cig_d & agegrp==0
replace smk_d=cig_d if cig_d!=. & cig_d!=0 &
smk_d
== . & smk_d_est2 != . & smk_c==. & smk_d_est2<=cig_d & agegrp==0

cap append using `finaltemp'
save `finaltemp', replace
}
invalid syntax
r(198);

OUTPUT
.
. ****Regression****
.
. set more off

. set matsize 5000

Current memory allocation

current memory usage
settable value description (1M = 1024k)
--------------------------------------------------------------------
set maxvar 5000 max. variables allowed 1.733M
set memory 500M max. data space 500.000M
set matsize 5000 max. RHS vars in models 191.154M
-----------
692.887M

.
. use `allsvys1', clear

. egen i=group(un_region2 sex_code) if excind!=1 & sex_code!=3 &
total_age_grou
>p!=1 & agegrp==0
(9973 missing values generated)

. summarize i

Variable | Obs Mean Std. Dev. Min Max
-------------+--------------------------------------------------------
i | 0

. local max1=r(max)

. gen error1=.
(9973 missing values generated)

. gen error2=.
(9973 missing values generated)

. save `allsvys', replace
(note: file C:\DOCUME~1\khannag\LOCALS~1\Temp\ST_0000001b.tmp not found)
file C:\DOCUME~1\khannag\LOCALS~1\Temp\ST_0000001b.tmp saved

.
.
. forvalues svy = 1/`max1' {
2. use `allsvys', clear
3. keep if i==`svy'
4.
. summarize smk_d
5. local maxsd=r(max)
6. if `maxsd'>=95 local maxsd=95
7. local minsd=r(min)
8. if `minsd'<=5 local minsd=5
9.
. summarize cig_d
10. local maxcd=r(max)
11. if `maxcd'>=95 local maxcd=95
12. local mincd=r(min)
13. if `mincd'<=5 local mincd=5
14.
. replace lgsmkd=ln((smk_d-(`minsd'-5))/((`maxsd'+5)-smk_d))
15. replace lgcigd=ln((cig_d-(`mincd'-5))/((`maxcd'+5)-cig_d))
16.
. replace agesmkd=midage*lgsmkd
17. replace agecigd=midage*lgcigd
18.
. list un_region2 in 1
19. list sex_code in 1
20.
. sw regress lgcigd lgsmkd midage agesmkd [aweight=sample_size] if
exci
>nd!=1 & sex_code!=3 & total_age_group!=1 & totalind!=1 & agegrp==0 &
>source_c
>ode!=101692 & source_code!=101731 & source_code!=101712 &
>source_code!=101014
> & source_code!=101711 & source_code!=101724, pr(0.05) lockterm1 hier
21. mark type if excind!=1 & sex_code!=3 & total_age_group!=1 &
totali
>nd!=1 & agegrp==0 & lgcigd!=. & lgsmkd!=. & source_code!=101014 &
>source_code
>!=101711 & source_code!=101724 & source_code!=101692 & source_code!=101731
>& source_code!=101712
22. predict yhat1
23. replace r2noc=e(r2_a)
24. replace cig_d_est2 =
((`mincd'-5)+(`maxcd'+5)*exp(yhat1))/(1+exp(y
>hat1))
25. replace cig_d=min(cig_d_est2, cig_c, smk_d) if cig_d == . &
cig_d
>_est2 != . & agegrp==0
26.
.
****************************************************************************
*
>*****************************************************
.
****************************************************************************
*
>*****************************************************
.
. sw regress lgsmkd lgcigd midage agecigd [aweight=sample_size] if
type
>==1, pr(0.05) lockterm1 hier
27. predict yhat3
28. replace r2nos=e(r2_a)
29. replace smk_d_est2 =
((`minsd'-5)+(`maxsd'+5)*exp(yhat3))/(1+exp(y
>hat3))
30.
. replace trace22="typeconv cigd->smkd ok" if cig_d!=. & cig_d!=0

&
>smk_d == . & smk_d_est2 != . & smk_c==. & smk_d_est2>cig_d & agegrp==0
31. replace smk_d=smk_d_est2 if cig_d!=. &
cig_d!=0
> & smk_d == . & smk_d_est2 != . & smk_c==. & smk_d_est2>cig_d & agegrp==0
32. replace trace22="typeconv cigd->smkd cigd" if cig_d!=. &
cig_d!=0
> & smk_d == . & smk_d_est2 != . & smk_c==. & smk_d_est2<=cig_d & agegrp==0
33. replace smk_d=cig_d if cig_d!=. &
cig_d!=0
> & smk_d == . & smk_d_est2 != . & smk_c==. & smk_d_est2<=cig_d & agegrp==0
34.
. cap append using `finaltemp'
35. save `finaltemp', replace
36. }
invalid syntax
r(198);

Thank you.

Gauri

_________________________________________________________________
Don't just search. Find. Check out the new MSN Search!
http://search.msn.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/
_________________________________________________________________
Express yourself instantly with MSN Messenger! Download today it's FREE! http://messenger.msn.click-url.com/go/onm00200471ave/direct/01/

*
* 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