Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.

# Re: st: Running loops with graphs

 From Sebastian Galarza To statalist@hsphsun2.harvard.edu Subject Re: st: Running loops with graphs Date Fri, 27 Apr 2012 09:38:01 -0400

```Eric,

Thanks for your suggestions. I will try to clarify my problem so that you or others can better address the issue. I currently have a panel database organized by oil wells (i) with monthly data (WOR) for 50 years.  I want to run the goodness of fit tests for each oil well just as you suggested but I am not sure where to include i unless it would be something along the lines of by `"`i'"' : nl …. however this is the reason I thought about using the for values command.

I hope this provides greater clarity and can serve others that encounter similar issues in the future.

Best,

Sebastian

On Apr 26, 2012, at 5:05 PM, Eric Booth wrote:

> <>
>
> You make reference to "i" without defining where it fits into your code, but assuming you know where "i" goes, here a (untested) code example that should do what you are asking with some adapting on your part.
>
> ******************
> forvalues i = 1/51000 {
>
> di `"`i'"' //where does `i' plug-in in the equations below??
> nl (WOR = {b0=0.1}*(1 * exp({b1=0.05}* (ageyear))))
> predict a
>
> nl log3: WOR ageyear
> predict b
>
> nl gom3: WOR ageyear
> predict c
>
> twoway (scatter a y) (line a y, sort), legend(off) aspect(1) name(g1, replace)
> twoway (scatter b y) (line b y, sort), legend(off) aspect(1) name(g2, replace)
> twoway (scatter c y) (line c y, sort), legend(off) aspect(1) name(g3, replace)
> gr combine g1 g2 g3, c(1) // change look of this combined graph with options
> gr export "graph`i'.eps", as(eps) replace
> }
> ******************
> ^ See -help forvalues- and -help graph combine- for more help.
>
>
>
> - Eric
>
> __
> Eric A. Booth
> Public Policy Research Institute
> Texas A&M University
> ebooth@ppri.tamu.edu
> +979.845.6754
>
>
>
>
>
> On Apr 26, 2012, at 12:18 PM, Sebastian Galarza wrote:
>
>> Hi,
>>
>> I am working on a large dataset and need to run different goodness of fit tests, graph these and save the graphs separately. I am not sure how to go about it.
>>
>> Basically, I want to run the following tests for up to 51000 values of i
>>
>> nl (WOR = {b0=0.1}*(1 * exp({b1=0.05}* (ageyear))))
>> predict a
>>
>> nl log3: WOR ageyear
>> predict b
>>
>> nl gom3: WOR ageyear
>> predict c
>>
>> I would ideally save this output for each value of i.
>>
>> twoway (scatter a y) (line a y, sort), legend(off) aspect(1)
>> twoway (scatter b y) (line b y, sort), legend(off) aspect(1)
>> twoway (scatter c y) (line c y, sort), legend(off) aspect(1)
>>
>> I would ideally save these three graphs as one image for each value of i.
>>
>> Any help in doing this would be greatly appreciated.
>>
>> Best,
>>
>> Sebastian
>>
>>
>> *
>> *   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/
>
>
> *
> *   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/

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