MasterIPDADmeta.R

In htx-r/GenericModelNMA:

# libraries
library(R2jags)
library(meta)
library(devtools)
install_github('htx-r/GenericModelNMA',force = TRUE)
library(GenericModelNMA)
source('functions needed for master simulations IPDADmeta.R')

#**** IPDmeta
# 1. (twostage) freq=true
nsim <- 50
twostageIPDmetasim <-  replicate(nsim,twostageIPDmeta(reg = FALSE))
colMeans(t(twostageIPDmetasim))


# 2. (onestage) freq=true
onestageIPDmetasim <-  replicate(nsim,onestageIPDmeta(reg = FALSE))
colMeans(t(onestageIPDmetasim))

# 3. bayesian = true
nsim <- 2
bayesIPDmetasim <-  replicate(nsim,bayesIPDmeta(reg = FALSE))
colMeans(t(bayesIPDmetasim))


#**** IPDmetareg

# 1. (twostage) freq=true
nsim <- 50
twostageIPDmetaregsim <-  replicate(nsim,twostageIPDmeta(reg = TRUE))
colMeans(t(twostageIPDmetaregsim))


# 2. (onestage) freq=true
onestageIPDmetaregsim <-  replicate(nsim,onestageIPDmeta(reg = TRUE))
colMeans(t(onestageIPDmetaregsim))


# 3. bayesian=true
nsim <- 2
bayesIPDmetaregsim <-  replicate(nsim,bayesIPDmetareg(reg = TRUE))
colMeans(t(bayesIPDmetaregsim))
#

## IPD-AD meta
#  freq
nsim <- 100
freqIPDADmetasim <- replicate(nsim,freqIPDADmeta(reg = FALSE))
colMeans(t(freqIPDADmetasim))

# bayes
nsim=2
bayesIPDADmetasim <-  replicate(nsim,bayesIPDADmeta(reg = FALSE))
colMeans(t(bayesIPDADmetasim))

## IPD-AD metareg

# freq
nsim <- 100
freqIPDADmetaregsim <- replicate(nsim,freqIPDADmetareg(reg = TRUE))
colMeans(t(freqIPDADmetaregsim))
# 20.05.2020: it does not give the true results for delta and beta0

# bayes
nsim <- 2
bayesIPDADmetaregsim <-  replicate(nsim,bayesIPDADmetareg(reg = TRUE))
colMeans(t(bayesIPDADmetaregsim))