网状meta回归求助,恳请懂得人解答
按照15年张天嵩老师发表的文章《基于频率学框架的网状Meta回归实现与应用》进行网状meta回归,但是在运行运行过程中出现以下报错,恳请各位大神们帮助
第一步
foreach trt in A B C D E F {
if " `trt' "=="A" continue
gen y`trt' = log(r`trt'/(n`trt'-r`trt')) - log(rA/(nA-rA))
gen S`trt'`trt' = 1/r`trt' + 1/(n`trt'-r`trt') + 1/rA + 1/(nA-rA)
foreach trt2 in A B C D E F {
if " `trt2' "=="A" continue
if " `trt2' ">" `trt' " gen S`trt'`trt2' = 1/rA + 1/(nA-rA) if !mi(r`trt') & !mi(r`trt2')
}
}
运行后出现
(5 missing values generated)
(6 missing values generated)
(6 missing values generated)
(5 missing values generated)
(6 missing values generated)
variable yB already defined
r(110);
最后一步
mvmeta y S x, bscov(prop P6) pbest(max,all zero gen(probm) reps(10000)) eform
运行后出现
bscovariance() option: P6 is not 2x2
r(198);
以下是我的数据