dat.lloyd2010.RdResults from 16 studies evaluating anti-TNF-alpha inhibitors in patients with rheumatoid arthritis.
dat.lloyd2010The data frame contains the following columns:
| author | character | first author |
| year | integer | publication year |
| mean.das | numeric | mean for outcome DAS-28 |
| lower.das | numeric | lower limit for outcome DAS-28 |
| upper.das | numeric | upper limit for outcome DAS-28 |
| mean.haq | numeric | mean for outcome HAQ |
| lower.haq | numeric | lower limit for outcome HAQ |
| upper.haq | numeric | upper limit for outcome HAQ |
| n | integer | sample size |
Lloyd et al. (2010) report results of a systematic review evaluating the effectiveness of anti-TNF-alpha inhibitors in the treatment of rheumatoid arthritis. The authors conducted separate meta-analyses for HAQ and DAS-28.
Lloyd, S., Bujkiewicz, S., Wailoo, A. J., Sutton, A. J., & Scott, D. (2010). The effectiveness of anti-TNF-\(\alpha\) therapies when used sequentially in rheumatoid arthritis patients: A systematic review and meta-analysis. Rheumatology, 49(12), 2313–2321. https://doi.org/10.1093/rheumatology/keq169
medicine, rheumatology, raw mean differences
### Show first five studies
head(dat.lloyd2010, 5)
#> author year mean.das lower.das upper.das mean.haq lower.haq upper.haq n
#> 1 Bennet 2005 -1.70 -2.11 -1.29 -0.31 -0.57 -0.05 26
#> 2 Bingham 2009 -1.60 -1.78 -1.42 -0.35 -0.44 -0.26 188
#> 3 Bombardieri 2007 -1.90 -2.00 -1.80 -0.48 -0.52 -0.44 810
#> 4 Buch 2007 -1.47 -1.79 -1.15 NA NA NA 72
#> 5 Cohen 2005 -1.87 -2.26 -1.48 NA NA NA 30
### Load meta package
suppressPackageStartupMessages(library(meta))
### Only consider studies providing data for both outcomes
lloyd5 <- subset(dat.lloyd2010, !is.na(mean.haq) & !is.na(mean.das))
### Univariate meta-analysis of the DAS-28 outcome
m.das <- metagen(mean.das,
lower = lower.das, upper = upper.das,
data = lloyd5, sm = "MD",
studlab = paste(author, year),
random = FALSE)
### Univariate meta-analysis of the HAQ outcome
m.haq <- metagen(mean.haq,
lower = lower.haq, upper = upper.haq,
data = lloyd5, sm = "MD",
studlab = paste(author, year),
random = FALSE)
### Forest plots
forest(m.das, test.overall = TRUE, hetstat = FALSE,
digits.TE = 2, digits.se = 2)
forest(m.haq, test.overall = TRUE, hetstat = FALSE,
digits.TE = 2, digits.se = 2)