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:
- A basic introduction to fixed-effect and random-effects models for meta-analysis
- Package ‘meta’ documentation
- Meta-analysis of prevalence
More from my site
This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License.
Leave a Comment