Hi,
I note that admetan has included the IVhet and quality effects models of meta-analysis which is great for the community of researchers. However a concern I have is that the use of the DL tau squared within these methods does not correspond to any actual probability distribution and therefore is not replaceable. However admetan does allow this to happen. For example take the dataset below:
The IVhet model can be run correctly as:
or wrongly as
I would consider that the second code above is no longer the IVhet model
Similarly, the quality effects model can be run correctly as:
or wrongly as:
Again the second is no longer the quality effects model
Any thoughts on this would be welcome
I note that admetan has included the IVhet and quality effects models of meta-analysis which is great for the community of researchers. However a concern I have is that the use of the DL tau squared within these methods does not correspond to any actual probability distribution and therefore is not replaceable. However admetan does allow this to happen. For example take the dataset below:
Code:
* Example generated by -dataex-. To install: ssc install dataex clear input str11 studyname float(es lo95ci hi95ci qi rescaled lnes lnlo lnhi) "Leenders" .99 .98 1 8 1 -.010050327 -.020202687 0 "Nagura" .91 .86 .95 8 1 -.09431065 -.1508229 -.05129331 "Sahyoun" .95 .83 1.08 7 .875 -.05129331 -.1863296 .07696108 "Strandhagen" .9 .81 .99 5 .625 -.10536055 -.21072103 -.010050327 "Tucker" .93 .84 1.03 6 .75 -.07257068 -.1743534 .029558774 "Whiteman" .84 .72 .99 6 .75 -.1743534 -.328504 -.010050327 "Zhang" .94 .91 .97 7.5 .9375 -.06187541 -.09431065 -.03045918 end
Code:
admetan lnes lnlo lnhi, re(ivhet) rr forestplot(astext(70) boxscale(50) spacing(1) leftjustify lcols(studyname) dp(2))
Code:
admetan lnes lnlo lnhi, re(reml,ivhet) rr forestplot(astext(70) boxscale(50) spacing(1) leftjustify lcols(studyname) dp(2))
Similarly, the quality effects model can be run correctly as:
Code:
admetan lnes lnlo lnhi, qe(rescaled) rr forestplot(astext(70) boxscale(50) spacing(1) leftjustify lcols(studyname) dp(2))
Code:
admetan lnes lnlo lnhi, re(reml,qe(rescaled)) rr forestplot(astext(70) boxscale(50) spacing(1) leftjustify lcols(studyname) dp(2))
Any thoughts on this would be welcome
Comment