## FIRST EXAMPLE: Single case-control study
## Linear trend estimation
## Inspect data
data("cc_ex")
## Fitting the model
mod1 <- dosresmeta(formula = logrr ~ dose, type = "cc", cases = case,
n = n, lb = lb, ub = ub, data= cc_ex)
summary(mod1)
## Results
predict(mod1, delta = 1)
## SECOND EXAMPLE: Multiple studies
## Linear and quadratic trend using random-effects meta-analysis
## Inspect data
data("alcohol_cvd")
## Linear trend
lin <- dosresmeta(formula = logrr ~ dose, type = type, id = id,
se = se, cases = cases, n = n, data = alcohol_cvd)
## Summarize the results
summary(lin)
predict(lin, delta = 1)
## Non-linear (quadratic) trend
quadr <- dosresmeta(formula = logrr ~ dose + I(dose^2), type = type, id = id,
se = se, cases = cases, n = n, data = alcohol_cvd)
## Summarize the results
summary(quadr)
## Graphical results
with(predict(quadr), {
plot(dose, pred, log = "y", type = "l",
xlim = c(0, 45), ylim = c(.4, 2))
lines(dose, ci.lb, lty = 2)
lines(dose, ci.ub, lty = 2)
rug(dose, quiet = TRUE)
})
Run the code above in your browser using DataLab