网状meta分析的meta回归，代码出错求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

参考张天嵩老师发表的“基于频率学框架的网状Meta回归实现与应用”的文章，使用相同的数据和代码，但是出现“could not calculate numerical derivatives，flat or discontinuous region encountered” 这个错误。
原代码：
foreach trt in A B C D E F G H I{
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 G H I{
if "`trt2'"=="A" continue
if "`trt2'">"`trt'" gen S`trt'`trt2' = 1/rA + 1/(nA-rA) if !mi(r`trt') & !mi(r`trt2')
}
}
format y* S* %6.2g
mat P8=0.5*(I(8)+J(8,8,1))

sum duration
gen x=duration-r(mean)
mvmeta y S x, bscov(prop P8) pbest(max,all zero gen(probm) reps(10000)) eform


求助！！万分感谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：实现与应用 与应用 原代码

“could not calculate numerical derivatives，flat or discontinuous region encountered”

数据不好
无法找到解

蓝色 发表于 2018-7-1 10:05
“could not calculate numerical derivatives，flat or discontinuous region encountered”

数据不好 ...

你好，数据不好是单纯从计算方面解释吗？但是别人做出来了。。能有其他的解决办法吗？