capture log close 
log using 771bs2, replace
set more off
set seed 123456789
local reps 1000

capture program drop muratio
program define muratio, rclass
	version 8.0
	syntax varlist(min=2 max=2)
	tempname ymu
	summarize `1', meanonly
	scalar `ymu' = r(mean)
	summarize `2', meanonly
	return scalar ratio = `ymu'/r(mean)
end

webuse auto, clear
tabstat price, by(foreign) stat(n mean semean)
g p_dom = price if foreign==0
g p_for = price if foreign==1
muratio p_dom p_for
return list

bootstrap "muratio p_dom p_for" r(ratio), reps(`reps') saving(771bs2) replace

use 771bs2,clear
graph drop _all
histogram _bs_1, normal name(g771bs2_r) xsize(9) ysize(7) ti("Bootstrap distribution of ratio of means")
graph export 771bs2.pdf, replace
log close
log2html 771bs2.smcl, replace