vignettes/metalr_document.R

In LaminJuwara/metalr: Likelihood ratio meta-analysis

## ----setup, include = FALSE----------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

## ---- fig.show='hold'----------------------------------------------------
if (!require("devtools"))install.packages("devtools")
devtools::install_github("laminjuwara/metalr")

## ---- fig.show='hold'----------------------------------------------------
# Clopidogrel vs Aspirin trial dataset
cases<-c(939,1021)
person_yrs<-c(17636,17519)
patients<-c(9599,9586)   

## ---- fig.show='hold'----------------------------------------------------
library(metalr)
ici.rr(cases = cases,patients = patients,person_yrs = person_yrs)

## ---- fig.show='hold'----------------------------------------------------
library(metalr)
data("sample_metarr_data")
head(sample_metarr_data)

## ---- fig.show='hold'----------------------------------------------------
metalr_rr(sample_metarr_data,refval=0,num_iter=3000,increm=0.001,method="random")

## ---- fig.show='hold'----------------------------------------------------
library(metalr)
data("statindata")
statindata

## ---- fig.show='hold'----------------------------------------------------
ici.or(idata = statindata[6,2:5])

## ---- fig.show='hold'----------------------------------------------------
metalr_or(statindata[,2:5],refval=0,num_iter=3000,increm=0.001,method="fixed")

## ---- fig.show='hold',fig.align="center",fig.height=3.5, fig.width=5.5----
# the metalr object
library(metalr)
metalr_obj<-metalr_or(idata=statindata[,2:5],refval=0,num_iter=3000,increm=0.001,method = "random")

# forest plot of the metalr object
forest_metalr(metalr_obj)