Hi All, I performed a random-effect summary meta-analysis where study was used as a random effect. Here, I used the rma.mv function using the Fisher's r-to-z transformed correlation coefficient effect size (see code below). CODE: results=rma.mv(yi,vi,data=datum,method="REML",random=~1|Study) I can use then use the below code to create funnel and forest plots. CODE: forest(results) CODE: funnel(results) The x axis for these plots would be my effect size (Fisher's r-to-z transformed correlation coefficient). My question is, what must I do in order to back transform my Fisher's r-to-z transformed correlation coefficient here in the plots so that they are displayed as the Pearson's correlation coefficient (r). Thanks, John
[R-meta] Back-transforming Fisher's r-to-z transformed correlation coefficient to r in funnel and forest plots results
4 messages · John Mahas, Wolfgang Viechtbauer
Hi All, Please disregard my previous message - I figured it out. I isolated my generated yi and vi from my dataset and ran the following code. #load z results datum=read.csv(file.choose()) Z=datum #transform z to r Z_r=transf.ztor(Z) Once this was done, I added this data to a csv to perform my analyses with the back-transformed r instead of my original Fisher's r-to-z transformed correlation coefficient
From: R-sig-meta-analysis <r-sig-meta-analysis-bounces at r-project.org> on behalf of John Mahas via R-sig-meta-analysis <r-sig-meta-analysis at r-project.org>
Sent: Wednesday, August 16, 2023 2:36 PM
To: R Special Interest Group for Meta-Analysis <r-sig-meta-analysis at r-project.org>
Cc: John Mahas <jwm0055 at auburn.edu>
Subject: [EXT] [R-meta] Back-transforming Fisher's r-to-z transformed correlation coefficient to r in funnel and forest plots results
Sent: Wednesday, August 16, 2023 2:36 PM
To: R Special Interest Group for Meta-Analysis <r-sig-meta-analysis at r-project.org>
Cc: John Mahas <jwm0055 at auburn.edu>
Subject: [EXT] [R-meta] Back-transforming Fisher's r-to-z transformed correlation coefficient to r in funnel and forest plots results
CAUTION: Email Originated Outside of Auburn.
Hi All,
I performed a random-effect summary meta-analysis where study was used as a random effect. Here, I used the rma.mv function using the Fisher's r-to-z transformed correlation coefficient effect size (see code below).
CODE: results=rma.mv(yi,vi,data=datum,method="REML",random=~1|Study)
I can use then use the below code to create funnel and forest plots.
CODE: forest(results)
CODE: funnel(results)
The x axis for these plots would be my effect size (Fisher's r-to-z transformed correlation coefficient). My question is, what must I do in order to back transform my Fisher's r-to-z transformed correlation coefficient here in the plots so that they are displayed as the Pearson's correlation coefficient (r).
Thanks,
John
[[alternative HTML version deleted]]
_______________________________________________
R-sig-meta-analysis mailing list @ R-sig-meta-analysis at r-project.org
To manage your subscription to this mailing list, go to:
https://nam11.safelinks.protection.outlook.com/?url=https%3A%2F%2Fstat.ethz.ch%2Fmailman%2Flistinfo%2Fr-sig-meta-analysis&data=05%7C01%7Cjwm0055%40auburn.edu%7C6e8d806eb6fc425e0c6008db9e901c35%7Cccb6deedbd294b388979d72780f62d3b%7C1%7C0%7C638278114034637278%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=6tPP73Ww0yRIqnF7zU7cXA1dx8jSZr4K5iLfOa9YTio%3D&reserved=0<https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis>
4 days later
Dear John, That would be one way of doing this, but then your analysis uses raw correlation coefficients. The forest() and funnel() functions have the 'atransf' arguments for back-transforming the x-axes. Best, Wolfgang
-----Original Message----- From: R-sig-meta-analysis [mailto:r-sig-meta-analysis-bounces at r-project.org] On Behalf Of John Mahas via R-sig-meta-analysis Sent: Wednesday, 16 August, 2023 22:05 To: R Special Interest Group for Meta-Analysis Cc: John Mahas Subject: Re: [R-meta] Back-transforming Fisher's r-to-z transformed correlation coefficient to r in funnel and forest plots results Hi All, Please disregard my previous message - I figured it out. I isolated my generated yi and vi from my dataset and ran the following code. #load z results datum=read.csv(file.choose()) Z=datum #transform z to r Z_r=transf.ztor(Z) Once this was done, I added this data to a csv to perform my analyses with the back-transformed r instead of my original Fisher's r-to-z transformed correlation coefficient
________________________________ From: R-sig-meta-analysis <r-sig-meta-analysis-bounces at r-project.org> on behalf of John Mahas via R-sig-meta-analysis <r-sig-meta-analysis at r-project.org> Sent: Wednesday, August 16, 2023 2:36 PM To: R Special Interest Group for Meta-Analysis <r-sig-meta-analysis at r- project.org> Cc: John Mahas <jwm0055 at auburn.edu> Subject: [EXT] [R-meta] Back-transforming Fisher's r-to-z transformed correlation coefficient to r in funnel and forest plots results Hi All, I performed a random-effect summary meta-analysis where study was used as a random effect. Here, I used the rma.mv function using the Fisher's r-to-z transformed correlation coefficient effect size (see code below). CODE: results=rma.mv(yi,vi,data=datum,method="REML",random=~1|Study) I can use then use the below code to create funnel and forest plots. CODE: forest(results) CODE: funnel(results) The x axis for these plots would be my effect size (Fisher's r-to-z transformed correlation coefficient). My question is, what must I do in order to back transform my Fisher's r-to-z transformed correlation coefficient here in the plots so that they are displayed as the Pearson's correlation coefficient (r). Thanks, John
Hi Wolfgang, Thank you for pointing this out! I also realized that issue soon after performing the analysis. I did not know about the "atransf" argument, but I have incorporated this into my analysis and was able to successfully back-transform the x-axes. Thanks, John
From: Viechtbauer, Wolfgang (NP) <wolfgang.viechtbauer at maastrichtuniversity.nl>
Sent: Monday, August 21, 2023 2:58 AM
To: R Special Interest Group for Meta-Analysis <r-sig-meta-analysis at r-project.org>
Cc: John Mahas <jwm0055 at auburn.edu>
Subject: [EXT] RE: Back-transforming Fisher's r-to-z transformed correlation coefficient to r in funnel and forest plots results
Sent: Monday, August 21, 2023 2:58 AM
To: R Special Interest Group for Meta-Analysis <r-sig-meta-analysis at r-project.org>
Cc: John Mahas <jwm0055 at auburn.edu>
Subject: [EXT] RE: Back-transforming Fisher's r-to-z transformed correlation coefficient to r in funnel and forest plots results
CAUTION: Email Originated Outside of Auburn. Dear John, That would be one way of doing this, but then your analysis uses raw correlation coefficients. The forest() and funnel() functions have the 'atransf' arguments for back-transforming the x-axes. Best, Wolfgang >-----Original Message----- >From: R-sig-meta-analysis [mailto:r-sig-meta-analysis-bounces at r-project.org] On >Behalf Of John Mahas via R-sig-meta-analysis >Sent: Wednesday, 16 August, 2023 22:05 >To: R Special Interest Group for Meta-Analysis >Cc: John Mahas >Subject: Re: [R-meta] Back-transforming Fisher's r-to-z transformed correlation >coefficient to r in funnel and forest plots results > >Hi All, > >Please disregard my previous message - I figured it out. I isolated my generated >yi and vi from my dataset and ran the following code. > >#load z results >datum=read.csv(file.choose()) >Z=datum > >#transform z to r >Z_r=transf.ztor(Z) > >Once this was done, I added this data to a csv to perform my analyses with the >back-transformed r instead of my original Fisher's r-to-z transformed correlation >coefficient > >________________________________ >From: R-sig-meta-analysis <r-sig-meta-analysis-bounces at r-project.org> on behalf >of John Mahas via R-sig-meta-analysis <r-sig-meta-analysis at r-project.org> >Sent: Wednesday, August 16, 2023 2:36 PM >To: R Special Interest Group for Meta-Analysis <r-sig-meta-analysis at r- >project.org> >Cc: John Mahas <jwm0055 at auburn.edu> >Subject: [EXT] [R-meta] Back-transforming Fisher's r-to-z transformed correlation >coefficient to r in funnel and forest plots results > >Hi All, > >I performed a random-effect summary meta-analysis where study was used as a >random effect. Here, I used the rma.mv function using the Fisher's r-to-z >transformed correlation coefficient effect size (see code below). > >CODE: results=rma.mv(yi,vi,data=datum,method="REML",random=~1|Study) > >I can use then use the below code to create funnel and forest plots. > >CODE: forest(results) >CODE: funnel(results) > >The x axis for these plots would be my effect size (Fisher's r-to-z transformed >correlation coefficient). My question is, what must I do in order to back >transform my Fisher's r-to-z transformed correlation coefficient here in the >plots so that they are displayed as the Pearson's correlation coefficient (r). > >Thanks, > >John