运行mvmeta命令时，遇到“file mvmeta.ster not found r(601)”的问题-经管之家官网！

各位站友，我在学习stata11的Worked example of glst command时遇到一下问题：
不知道是因为什么原因，在“capture estimates save mvmeta, replace”这句命令中，这个.ster文件始终无法保存，以至于在"estimates use mvmeta" 命令里出现“file mvmeta.ster not found r(601)”的问题.
请问这是哪里出问题了，怎么解决啊？
谢谢了！
一下附上这个.do文件里面的命令：
version 11
use http://nicolaorsini.altervista.org/data/bmi_rc, clear
// Center each study to the baseline reference dose (first row within each study)
// NB: This is important to avoid bias estimates
// because glst is assuming the fitted line is going through the origin (lnRR=0,dose=0)
bysort id: gen bmic = bmi - bmi[1]
// create spline transformations for the centered exposure
mkspline bmics = bmic, nk(3) cubic
// create spline transformations for the original exposure to be used for predictions
mkspline bmis = bmi, nk(3) cubic
// First stage - Study specific estimates using splines
mvmeta_make glst logor bmics1 bmics2, cov(n case) se(selogor) pfirst(id studyt) saving(ssest_spline) replaceby(id)names(b V)
// Second stage - Pool study specific estimates using random-effects
preserve
use ssest_spline, clear
mvmeta b V, mm i2
capture estimates save mvmeta, replace
restore
// Obtain and plot the dose-response relationship back on the original scale
estimates use mvmeta
// Get values of the splines at the chosen reference value (min)
su bmis1 bmis2 if bmi == float(19.3)
predictnl logor_sp = _b[bbmics1]*(bmis1-19.3) + _b[bbmics2]*(bmis2-0) , ci(los his)
gen ors = exp(logor_sp)
gen lbs = exp(los)
gen ubs = exp(his)
// Compare with a linear trend using a two-stage random-effect
glst logor bmic, cov(n case) se(selogor) pfirst(id studyt) ts(r)
// Obtain and plot the dose-response relationship back on the original scale
predictnl logor_l= _b[bmic]*(bmi-19.3)
gen orl= exp(logor_l)
// Overlay predictions from spline and linear models
twoway///
(line orl bmi, sortlc(black) lp(-) ) ///
(line orslbs ubsbmi, sort ///
lw(thick thick thick) lc(black black black) ///
lp(l longdash longdash)) , ///
yscale(log)///
ytitle("Relative Risk") xtitle("Body Mass Index") ///
scheme(s1mono) ///
xlabel(19(1)35) xmtick(19(.5)35) ///
ymtick(1(.5)8) ylabel(1 2 4 8 , format(%3.2fc) angle(horiz)) ///
legend(label(1 "Linear Model")label(2 "Spline Model") ///
order(1 2) ring(0) pos(11) col(1) ) ///
plotregion(style(none))