Learn R Programming

mvmeta (version 1.0.3)

berkey98: Five Published Trials on Periodontal Disease

Description

The dataset contains the results of 5 published trials comparing surgical and non-surgical treatments for medium-severity periodontal disease, one year after treatment. The 2 estimated outcomes are average improvements (surgical minus non-surgical, in mm) in probing depth (PD) and attachment level (AL).

Usage

berkey98

Arguments

Format

A data frame with 5 observations on the following 7 variables:

pubyear

publication year of the trial.

npat

number of patients included in the trial.

PD

estimated improvement of surgical versus non-surgical treatments in probing depth (mm).

AL

estimated improvement of surgical versus non-surgical treatments in attachment level (mm).

var_PD

variance of the estimated outcome for PD.

cov_PD_AL

covariance of the estimated outcomes for PD and AL.

var_AL

variance of the estimated outcome for AL.

Row names specify the author of the paper reporting the results of each trial.

Examples

Run this code
# NOT RUN {
### REPRODUCE THE RESULTS IN BERKEY ET AL. (1998)

# INSPECT THE DATA
berkey98

# FIXED-EFFECTS
year <- berkey98$pubyear - 1983
model <- mvmeta(cbind(PD,AL)~year,S=berkey98[5:7],data=berkey98,method="fixed")
print(summary(model),digits=3)

# GLS MODEL (VARIANCE COMPONENTS) 
model <- mvmeta(cbind(PD,AL)~year,S=berkey98[5:7],data=berkey98,method="vc",
  control=list(vc.adj=FALSE))
print(summary(model),digits=3)
round(model$Psi,3)

# ML MODEL
model <- mvmeta(cbind(PD,AL)~year,S=berkey98[5:7],data=berkey98,method="ml")
print(summary(model),digits=3)
round(model$Psi,3)
# }

Run the code above in your browser using DataLab