Learn R Programming

mvmeta (version 1.0.3)

smoking: Meta-Analysis of Interventions to Promote Smoking Cessation

Description

The dataset contains the results of 24 trials comparing four alternative interventions to promote smoking cessation. The trials have different designs, comparing two or three different interventions. The data consist of the number of successes out of the total participants, and the estimated log-odds ratio for arms B, C, and D relative to arm A, as well as the (co)variance matrix of these three estimates.

Usage

smoking

Arguments

Format

A data frame with 24 observations on the following 19 variables:

  • study: study ID.

  • design: design of the trial, reporting the interventions being compared.

  • dA, dB, dC, dD: number of successes for each intervention.

  • nA, nB, nC, nD: number of participants for each intervention.

  • yB, yC, yD: estimated log-odds ratios for interventions B, C and D versus intervention A.

  • SBB, SBC, SBD, SCC, SCD, SDD: variances and co-variances of the estimated log-odds ratios for interventions B, C and D versus intervention A. The order corresponds to the lower triangular elements of the (co)variance matrix taken by column.

Details

Intervention A is chosen as the reference category. Trials without an arm A (trials 2 and 21-24) are augmented with an arm A with 0.01 individuals and 0.001 successes. Trials containing zero cells (trials 9 and 20) have 1 individual with 0.5 successes added to each intervention. Details on the data augmentation and estimation of (co)variances of the log-odds ratios are provided by White (2011).

Examples

Run this code
# NOT RUN {
### REPRODUCE THE RESULTS IN WHITE (2011)

# INSPECT THE DATA
head(smoking)
names(smoking)

# UNSTRUCTURED BETWEEN-STUDY (CO)VARIANCE
y <- as.matrix(smoking[11:13])
S <- as.matrix(smoking[14:19])
model <- mvmeta(y,S)
summary(model)

# STRUCTURED BETWEEN-STUDY (CO)VARIANCE (PROPORTIONAL)
model <- mvmeta(y,S,bscov="prop",control=list(Psifix=diag(3)+1))
summary(model)

# SEE help(mvmetaCovStruct) for additional info and examples
# }

Run the code above in your browser using DataLab