SAS Meta Analysis, Funnel plot

[mzfirat]

Hi all,

I have the following simple data set and I have been trying to do Meta analysis and obtain Funnel plot. I would be grateful if you could help me solve this problem.

Best regards

data mydata1 ;
input Record Subgroup $3-25 N_Gin MMP_Mean_Gin MMP_Std_Gin N_Healthy MMP_Mean_Healthy MMP_Std_Healthy;
datalines;
1 Yucel et al. 2020                20 609.77 174.13 23 625.74 163.1
2 Rai et al. 2008                   18 312 301.8 15 95.2 70.2
3 Ebersole et al. 2015         43 199 29.1 65 130.7 14.50
4 Syndergaard et al. 2014 40 208.2 194.2 40 156.7 121.9
5 Bostancı et al. 2021        31 272.4 403.2 36 76.5 53.7
6 Zhang et al. 2021 24 603.2 220.7 25 435.8 180.6
;
run;

[ballardw]

Which variable is used for which part of the plot? Any plot in two dimensions requires knowing which is an X axis value and which is a Y axis value. Other variables may be used to label or adjust display in some manner. But you haven't given us much to go on.

You should paste code such as a data step into a text box opened with the </> icon that appears above the message window. Your code as shown gets invalid data messages because there is not enough space, possibly because this forum re-formats text pasted into the main window, between the Subgroup and N_Gin values. So N_gin on the the last record becomes part of Subgroup and there are not enough values for the other variables and all of the ones after N_gin are in the wrong variable.

---

@mzfirat wrote:

Hi all,

I have the following simple data set and I have been trying to do Meta analysis and obtain Funnel plot. I would be grateful if you could help me solve this problem.

Best regards

data mydata1 ;
input Record Subgroup $3-25 N_Gin MMP_Mean_Gin MMP_Std_Gin N_Healthy MMP_Mean_Healthy MMP_Std_Healthy;
datalines;
1 Yucel et al. 2020                20 609.77 174.13 23 625.74 163.1
2 Rai et al. 2008                   18 312 301.8 15 95.2 70.2
3 Ebersole et al. 2015         43 199 29.1 65 130.7 14.50
4 Syndergaard et al. 2014 40 208.2 194.2 40 156.7 121.9
5 Bostancı et al. 2021        31 272.4 403.2 36 76.5 53.7
6 Zhang et al. 2021 24 603.2 220.7 25 435.8 180.6
;
run;

---

[mzfirat]

The meta-analysis is used to assess the impact by comparing the average concentration of 8-OHdG in gingivitis patients with that in a healthy control group. I want to assess the heterogeneity through the Q test's p-value and I2. I can easily obtain the Forest plots. But I need a  funnel plot depicting publication bias, and calculations including Egger's regression intercept, and its associated p-value. There will be mean difference on the X axis and standard error on the y axis in my funnel plot. 

proc format;
   value std_fmt
   1='Cross sectional'
   2='Case control';
run;
data mydata1 ;
input Record Subgroup $3-25 StudyType N_Gin	MMP_Mean_Gin	MMP_Std_Gin	N_Healthy	MMP_Mean_Healthy	MMP_Std_Healthy;
format StudyType std_fmt.;
datalines;
1 Yucel et al. 2020       1 20	609.77	174.13	23	625.74	163.1
2 Rai et al. 2008         1 18	312	    301.8	15	95.2	70.2
3 Ebersole et al. 2015    2 43	199	    29.1	65	130.7	14.50
4 Syndergaard et al. 2014 2 40	208.2	194.2	40	156.7	121.9
5 Bostancı et al. 2021    1 31	272.4	403.2	36	76.5	53.7
6 Zhang et al. 2021       1 24	603.2	220.7	25	435.8	180.6
;
run;

run;

>