Nothing
## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>"
)
R <- function() knitr::include_graphics("Rlogo.png", dpi = 5000)
## ----echo=FALSE---------------------------------------------------------------
knitr::include_graphics("hex-hce.png", dpi = 1000)
## ----eval=TRUE----------------------------------------------------------------
library(hce)
packageVersion("hce")
## ----eval=FALSE, include=FALSE------------------------------------------------
# devtools::load_all()
## -----------------------------------------------------------------------------
ls("package:hce")
## -----------------------------------------------------------------------------
args("hce")
## -----------------------------------------------------------------------------
set.seed(2022)
n <- 100
dat <- hce(GROUP = rep(x = c("I", "II", "III"), each = 100),
TRTP = sample(x = c("Active", "Control"), size = n*3, replace = TRUE))
class(dat)
## -----------------------------------------------------------------------------
head(dat)
## -----------------------------------------------------------------------------
unique(dat[, c("GROUP", "GROUPN")])
## -----------------------------------------------------------------------------
set.seed(2022)
n <- 100
dat <- hce(GROUP = rep(x = c("I", "II", "III"), each = 100),
TRTP = sample(x = c("A", "P"), size = n*3, replace = TRUE), ORD = c("III", "II", "I"))
unique(dat[, c("GROUP", "GROUPN")])
## -----------------------------------------------------------------------------
data(HCE1)
unique(HCE1$GROUP)
## -----------------------------------------------------------------------------
HCE <- hce(GROUP = HCE1$GROUP, TRTP = HCE1$TRTP, AVAL0 = HCE1$AVAL0,
ORD = c("TTE1", "TTE2", "TTE3", "TTE4", "C"))
class(HCE)
head(HCE)
## -----------------------------------------------------------------------------
data(HCE1, package = "hce")
class(HCE1)
head(HCE1)
## -----------------------------------------------------------------------------
dat1 <- as_hce(HCE2)
str(dat1)
## -----------------------------------------------------------------------------
args("simHCE")
## -----------------------------------------------------------------------------
Rates_A <- c(1.72, 1.74, 0.58, 1.5, 1)
Rates_P <- c(2.47, 2.24, 2.9, 4, 6)
dat3 <- simHCE(n = 2500, n0 = 1500, TTE_A = Rates_A,
TTE_P = Rates_P,
CM_A = -3, CM_P = -6,
CSD_A = 16, CSD_P = 15,
fixedfy = 3, seed = 2023)
## -----------------------------------------------------------------------------
class(dat3)
head(dat3)
## -----------------------------------------------------------------------------
methods(class = "hce")
## -----------------------------------------------------------------------------
HCE <- hce(GROUP = HCE3$GROUP, TRTP = HCE3$TRTP,
ORD = c("TTE1", "TTE2", "TTE3", "TTE4", "C"))
calcWO(HCE)
summaryWO(HCE)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.