doseresNMA antidep/random.R

In htx-r/doseresNMA:

# data
source('data.prep.R') # prepare data
antidep <- data.prep()

# number of treatments in the network: all agent-specific doses (330+ 1= placebo)
antidep%>%
  group_by(drug)%>%
  group_map(~length(unique(.x$dose)))%>%
  unlist %>%
  sum

# study with the largest sample size 
antidep[which.max(antidep$n),]