library("mvmeta")
library("dosresmeta")
lin <- dosresmeta(formula = logHR ~ dose, id = study,
type = type, se = se, cases = cases,
n = personyear, data = breast)
如果(delta〈tol)中断时出错:需要TRUE/FALSE时缺少值否则:警告信息:对数(Ax [v!= 0]):生成NaN
1条答案
按热度按时间bihw5rsg1#
breast$cases
中至少有一个值是负数。您可以使用以下命令进行检查:dosresmeta
函数将调用covar.logrr
函数,covar.logrr
函数将调用grl
函数(在得到错误后尝试运行traceback()
)。grl
中的基本代码为:实际上,该函数获取
cases
列,将其赋给Ax
,然后尝试运行log(Ax[v!=0])
(请参见您的错误)。我们知道,取负数的对数将得到
NaN
,如您的错误所述:因此,我认为您的
cases
列包含非正数。