原程序代码:
usehttp://nicolaorsini.altervista.org/data/bmi_rc, clear
bysort id: gen bmic = bmi - bmi[1]
mkspline bmics = bmic, nk(3) cubic
mkspline bmis = bmi, nk(3) cubic
mvmeta_makeglst logor bmics1 bmics2, cov(n case) se(selogor) pfirst(id studyt) saving(ssest_spline)replace by(id) names(b V)
preserve
use ssest_spline, clear
mvmeta b V , mm i2
capture estimates savemvmeta, replace
restore
estimates use mvmeta
su bmis1 bmis2 if bmi == float(19.3)
predictnl logor_sp =_b[bbmics1]*(bmis1-19.3) + _b[bbmics2]*(bmis2-0) , ci(los his)



雷达卡







京公网安备 11010802022788号







