st: Re: monte carlo simulation on panel data

["Joseph Coveney" <stajc2@gmail.com>]

Ghazal Zulfiqar wrote:

I need help with a Monte Carlo simulation on panel data with 500 observations
(the panel is organized by year and country_code). I need to run a regression
and then collect fitted values in 1000 runs, so that I have a new dataset of
500,000. The program I wrote only gives me 1000 observations, that is, it is
giving me one statistic per run. Can someone please help?

set seed 3213799
gen random=.
gen diversity=.	
tempname sim
postfile `sim' xb using results, replace
quietly {
forvalues i=1/1000 {
replace random=runiform()
replace diversity=median(conc_diversity)+random
regress stddev diversity asia latin_am europe africa_mideast log_gdp_capita
log_gdp_capita_1 log_gdp_capita_2 pop_g_rate barro_ipol investment_rate t t_sqrd
predict xb
post `sim' (xb) 
drop xb
}
}
postclose `sim'
summarize

--------------------------------------------------------------------------------

After -postclose `sim'-, and before -summarize-, you need to
-use `sim', clear-

Joseph Coveney


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