st: RE: RE: Leave one out

["Scott Merryman" <smerryman@kc.rr.com>]

Nick Cox kindly pointed out to me that -regress- automatically returns the
root mean square error, so this could be rewritten as:

sysuse auto,clear
count
tempvar rmse
gen `rmse' = .
qui {
forv i = 1/`=r(N)' {
	reg mpg gear if _n != `i'
	replace `rmse' = e(rmse) in `i'
	}
 sum `rmse'
}	

disp "RMSE (LOOCV) = "r(mean)

Scott





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