FunnelPlotR
FunnelPlotR copied to clipboard
Allow choice of tau2 estimation method
An interesting additional feature could be the ability to use REML estimation to obtain the tau2
term, which could also allow for additional levels of variance terms for inflating the control limits.
A possible workflow could be:
- User calls
tau2_method="REML"
in thefunnel_plot
call (where thetau2_method
option would default toDL
) - Construct unadjusted z-scores
- Truncate data based on unadjusted z-scores
- Pass truncated data to
metafor::rma.mv
, using the transformed outcome and standard errors - Extract estimated variance
Alternatively, metafor
's escalc
function could be used directly for constructing the standard errors.
Let me know what you think, happy to work this up and open a PR
That's a really good option. I'll think on that. Happy to collaborate on it, or receive a PR. Thanks!