R/ModelIPDADmetareg.R

In htx-r/GenericModelNMA:

modelIPDADmetareg <-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]+beta0*x[i,j]+beta_w*(x[i,j]*treat[i,j]-xbar[i])+beta_b*xbar[i]
    }
  }

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


  # Model for AD
  for (i in 1:nAD) {
    rc[i] ~ dbin(pc[i],nc[i])
    rt[i] ~ dbin(pt[i],nt[i])
    logit(pc[i]) <- u.a[i]
    logit(pt[i]) <- u.a[i]+delta[i+nIPD]+beta_b*xbar.a[i]
  }

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

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

  # prior
  mu~dnorm(0,0.001)
  inv.tau.sq <- 1/(tau.sq)
  tau.sq <- tau*tau
  tau~dnorm(0,1)%_%T(0,)

  beta_w~dnorm(0,0.001) # interaction coefficients within trial
  beta_b~dnorm(0,0.001) # interaction coefficients between trial
  beta0~dnorm(0,0.001)  # covariate effect coefficient


}