Dear Statalist, I have recently done a meta-analysis of the correlation coefficient, but I encountered such an error in the process of analyzing a set of data with stata software,I tried it many times, but it failed. Please help me correct the error and rewrite the running code. Thank you!
My operation process was as follows:
First I run the code
input str10 study year r n
Bill 2002 .34 39
Amanda 2004 .41 80
Kevin 2004 .55 40
Geroge 2006 .45 200
Jacob 2008 .38 98
lin 2009 .35 100
Fang 2010 .24 58
end
generate z = atanh(r)
generate sez = sqrt(1/(n - 3))
metan z sez, label(namevar = study, yearvar = year)
admetan z sez
generate _USE = 1
generate lb = tanh(_LCI)
generate ub = tanh(_UCI)
generate _LABELS = study + " (" + 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) nowt
Then I got the result like this (as shown below. )
--------------------------------------------------------------------
Study | Effect [95% Conf. Interval] % Weight
---------------------+----------------------------------------------
1 | 0.354 0.027 0.681 6.06
2 | 0.436 0.212 0.659 12.96
3 | 0.618 0.296 0.941 6.23
4 | 0.485 0.345 0.624 33.16
5 | 0.400 0.199 0.601 15.99
6 | 0.365 0.166 0.564 16.33
7 | 0.245 -0.020 0.509 9.26
---------------------+----------------------------------------------
Overall effect | 0.424 0.343 0.504 100.00
--------------------------------------------------------------------
Test of overall effect = 0: z = 10.322 p = 0.000
Heterogeneity Measures
---------------------------------------------------------
| Value df p-value
---------------------+-----------------------------------
Cochran's Q | 4.46 6 0.614
I?(%) | 0.0%
Modified H? | 0.000
tau? | 0.0000
---------------------------------------------------------
I?= between-study variance (tau? as a percentage of total variance
Modified H?= ratio of tau?to typical within-study variance
option lalign() not allowed
Error in forestplot
r(198);
Related Posts with How to draw a correlation forest map with the admetan command?
Doubt of the Cox modelHello everyone, I am doing my thesis which seeks to find the variables that influence the survival …
stcrreg competing risk model with random interceptDear Stata community, I have run the Fine-Gray competing risk model using stcrreg and worked perfec…
about xtqreg commandhello stata community, Can the xtqreg command be applied to a sample of 12 banks over 13 quarters (…
Q about generating survey weight in combined NHANES cycleHello everyone. I have a question about creating/generating NHANES survey weight for the 2011-2012 a…
Parsing a string variable to extract a particular wordHello everyone, I need to extract the word INC or CORPORATION or INCORPORATED from the string varia…
Subscribe to:
Post Comments (Atom)
0 Response to How to draw a correlation forest map with the admetan command?
Post a Comment