R/ModelIPDADmeta.R

In htx-r/GenericModelNMA:

modelIPDADmeta <-function(){

  #** Model for IPD
  for (i in 1:nIPD) { # loop through studies
    for(j in 1: npIPD[i]){ # loop through participants

      # binomial likelihood
      y[i,j]  ~ dbern(p[i,j])
      # model linear predictor
      logit(p[i,j]) <- u[i] + delta[i]*treat[i,j]
    }
  }

  # IPD-control arm: common effect
  for (i in 1:nIPD) {
    u[i]~dnorm(0, 1e-6)
  }

  #** Model for AD
  for (i in 1:nAD) {
    # likelihood on both arms
    rc[i] ~ dbin(pc[i],nc[i])
    rt[i] ~ dbin(pt[i],nt[i])

    # logistic transformation on both arms
    logit(pc[i]) <- u.a[i]
    logit(pt[i]) <- u.a[i]+delta[i+nIPD]
  }

  # AD-control arm: common effect
  for (i in 1:nAD) {
    u.a[i] ~dnorm(0,0.001)
  }

  # random effect model to combine all treatment effect estimates
  for (i in 1:(nIPD+nAD) ) {
    delta[i]~dnorm(mu,inv.tau.sq)
  }

  # prior distributions
  mu~dnorm(0,1e-6)
  inv.tau.sq <- 1/(tau.sq)
  tau.sq <- tau*tau
  tau~dunif(0,0.1)
}