Ben Jann replied:
> Clive, the program should be:
>
> program define _estadd_rmse, eclass
> local rmse = sqrt(e(rss)/(e(N) - e(df_m) - e(df_a) -1))
> ereturn scalar rmse = `rmse'
> end
>
> The following would work too:
>
> program define _estadd_rmse, eclass
> ereturn scalar rmse = sqrt(e(rss)/(e(N) - e(df_m) - e(df_a) -1))
> end
Excellent! I ran the first formulation (successfully), and after running
-estadd- in the form you suggested, it worked a treat. Thanks!
CLIVE NICHOLAS |t: 0(044)7903 397793
Politics |e: [email protected]
Newcastle University |http://www.ncl.ac.uk/geps
*
* 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/
