풀링 유병률 메타분석 R로 하는 방법

Pooled prevalence meta-analysis with R

이미 Systematic review 방식으로 논문을 모아놓은 것을 가정하고.

다음과 같이 csv 파일로 데이터를 정리한다. (data extraction)

RStudio에서 우선 meta 패키지를 로딩하고

library(meta)

Session –> Set Working Directory로 들어가서 파일이 저장된 경로를 설정해준다.

CSV 파일을 읽어들인다.

df1 <- read.csv(“meta.csv”, header=T)

single proportion data를 분석할 것이므로 metaprop() 함수를 이용한다.

meta1 <- metaprop(diseased.n, total.n, sm=”PLN”, data=df1, studlab=paste(study,year), comb.fixed=F)

  • random effect model을 이용할 것이므로 comb.fixed = False로 설정했다. (heterogeneity가 높음. 연구마다 유병률이 꽤 다르고 시대적 특징 등 인구집단에 차이가 있음)
  • logit normalization을 이용하므로 sm=”PLN”으로 설정.

이렇게 하면 meta1에는 메타분석 결과가 정리된다.

논문에 쓸 forest plot도 간단하게 뽑아낼 수 있다.

forest(meta1)

보고하는 내용

  • The pooled prevalence was 23% (95% CI, 11-47%) with high between-study heterogeneity (Q=240, τ^2=0.54, I^2=99%, P<0.01).

 

More reads:

CC BY-NC-SA 4.0 This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License.

Leave a Comment

Your email address will not be published. Required fields are marked *

Time limit is exhausted. Please reload CAPTCHA.