我对meta回归有疑问。我正在使用变量执行metareg:'sumvariable',其值可以是(0-6)。我想看看这个变量对死亡率的影响,表示为RD或OR。 这些是我生成的代码:
gen or = (death1/nondeath1)/(death2/nondeath2)
gen logor = log(or)
gen selogor = sqrt((1/death1)+(1/nondeath1)+(1/death2)+(1/nondeath2))
gen rd = (death1/Npat1) - (death2/Npat2)
gen risk1= death1/Npat1
gen risk2= death2/Npat2
gen serd= sqrt(risk1*(1-risk1)/Npat1) + (risk2*(1-risk2)/Npat2)
当我为RD执行metareg时:“metareg rd sumvariable,wsse(serd)”我收到错误:零或负wsse()不允许。所以我在没有以下代码的情况下进行了新的metareg分析:
gen risk1= death1/Npat1
gen risk2= death2/Npat2
gen serd= sqrt(risk1*(1-risk1)/Npat1) + (risk2*(1-risk2)/Npat2)
然后metareg会给出一个输出,但是现在它具有与我执行带有OR作为结果的metareg时相同的标准错误(参见输出)。当我使用相同的变量执行线性回归时,RD和OR之间的标准误差不同(您会期望的)。为什么metareg为OR和RD提供相同的SE?
Stata输出我得到了类似的SE: 在一些研究中,由于零细胞,RD和OR之间的观察数量不同 当我使用CVA而不是死亡率(使用以下代码生成)执行metareg时,metareg确实为RD和OR提供了2个不同的SE:
gen or = (cerebro1/noncerebro1)/(cerebro2/noncerebro2)
gen logor = log(or)
gen selogor = sqrt((1/cerebro1)+(1/noncerebro1)+(1/cerebro2)+(1/noncerebro2))
gen rd = (cerebro1/Npat1) - (cerebro2/Npat2)
CVA的STATA输出: