R/ModelIPDmetareg.R

In htx-r/GenericModelNMA:

modelIPDmetareg <- function(){

  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]
    }
  }

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

  # tretament arm: random effect
  for (i in 1:(nIPD) ) {
    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,)

  # interaction coefficient
  beta_w~dnorm(0,0.001)
  beta_b~dnorm(0,0.001)

  # covariate coefficient
  beta0~dnorm(0,0.001)
}
#