[David Fisher]

Hi Suhail,
I think the metan option extraline(yes) should do what you need.
Thanks,
David.

[Suhail Doi]

Hi David,

Another small suggestion for a future update - for the forestplot, if the het stats are by default on the next line after overall then it makes formatting easier. Currently they are on the same line and need adding a column before the CI plot to move them to the next line

Regards
Suhail

[Suhail Doi]

Hi David,

Many thanks for the update. Something has changed with the "range" option in forestplot (only for proportions). For example if I run the below, it used to work perfectly in the previous version but now I cannot use range to get a well fitted plot no matter what I try and even if I remove the denom option

input str3 studyname long n int cases byte qi
"S1"  217154 422 1
"S10"  16557  32 1
"S13"    676   1 1
"S18"     44   1 1
"S26"     29   1 1
end


metan cases n  , pr qe(qi) transform(ftukey, iv) study( studyname ) counts forestplot(astext(85) boxscale(35) spacing(0.99) leftjustify range(0.1 20) dp(3)) denom(1000)

Please advise
Regards
Suhail

[David Fisher]

Luis Furuya-Kanamori : Unfortunately, upon opening this thread to post the above, I realise that I completely missed your post #16. Really sorry about that; this will be top of the list for the next update! In the meantime, here is a work-around which I hope is acceptable:

preserve
    metan mild n, pr model(re \ 0.50 0.40 0.60) transform(ftukey, iv) counts clear nograph
    foreach x in ES LCI UCI {
        replace _Prop_`x' = _`x' if _LABELS=="Overall, User"
    }
    forestplot, useopts
restore

Transformed proportion data is currently handled as a unique case in metan, because the Freeman-Tukey transform and back-transform are complex. The exponential/log transform is simple, so I only need to store the log-transformed values and exponentiate them as and when needed. But with transformed proportion data I need to maintain two different sets of effect sizes. The bug you have found is a result of metan looking in the wrong place for the user-defined values. The code above copies these values and pastes them in the correct place.

Apologies again, and best wishes,

David.

[David Fisher]

Dear all,

With thanks as ever to Kit Baum, an updated version of metan (v4.03 28apr2021) is now available via SSC.

The following bugs have (hopefully) been fixed:

• cumulative option with proportions (post #12 in this thread)
• ipdmetan with aggregated data (post #13 in this thread)
• incorrect Mantel-Haenszel Risk Ratios in the presence of zero cells (post #52 in https://www.statalist.org/forums/for...lot-from-metan)
• use of lalign as a plot option in Stata 15.0 (https://www.statalist.org/forums/for...gn-not-allowed)

Thanks to all users who have reported issues; please continue to do so!

Best wishes,

David.


P.S. Please note that, as a result of this update, users of ipdmetan will also need to update that package (via SSC). The latest version of ipdmetan is now v4.02 20apr2021.

[Suhail Doi]

Luis Furuya-Kanamori and David Fisher I don't think we should open up user defined tau2 for other models especially QE or IVhet as they only work with DL values and thus would open this up to wrong implementation and then another series of problems....

Regards
Suhail

[Luis Furuya-Kanamori]

Hi David Fisher,
I think I found a bug - when I do a proportion meta-analysis using the options 'model' and 'user_spec', in the forest plot, the diamond does not appear and ‘(Insufficient data)’ appears instead of the pooled estimate.
These options seem to work fine with standard (non-proportion) meta-analysis.

Attached is the dataset and the command I used

metan mild n, pr model(re \ 0.50 0.40 0.60) transform(ftukey, iv) counts

Also, would it be possible to active the option 'sa [, isq(real) tausq(real)]' for non-standard models or modifications to standard models (at the moment it only works with standard inverse-variance random-effects model).

Regards,
Luis

[David Fisher]

Nicoletta Riva Thanks for the feedback, great to hear!

Carole Khairallah :
Thanks for the bug report. I think I have identified and fixed it (privately), but of course I am only working with my own data and not with yours. Also, I am still unsure what is going on when you say "I could get estimates from that particular study before". Therefore, could you please do one of the following things for me?
(1) Send me (via Private Message) a copy of your working data and code which reproduces the error. Anonymised and/or sampled if necessary, just so long as the error is present.
or (2) Send me a Private Message with a contact email address for you, so that I can send you a corrected version of the program which you could run on your data and feed-back on whether it works or not.

Many thanks, and best wishes,

David.

[Nicoletta Riva]

Thanks a lot David!!

metan [varlist], proportion transform(ftukey, iv) works great!

Best regards,

Nicoletta

[Carole Khairallah]

Thank you for the update.
However, since I updated my ado files, I get a warning when using ipdmetan with aggregated data from 1 study, while I could get estimates from that particular study before.

Note: 1 study with missing or insufficient data found;  use the keepall option to include this study in the summary table and forest plot

Also, adding the option "keepall" to show the study that was excluded from the model returns an error

>1 invalid name
Error in ipdmetan.metan_setup

Is there a way to run an earlier version of ipdmetan? It may not be related to the update, but I can't think of anything else...I am using Stata v16.1.
Please see below the syntax I use for ipdmetan

 ipdmetan, study(study) ad(Data\temp\nigeria_bin21.dta, byad vars(_ES _lci _uci)) saving(Data\temp\bin21.dta, replace): glm bin21 arm01 i.Grav3Gr, fam(bin) link(log) eform

Would you have any suggestion? Thanks in advance for your help.
Best regards,
Carole

[Suhail Doi]

Hi David,

Just noted that the cumulative command on the proportions forest plot in metan and on the table output only changes the weights but not the rest of the plot or the table output.

Regards
Suhail

Edit: It does work as intended if we omit the tr(xyz , x) command

[Luis Furuya-Kanamori]

Many thanks David, it is working well now

Luis

[David Fisher]

With thanks as ever to Kit Baum, an updated version of metan (v4.02)is now available via SSC.

The bugs identified by Nicoletta Riva and by Luis Furuya-Kanamori have now (hopefully!) been fixed.

Many thanks to daniel klein for suggesting a solution to the first bug and for help with testing under Stata 11.2 and 12.1, and to Suhail Doi for helping me to track down the second bug.

Additionally, the method described in post #5 is now implemented as: metan [varlist], proportion transform(ftukey, iv)

Best wishes to all,

David.

[Luis Furuya-Kanamori]

Hi David,
Thanks for updating metan.
I tried to use the Freeman-Tukey double-arcsine transformation and there may be a glitch.

use http://fmwww.bc.edu/repec/bocode/m/metan_example_data, clear
metan cdeath csample, proportion transform(ftukey, iv) label(namevar=id)

gives the error message " 'ivariance' found where number expected "


Kind regards,
Luis

[Nicoletta Riva]

Hi David,

Thanks for your prompt reply.

Temporarily, I went back to the previous commands (admetan and metaprop, respectively) and the forest plots are created correctly. I'm looking forward to see the update.

Thanks again for your help!!

Nicoletta