test/c212.BB/default_parameters/test/run.r

In rcarragh/c212: Methods for Detecting Safety Signals in Clinical Trials Using Body-Systems (System Organ Classes)

RNGkind(sample.kind="Rounding")
library(c212)

set.seed(5594)

trial_data <- read.table("../data/c212.trial.data.dat", header = T, as.is = T)

raw = c212.BB(trial_data)
conv = c212.convergence.diag(raw)
sink("conv.dat")
c212.print.convergence.summary(conv)
sink()

summ = c212.summary.stats(raw)
sink("summary.dat")
c212.print.summary.stats(summ)
sink()


ptheta = c212.ptheta(raw)
write.table(ptheta, "ptheta.dat")

ptheta80 = ptheta[ptheta$ptheta > 0.80,]
write.table(ptheta80, "ptheta80.dat")

# RAISED RATES
z <- order(ptheta80$AE)
ptheta80 <- ptheta80[z,]

write.table(ptheta80[, c("B", "AE")], "RAISED_AES", row.names = FALSE, quote = FALSE)