Dear Statalist
I really need your help to correct and run code. I'm currently working on meta-analysis (correlation coefficients) using a random-effects approach. I would have to Fisher's z-transform r.
*I am using stata 16/SE 16.
My operation code was as follow :
clear
import excel file//
ssc install admetan
generate z = atanh(r)
generate sez = sqrt(1/(n - 3))
metan z sez, label(namevar = author, yearvar = year)
admetan z sez
generate _USE = 1
generate lb = tanh(_LCI)
generate ub = tanh(_UCI)
generate _LABELS = author + " (" + string(year, "%02.0f") + ")"
label var n "Sample size"
local new = _N + 1
set obs `new'
replace _LABELS = "{bf:Overall}" if _n == _N
replace r = tanh(r(eff)) if _LABELS == "{bf:Overall}"
replace lb = tanh(r(eff) - (1.96 * r(se_eff))) if _LABELS == "{bf:Overall}"
replace ub = tanh(r(eff) + (1.96 * r(se_eff))) if _LABELS == "{bf:Overall}"
replace _USE = 5 if _LABELS == "{bf:Overall}"
format %-9.0g n
forestplot r lb ub, nonull effect("Correlation") rcol(n) leftjustify nowt
Then, I got a nice forest plot, however, the forest plot used fixed-effect inverse-variance model NOT random-effects modal.
My questions are,
- Is there perhaps any other Stata command for meta-analyses that would do Random-effect modal for Fisher’s transformed correlation coefficient?
- Can I have a correlation forest plot for Subgroup Analysis?
I really appreciated your help and thanks in advance.
Best wishes,
Areej
Comment