mu$styles()              # define HTML styles, functions

# Function to generate data for one study
gendat <- function(study, n) {
f <- function(p) rbinom(n, 1, p)
d <-
  data.frame(study=study,
               subjid=1:n,
             site=sample(1:10, n, TRUE),
             tx=sample(c('A','B'), n, TRUE),
                         age=rnorm(n, 50, 10),
                         sex=sample(c('female', 'male'), n, TRUE),
                         edate=as.Date('2019-10-01') + round(runif(n, 0, 100)),
                         sbp=rnorm(n, 120, 7),
             dbp=rnorm(n,  80, 6),
             race=sample(c('Asian', 'Black/AA', 'White'), n, TRUE),
             meda=f(0.3),
             medb=f(0.2),
                         td=round(runif(n, 0, 30)),
                         death=rbinom(n, 1, 0.4),
                         e1     = f(.03),
             e2     = f(.03),
             e3     = f(.03),
             e4     = f(0.05),
             tested = rbinom(n, 1, .75),
             e6     = f(.03),
                 e7     = f(0.04)
                         )
d <-
  upData(d,
        e5     = ifelse(tested, f(.06), NA),    # don't randomize if exclusion present:
        rdate  = ifelse(e1|e2|e3|e4|ifelse(is.na(e5), 0, e5)|e6|e7, NA,
                    edate + round(runif(n, 0, 30))),
        rdate  = as.Date(rdate, origin='1970-01-01'),   # ifelse loses Date class
        tx     = ifelse(is.na(rdate), NA, tx),
        labels=c(age='Age',
                 rdate='Randomization Date',
                         edate='Enrollment Date',
                   e1='Prior MI', e2='History of Asthma',
             e3='History of Upper GI Bleeding',
             e4='No Significant CAD', e5='Inadequate Renal Function',
             e6='Pneumonia within 6 weeks', e7='Hostile',
             race='Race', sex='Sex', tx='treatment',
             sbp='Systolic BP', dbp='Diastolic BP',
             meda='Comcomitant medication A',
                         medb='Concomitant medication B',
                         td='Follow-up Time', death='Death'),
             units=c(sbp='mmHg', dbp='mmHg', age='years', td='day'), 
       print=FALSE)

denom <- c(enrolled=nrow(d), randomized=sum(! is.na(d$rdate)),
           table(subset(d, ! is.na(rdate))$tx))
sethreportOption(tx.var='tx', denom=denom, study=study)
d
}
set.seed(1)
d1 <- gendat('Study 1', 500)
d2 <- gendat('Study 2', 250)
d3 <- rbind(d1, d2)
denom <- c(enrolled=nrow(d3), randomized=sum(! is.na(d3$rdate)),
           table(subset(d3, ! is.na(rdate))$tx))
sethreportOption(tx.var='tx', denom=denom, study='Study 1+2')