Hello, I am using the metan command to do a meta-analysis but now I need to display p-values for overall effect (to the right) and the forest plot I am getting only shows p value for heterogeneity (to the left).
Is there a way or code to still have the heterogeneity p-value on the left-hand side of the plot, and also add the p-value for overall effect on the right-hand side underneath the overall OR?
I am using: gen logor= log((ExpOa/NexpOC)/(ExpNOb/NExPNOD))
gen selogor=sqrt((1/ExpOa)+(1/ExpNOb)+(1/NexpOC)+(1/NExPNOD))
metan logor selogor, eform fixed label (namevar=NAME, yearvar=YEAR) xlabel (0.1, 0.2, 0.5, 1, 2, 3, 5, 7, 10) force effect("Odds ratio")
metan ExpOa ExpNOb NexpOC NExPNOD, rr lcols (NAME YEAR) xlabel (0.1, 0.2, 0.5, 1, 2, 3, 5, 7, 10) title (¨SDNN MT¨) xtitle (¨Risk ratio¨) textsize (175) boxsca (80) notable xsize(10) ysize (6)
to get this forest plot:
Is there a way or code to still have the heterogeneity p-value on the left-hand side of the plot, and also add the p-value for overall effect on the right-hand side underneath the overall OR?
I am using: gen logor= log((ExpOa/NexpOC)/(ExpNOb/NExPNOD))
gen selogor=sqrt((1/ExpOa)+(1/ExpNOb)+(1/NexpOC)+(1/NExPNOD))
metan logor selogor, eform fixed label (namevar=NAME, yearvar=YEAR) xlabel (0.1, 0.2, 0.5, 1, 2, 3, 5, 7, 10) force effect("Odds ratio")
metan ExpOa ExpNOb NexpOC NExPNOD, rr lcols (NAME YEAR) xlabel (0.1, 0.2, 0.5, 1, 2, 3, 5, 7, 10) title (¨SDNN MT¨) xtitle (¨Risk ratio¨) textsize (175) boxsca (80) notable xsize(10) ysize (6)
to get this forest plot:
