I am using the command 'metatrim' in stata.
I am doing a meta-analysis of effective reproduction number. The LFK index shows a major asymmetry and publication bias. So, I try to trim and fill in the effect size.
But the result states: "No trimming was performed: data unchanged". I am not sure why it is not working. Can anyone please help?
Code:
metatrim rt lci uci, reffect funnel
I expect to get some trimming and filling. However, in the result, state states:
Code:
Note: no trimming performed; data unchanged
For some weird reason, I can reproduce my result using some fake data.
Here I am just manipulating the automiss data, to have an effect size- mpg with its lower and upper 95% confidence lintervals- lci and uci
Code:
webuse automiss.dta, clear keep make mpg rep78 replace rep78 = 2 if rep78 == . generate lci = mpg - rep78 generate uci = mpg + rep78 drop rep78 drop if mpg == .
package st0061 from http://www.stata-journal.com/software/sj4-2
Code:
net sj 4-2 st0061
package sbe19_6 from http://www.stata-journal.com/software/sj9-2
Code:
net sj 9-2 sbe19_6
Code:
metafunnel mpg lci uci, ci name(funnel_ci, replace) // funnel plot metabias mpg se, egger //
Code:
generate se = (uci-lci)/ 3.92 generate logmpg = log(mpg) generate loglci = log(lci) generate loguci = log(uci) generate logse = (loguci-loglci)/ 3.92 admetan logmpg loglci loguci, model(reml) effect("mpg") eform
Code:
net describe lfk, from(http://fmwww.bc.edu/RePEc/bocode/l)
Code:
lfk _ES _seES
sbe39_2 from http://www.stata.com/stb/stb61
Code:
net describe sbe39_2, from(http://www.stata.com/stb/stb61)
Code:
metatrim mpg lci uci, ci funnel reffect
Putting in all the same codes as above together for easier copy-pasting in stata:
Code:
webuse automiss.dta, clear keep make mpg rep78 replace rep78 = 2 if rep78 == . generate lci = mpg - rep78 generate uci = mpg + rep78 drop rep78 drop if mpg == . generate se = (uci-lci)/ 3.92 generate logmpg = log(mpg) generate loglci = log(lci) generate loguci = log(uci) generate logse = (loguci-loglci)/ 3.92 metafunnel mpg lci uci, ci name(funnel_ci, replace) // funnel plot metabias mpg se, egger admetan logmpg loglci loguci, model(reml) effect("mpg") eform lfk _ES _seES metatrim mpg lci uci, ci funnel reffect
0 Response to Trim and fill using metatrim
Post a Comment