[R-meta] [SPAM](FN) Re: Is this R script correct for conducting a meta-analysis of AOR for hypertension prevalence?

[Dr. Gerta Rücker]

Hi Sancho and Michael,

There is clearly an error in Sancho's code. The first screenshot shows that he declared sm = "AOR" - which doesn't make sense. The possible values of sm (summary measure) are from a fixed list, e.g., "RD", "RR", "OR", "ASD", "HR", "MD", "SMD", or "ROM". Consequently, as "AOR" is not a possible value, it is ignored. 

Now, as no summary measure is specified, metagen() assumes the mean difference (MD) to be the effect measure. This also explains why the confidence intervals (that include 1) correspond to highly significant p-values (p < 0.0001) in the second screenshot: As  implicitly mean differences are assumed, the null hypothesis assumes a mean of 0, not 1. And the CIs don't include 0.

Solution: You have to set sm = "OR".

Best,
Gerta




UNIVERSITÄTSKLINIKUM FREIBURG
Institute for Medical Biometry and Statistics

Dr. Gerta Rücker
Guest Scientist

Stefan-Meier-Straße 26 · 79104 Freiburg
gerta.ruecker using uniklinik-freiburg.de

https://www.uniklinik-freiburg.de/imbi-en/employees.html?imbiuser=ruecker

-----Ursprüngliche Nachricht-----
Von: R-sig-meta-analysis <r-sig-meta-analysis-bounces using r-project.org> Im Auftrag von Michael Dewey via R-sig-meta-analysis
Gesendet: Mittwoch, 31. Juli 2024 18:30
An: Sancho Xavier <sanchoxavierxavier using gmail.com>; R Special Interest Group for Meta-Analysis <r-sig-meta-analysis using r-project.org>
Cc: Michael Dewey <lists using dewey.myzen.co.uk>
Betreff: [SPAM](FN) Re: [R-meta] Is this R script correct for conducting a meta-analysis of AOR for hypertension prevalence?

Dear Sancho

Responses in line

On 31/07/2024 15:27, Sancho Xavier wrote:
> Dear Michael,
> 
> Thank you for your response. I did not apply any moderator for the 
> pooled OR in the meta-analysis. Could you clarify something for me?
> In a meta-analysis, is it always necessary to calculate the pooled OR 
> including a moderator variable?

No, you include a moderator if your scientific hypothesis demands it.

> *Another question:*
> Was the pooled OR supposed to be significant in the following conde (image):

Well the confidence interval clearly includes the null. I am not sure 
what you mean by supposed to be here.

If you want to include output it is better to paste it in text format as 
including screenshots mangles the post.

Michael

> 
> 
> *Results:*
> 
> *Dataset:*
> *
> *
> **
> *
> *
> Do you have any suggestions, ideas, or observations about the scripts 
> and results based on the dataset?
> 
> Best wishes,
> 
> *Sancho Xavier*
> Candidate for Postgraduate Programme in Epidemiology in Public Health
> MBA in Data Science and Analytics at USP (ongoing)
> *Linkedin:* linkedin.com/in/sancho-xavier-59565a18a
> *Orcid iD:* https://orcid.org/0000-0001-9493-4098
> *Web of Science (Clarivate):* 
> https://www.webofscience.com/wos/author/record/GON-5675-2022 
> <https://www.webofscience.com/wos/author/record/GON-5675-2022>
> *Reference Citation Analysis (RCA): 
> *https://referencecitationanalysis.com/00001471 
> <https://referencecitationanalysis.com/00001471>
> *Whatsapp: *+258 852255846
> ------------------------------------------------------------------------------
> ------------------------------------------------------------------------
> *De:* Michael Dewey <lists using dewey.myzen.co.uk>
> *Enviado:* quarta-feira, 31 de julho de 2024 05:05
> *Para:* R Special Interest Group for Meta-Analysis 
> <r-sig-meta-analysis using r-project.org>
> *Cc:* Sancho Xavier <sanchoxavierxavier using gmail.com>
> *Assunto:* Re: [R-meta] Is this R script correct for conducting a 
> meta-analysis of AOR for hypertension prevalence?
> Dear Sancho
> 
> Are you intending to use the variable you have named variable as a
> moderator? If not the code looks fine to me although I do not usually
> use that package (I assume it is meta).
> 
> Michael
> 
> On 27/07/2024 00:18, Sancho Xavier via R-sig-meta-analysis wrote:
>> *Hi, community!*
>> I am conducting a meta-analysis and meta-regression to calculate the 
>> */pooled adjusted odds ratio (AOR)/* and */I² with the p-value/* for 
>> variables associated with hypertension prevalence. I have a dataset with 
>> the following columns:
>> 
>>   * *study_id:* Study identifier
>>   * *AOR:* Adjusted odds ratio (each independent variable)
>>   * *lower_CI:* Lower confidence interval (each independent variable)
>>   * *upper_CI:* Upper confidence interval (each independent variable)
>>   * *variable:* Independent variable
>> 
>> Here is the script I am using in *R*:
>> 
>> *# Calculate the standard error (seTE)*
>> *
>> *
>> data$seTE <- (log(data$upper_CI) - log(data$lower_CI)) / (2 * 1.96)
>> 
>> *# Meta-analysis*
>> *
>> *
>> meta_analysis <- metagen(TE = log(data$AOR),
>> seTE = data$seTE,
>> data = data,
>> studlab = data$study_id,
>> comb.fixed = TRUE,
>> comb.random = TRUE)
>> 
>> *I would like to know if this script is correct for achieving my goals of:*
>> 
>>   * Calculating the pooled AOR and the 95% confidence interval.
>>   *
>>     Calculating the I² and the p-value to assess heterogeneity.
>> 
>> Look at the dataset:
>> 
>> Any feedback or suggestions for improvement would be greatly appreciated!
>> 
>> Thank you!
>> 
>> *Sancho Xavier*
>> Candidate for Postgraduate Programme in Epidemiology in Public Health
>> MBA in Data Science and Analytics at USP (ongoing)
>> *Linkedin:* linkedin.com/in/sancho-xavier-59565a18a
>> *Orcid iD:* https://orcid.org/0000-0001-9493-4098 
> <https://orcid.org/0000-0001-9493-4098>
>> *Web of Science (Clarivate):* 
>> https://www.webofscience.com/wos/author/record/GON-5675-2022 
> <https://www.webofscience.com/wos/author/record/GON-5675-2022>
>> <https://www.webofscience.com/wos/author/record/GON-5675-2022 
> <https://www.webofscience.com/wos/author/record/GON-5675-2022>>
>> *Reference Citation Analysis (RCA): 
>> *https://referencecitationanalysis.com/00001471 
>> <https://referencecitationanalysis.com/00001471 
> <https://referencecitationanalysis.com/00001471>>
>> *Whatsapp: *+258 852255846
>> ------------------------------------------------------------------------------
>> 
>> <http://www.avg.com/email-signature?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=emailclient <http://www.avg.com/email-signature?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=emailclient>>  Virus-free.www.avg.com <http://www.avg.com/email-signature?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=emailclient <http://www.avg.com/email-signature?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=emailclient>>
>> 
>> <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
>> 
>> _______________________________________________
>> R-sig-meta-analysis mailing list @ R-sig-meta-analysis using r-project.org
>> To manage your subscription to this mailing list, go to:
>> https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis 
> <https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis>
> 
> -- 
> Michael

-- 
Michael

_______________________________________________
R-sig-meta-analysis mailing list @ R-sig-meta-analysis using r-project.org
To manage your subscription to this mailing list, go to:
https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis