# Use a subset of 3 studies from anesthesia data
anesthesia_subset <-
subset(anesthesia, study %in% unique(anesthesia$study)[1:3])
# Prepare data for DR-NMA
dat <- pairwise(
agent = list(agent1, agent2, agent3),
event = list(event1, event2, event3),
n = list(n1, n2, n3),
dose = list(dose1, dose2, dose3),
data = anesthesia_subset,
studlab = study,
append = FALSE)
# DR-NMA with linear dose-response function
dr1 <- netdose(TE, seTE, agent1, dose1, agent2,
dose2, studlab,
data = dat)
# Dose-response plot
dotplot(dr1)
Run the code above in your browser using DataLab