Fwd: st: loop does not load file as stated

[Jin-Hyun Bae <jinhyun.bae@stern.nyu.edu>]

I did as you said, but Stata didn't fail on any line. It just ran
through. Still have the same problem though.

Looking at the results, the only weird thing happening is that the
third loop isn't working. It isn't producing estimates (or outputting
them)...

But, the code I have been using for that loop is equivalent to the
other code doing essentially the same thing but with different
variables.
----------------------

set a read-only flag on the data file in the OS, and rerun. Then see
the line which fails in Stata.

On Fri, Aug 2, 2013 at 12:49 PM, Jin-Hyun Bae <jinhyun.bae@stern.nyu.edu> wrote:
> Hello!
>
> I have two levels of loops below. The top level cycles through years.
> Within this loop, the file Distance.dta is loaded, and then all drops
> all irrelevant year observations are dropped. When the second level
> loops have finished running, it should save this file as "year".dta
> and clear all.
>
> I think every time that the first loop goes back to the beginning the
> Distance.dta should be the same because I have not asked for any
> changes to be saved to Distance.dta. But when I run this code, it
> looks like all the observations that were dropped in the earlier loop
> remained dropped.
>
> I am almost certain that there is a simple explanation for this but I
> am completely baffled!!
>
> Jin
>
> --------------
>
> forvalues y = 1970/2004{
>         use Distance.dta, clear
>         keep if Year == `y'
>         gsort -usa
>         local N = _N
>         matrix V1 = (POLCON[1] \ WBgovexp[1] \ WBdomcred[1])
>         matrix rownames V1 = "POLCON WBgovexp WBdomcred"
>         forvalues i = 2/`N'{
>                 matrix V`i' = (POLCON[`i'] \ WBgovexp[`i'] \ WBdomcred[`i'])
>                 matrix rownames V`i' = "POLCON WBgovexp WBdomcred"
>                 mahascore2 POLCON WBgovexp WBdomcred, point1(V1) point2(V`i')
> covarpop(one) compute
>                 replace poldist = r(mahascore_sq) in `i'
>                 }
>                 matrix drop _all
>
>         matrix V1 = (WBpcGDP[1] \ Wbinfln[1] \ WBX[1] \ WBM[1])
>         matrix rownames V1 = "WBpcGDP Wbinfln WBX WBM"
>         forvalues i = 2/`N'{
>                 matrix V`i' = (WBpcGDP[`i'] \ Wbinfln[`i'] \ WBX[`i'] \ WBM[`i'])
>                 matrix rownames V`i' = "WBpcGDP Wbinfln WBX WBM"
>                 mahascore2 WBpcGDP Wbinfln WBX WBM, point1(V1) point2(V`i')
> covarpop(one) compute
>                 replace econdist = r(mahascore_sq) in `i'
>                 }
>         matrix drop _all
>
>         matrix V1 = (Power[1] \ Uncert[1] \ Indlm[1] \ Masc[1] \ Relig[1])
>         matrix rownames V1 = "Power Uncert Indlm Masc Relig"
>         forvalues i = 2/`N'{
>                 matrix V`i' = (Power[`i'] \ Uncert[`i'] \ Indlm[`i'] \ Masc[`i'] \ Relig[`i'])
>                 matrix rownames V`i' = "Power Uncert Indlm Masc Relig"
>                 mahascore2 Power Uncert Indlm Masc Relig, point1(V1) point2(V`i')
> covarpop(one) compute
>                 replace cultdist = r(mahascore_sq) in `i'
>                 }
>         matrix drop _all
>         saveold `y'.dta, replace
> }
> *
> *   For searches and help try:
> *   http://www.stata.com/help.cgi?search
> *   http://www.stata.com/support/faqs/resources/statalist-faq/
> *   http://www.ats.ucla.edu/stat/stata/
*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/faqs/resources/statalist-faq/
*   http://www.ats.ucla.edu/stat/stata/
*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/faqs/resources/statalist-faq/
*   http://www.ats.ucla.edu/stat/stata/