[R-meta] Fwd: How to add the regression line onto Deek's funnel plot

[Paul Chang]

To whom this may concern,

I've been learning how to perform meta-analysis of diagnostic test accuracy
(DTA).
I used mada package maintained by Prof. Doebler.
I'm aware that publication bias of DTA is better analyzed with the method
proposed by Deek et al in 2005.
The Deek's funnel plot plots the lnDOR against the inverse of the square
root of effective sample size (ESS).
Test of asymmetry is then conducted by regression of lnDOR with 1/root(ESS).
Please see the attachment "deeks2005" for more details.

Since Deek's funnel plot is not included in the mada package, and after
discussion with Prof. Doebler, I adopted his suggestion that I use funnel()
function within the metafor package to customize what I need to plot Deek's
funnel plot.
I've got so close except for the final step: adding the regression line to
the funnel plot.
I used the Child Behavior Checklist (CBCL) scale in the Chang et al study
to practice.
Please see the attachment for the dataset ("chang_cbcl.csv"), my code
("Deek funnel plot.R") and the results after running in the console
("console.txt").
The regression to test the asymmetry yielded same P-value as in the Chang
et al. study.
Please see the attachments of their funnel plot and mine
("chang_cbcl_deeksfunnel.png" and "myplot.png", respectively).

As I mentioned earlier, I couldn't find a way to correctly add the
regression line onto the funnel plot.
Is there any solution you suggest?
Please let me know.

Thank you very much. Best wishes.

-- 
張君毓
Chang, Chun-Yu (Paul)
慈濟大學醫學系100級
Class 2018, School of Medicine, Tzu Chi University
台北慈濟醫院PGY
Post-graduate-year doctor, Taipei Tzu Chi Hospital, Buddhist Tzu Chi
Medical Foundation, New Taipei, Taiwan
E-mail: paulchang1231 using gmail.com <paulchan1231 using gmail.com>
Cell: 0978000933

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://stat.ethz.ch/pipermail/r-sig-meta-analysis/attachments/20191023/efb2696a/attachment-0001.html>

-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: console.txt
URL: <https://stat.ethz.ch/pipermail/r-sig-meta-analysis/attachments/20191023/efb2696a/attachment-0001.txt>

-------------- next part --------------
A non-text attachment was scrubbed...
Name: chang_cbcl.csv
Type: text/csv
Size: 1408 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-meta-analysis/attachments/20191023/efb2696a/attachment-0001.csv>

-------------- next part --------------
A non-text attachment was scrubbed...
Name: Deek funnel plot.R
Type: application/octet-stream
Size: 1421 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-meta-analysis/attachments/20191023/efb2696a/attachment-0001.obj>

-------------- next part --------------
A non-text attachment was scrubbed...
Name: chang_cbcl_deeksfunnel.PNG
Type: image/png
Size: 31598 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-meta-analysis/attachments/20191023/efb2696a/attachment-0002.png>

-------------- next part --------------
A non-text attachment was scrubbed...
Name: myplot.png
Type: image/png
Size: 5566 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-meta-analysis/attachments/20191023/efb2696a/attachment-0003.png>

-------------- next part --------------
A non-text attachment was scrubbed...
Name: deeks2005.pdf
Type: application/pdf
Size: 2069940 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-meta-analysis/attachments/20191023/efb2696a/attachment-0001.pdf>