Learn R Programming

mixmeta (version 1.2.0)

thrombolytic: Randomized Trials of Thrombolytic Therapy

Description

The dataset contains the data on 20 randomized trials of thrombolytic therapy, which evaluated effect on short-term mortality after a myocardial infarction (up to 35 days) in 50,246 patients in relation to treatment delay. The hypothesis is that the thrombolytic therapy reduces the mortality risk following the myocardial infarction, and that the benefit is particularly substantial for very early treatment. Some of the trials report separate results according to treatment delay, generating 38 observations from full trials or subgroups of trials. Effect sizes were reported as absolute risk reduction computed as the difference between treated and control groups in each trial or subgroup.

Usage

thrombolytic

Arguments

Format

A data frame with 38 observations on the following 10 variables:

  • trial: label identifying the trial.

  • time2treat: treatment delay after the onset of the symptoms of a myocardial infaction, reported in hours.

  • dtreat, ntreat: number of deaths and total patients in the treated group, respectively.

  • dcontr, ncontr: number of deaths and total patients in the control group, respectively.

  • risktreat, riskcontr: risk of death in the treatment and control groups, respectively.

  • absrisk: absolute risk difference of death between the treatment and control groups. See Details.

  • var: variance of the absolute risk difference. See Details.

Details

The absolute risk is simply the difference in risk, which is computed empirically as ratio of the number of deaths and the number of total patients in treated and control groups (\(p_1=d_1/N_1\) and \(p_0=d_0/N_0\), respectively). The variance of the absolute risk difference is computed as \(p_0(1-p_0)/N_0+p_1(1-p_1)/N_1\). See Thompson and colleagues (2001) for details.

Examples

Run this code
# NOT RUN {
### REPRODUCE THE RESULTS IN THOMPSON ET AL (2001), TABLES 2, 3, AND 4

# STANDARD FIXED-EFFECTS META-ANALYSIS
mod1 <- mixmeta(absrisk, var, data=thrombolytic, method="fixed")
print(summary(mod1), digits=5)

# STANDARD RANDOM-EFFECTS META-ANALYSIS
subtrial <- seq(nrow(thrombolytic))
mod2 <- mixmeta(absrisk, var, random= ~ 1|subtrial, data=thrombolytic)
print(summary(mod2), digits=5)

# TWO-LEVEL RANDOM-EFFECTS META-ANALYSIS
mod3 <- mixmeta(absrisk, var, random= ~ 1|trial/subtrial, data=thrombolytic)
print(summary(mod3), digits=5)

# TWO-LEVEL RANDOM-EFFECTS META-REGRESSION
mod4 <- mixmeta(absrisk~time2treat, var, random= ~ 1|trial/subtrial,
  data=thrombolytic)
print(summary(mod4), digits=5)

# TWO-LEVEL RANDOM-EFFECTS META-REGRESSION WITH NON-LINEAR TERM
mod5 <- mixmeta(absrisk ~ time2treat + I(1/time2treat), var,
  random= ~ 1|trial/subtrial, data=thrombolytic)
print(summary(mod5), digits=5)

### SEE help(school) FOR A COMPLEMENTARY EXAMPLE
# }

Run the code above in your browser using DataLab